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Der induzierte Seitenwind von Flugzeugen.* 
Von W, Jacobs und E. Truckenbrodt. 


1, Einleitung. In dem Bestreben, die Flugeigenschaften eines Flugzeuges rechnerisch zu be- 
herrschen, ist man schon seit langerer Zeit bemiiht, die hierfiir erforderlichen aerodynamischen 
Beiwerte, die sog. Stabilitatsbeiwerte, nicht nur aus Windkanalmessungen zu beschaffen, sondern 
méglichst auch theoretisch zu ermitteln1*. Hierbei spielt die gegenseitige Beeinflussung der 
Flugzeugteile, die sog. Interferenz, eine ausschlaggebende Rolle 3451415, Fir die Stabilitat um die 
Querachse (Langsstabilitat) liegt eine besonders starke Interferenz zwischen Fliigel, Rumpf und 
Héhenleitwerk vor, die vom Fliigelabwind herrithrt. Durch diesen Abwind wird in bekannter 

_ Weise der Beitrag des Héhenleitwerks zum Moment um die Querachse stark abgemindert. Dieses 
Problem ist theoretisch und experimentell weitgehend untersucht worden 78%, Neuerdings ist 
nun eine ganz ahnliche Beeinflussung des Seitenleitwerks und seines Beitrages zur Richtungs- 
stabilitat durch den Fliigel und Rumpf aufgedeckt worden *1, Je nach der Lage des Fliigels 
relativ zum Rumpf (Tief-, Mittel-, Hochdecker) wird beim Schieben des Flugzeuges am Ort des 
Seitenleitwerks ein induzierter Seitenwind erzeugt, der den Beitrag des Seitenleitwerks zum Gier- 
moment und damit zur Richtungsstabilitat vergréBert bzw. verringert. Eine physikalische Deu- 
tung dieses Einflusses wird im folgenden Abschnitt 2 gegeben. Im Abschnitt 3 wird erstmalig 
versucht, diese recht verwickelten Strémungsvorgange in stark vereinfachter Form einer rech- 
nerischen Behandlung zuganglich zu machen. Experimentelle Untersuchungen zu diesem Pro- 
blemkreis liegen bereits vor 4 13, Im Abschnitt 4 wird dann tiber umfangreiche Windkanal- 
messungen an verschiedenen Fliigel-Rumpf-Anordnungen beim Schieben berichtet. Die Mes- 
sungen werden mit den theoretischen Rechnungen verglichen. 


2. Physikalische Grundlagen. Eine qualitative Erklarung des induzierten Seitenwindes wurde 
erstmalig von H. Multhopp ? und gleichzeitig auch von H. Schlichting und W. Frenz!° angegeben. 
Die vom Schieben herriithrende Querumstrémung des Rumpfes erzeugt am Fliigel eine anti- 
metrische Auftriebsverteilung, welche auch die Ursache fiir das durch den Rumpf stark beein- 
fluBte Schieberollmoment ist. Entscheidend fiir das Vorzeichen dieser zusadtzlichen Auftriebs- 


* Die vorliegende Untersuchung wurde in den Jahren 1942 bis 1944 am damaligen Aerodynamischen 
Institut der Technischen Hochschule Braunschweig (Leitung: Prof. Dr. H. Schlichting) ausgefiihrt. Sie ist 
unter dem Titel ,,Theoretische und experimentelle Untersuchungen tiber den induzierten Seitenwind von 
Flugzeugen“‘ im Jahre 1944 von der Fakultat Maschinenwesen der Technischen Hochschule Braunschweig 
als Habilitationsschrift des erstgenannten Verfassers angenommen worden. (Referent Prof. Dr. H. Schlich- 
ting, Korreferent Prof. Dr. H. Blenk.) Die hier vorliegende gekiirzte Fassung wurde vom zweitgenannten 
Verfasser angefertigt. Von ihm stammt auch die vereinfachte Neufassung des theoretischen Teiles (Ab- 
schnitt 3). 

1 G, Mathias, Die Seitenstabilitat des ungesteuerten Normalfluges und ihre technischen Vorbedingungen. 
DVL-Jabrb. 1932, S. VIII 17. 

2 C. H. Zimmermann, An analysis of lateral stability in power-off flight with charts for use in design, 
NACA Rep. No. 589, 1937. 

3 H. Multhopp, Luftfahrt-Forsch. 18 (1941), S. 51. 

4 F. Riegels und W. Liess, Jahrbuch 1942 der Deutschen Luftfahrtforschung, S. I 366. 

5 H. Schlichting, Jahrbuch 1940 der Deutschen Luftfahrtforschung, S. I 113. 

6 I. Fliigge-Lotz und D. Kiichemann, Jahrbuch 1938 der Deutschen Luftfahrtforschung, 5. I 172. 

7 H. Multhopp, Luftfabrt-Forsch. 15 (1938), S. 463. 

8 E. Truckenbrodt, Ing.-Arch. 18 (1950), S. 233. 

9 H. Trienes und E. Truckenbrodt, Ing.-Arch. 20 (1952), S. 26. 

10 Hf, Schlichting und W. Frenz, Jahrbuch 1941 der Deutschen Luftfahrtforschung, 5. I 300. 

11 R,O. House und A. R. Wallace, Wind-tunnel investigation of the effect of interference on lateral 
stability-characteristics of four NACA 23012 wings, an elliptical und a circular fuselage, and vertical fins, 
NACA Rep. No. 705, 1941. 

12 M. J. Bamber und R. 0. House, Wind-tunnel investigation of effect of yaw on lateral stability cha- 
racteristics, NACA Tech. Note No. 703 und No. 730, 1939. 

138 H. Schlichting und W. Frenz, Jahrbuch 1943 der Deutschen Luftfahrtforschung. Vorabdruck in Techn. 
Berichte, Bd. 11 (1944), H. 6. 

14 W. Jacobs, Lift and moment changes due to the fuselage for a yawed aeroplane with unswept 
and swept wings. Flygtekniska Férséksanstalten, Stockholm, Medd. 34, 1950. 

15 W. Jacobs, Untersuchungen iiber die Interferenz zwischen Seitenleitwerk und Deltafliigel beim 
Schieben. Flygtekniska Férséksanstalten, Stockholm, Institutsbericht, 1953. 
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verteilung ist die relative Lage des Fliigels zum Rumpf. Ist z. B. der Fliigel in einer Schulter- 


decker-Anordnung mit dem Rumpf verbunden, so gibt es am schiebenden Fliigel an der vor- | 


gehenden Fliigelhalfte zusitzliche Aufwartsgeschwindigkeiten und an der zuriickliegenden 
Fliigelhalfte entsprechende Abwartsgeschwindigkeiten (Abb.1). Bei einer Tiefdecker-Anordnung 


liegen die Verhaltnisse gerade umgekehrt. Praktisch beschrankt sich die durch den Rumpf 


hervorgerufene Auftriebsanderung auf einen kleinen Be- 


starke Spitzen, so da® in diesem Gebiet am Fliigel ein 
starker Druckgradient in Spannweitenrichtung entsteht, 
der am Fliigel groBe induzierte Seitengeschwindigkeiten 
erzeugt. Da die am Fliigel herrschenden induzierten 
Seitenwinde sich naturgemaB hinter dem Fligel fort- 
setzen, treffen sie auf das Seitenleitwerk. GemaS ihrer 
Entstehung sind diese induzierten Seitengeschwindig- 
keiten proportional dem Schiebewinkel / und unab- 
hangig vom Anstellwinkel «. Unterhalb und oberhalb der 
von der Fliigelhinterkante abgehenden Wirbelflache haben 
sie verschiedenes Vorzeichen. Je nach der Lage der Wirbel- 
schicht zum Seitenleitwerk wird also der EinfluB mehr 
oder weniger stark sein. Nimmt man an, dai das Seiten- 
leitwerk oberhalb der Wirbelflache liegt, so ergibt die 
induzierte Seitenstrémung beim Tiefdecker eine Erhéhung 
der Richtungsstabilitat durch das Seitenleitwerk und beim 
Schulterdecker eine Verminderung. 

Durch die Anstellung des Rumpfes werden aufBer den 
eben behandelten antimetrischen Verteilungen noch sym- 
metrische Auftriebsinderungen lings Spannweite erzeugt, 
mit denen wir uns aber nur insofern befassen wollen, als 
sie zur Bestimmung der Lage der Wirbelschicht von Be- 
deutung sind. 

Die bei einem ungepfeilten Fliigel allein (mit und 
ohne V-Form) infolge Schiebens entstehende Auftriebs- 
verteilung langs Spannweite (Abb.2) weist gegeniiber dem 
Rumpfeinflu8 nur sehr geringe antimetrische Betrage 
auf und tragt infolgedessen zum induzierten Seitenwind 
nichts Wesentliches bei. 

Ist der induzierte Seitenwind v, bekannt, so berechnet 
sich der Anstrémungswinkel des Seitenleitwerkes nach 
Abb.1 zu 


os = B—Ods, (1) 


Abb. 1. Erlauterungsskizze zur Entstehung der ‘ . : : 
induzierten Seitengeschwindigkeit. is wobei p der Schiebewinkel und 


I. Fliigel allein; II, Fliigel mit Rumpf. v 
os = — e (2) 
der induzierte Seitenwindwinkel ist. Aus (1) und (2) folgt sofort 
00 « 06 
sees Gg Seg ye eee dacs, (3) 
ap ap V Op 
Fiir das Giermoment des ganzen Flugzeuges gilt 
TsFs qs 
Cn = €No sFq °S° (4) 


Hierin ist cy = N/q F's der Giermomentenbeiwert des ganzen Flugzeuges, und cy, stellt den 
Anteil des Flugzeuges ohne Seitenleitwerk dar; rg ist der Abstand des Seitenleitwerks von der 
Momentenbezugsachse, F's die Seitenleitwerksflache, gs der Staudruck am Ort des Seitenleit- 
werkes und c,s der auf die Seitenleitwerksflache bezogene Auftriebsbeiwert; s = b/2 ist die halbe 
Fliigelspannweite, F die Fliigelflache und q der Anstrémungsstaudruck (vgl. Abb. 1). 

Fiir den Beiwert der Richtungsstabilitat folgt aus (4): 


dcn _ _9€No ssc (aes 
Of nr OB a ea Add sos? (5) 


reich in Rumpfnahe. Sie erreicht dort allerdings ziemlich 
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wobei 
his OF ESS O65 Le 
boop (Ges) 1 ae ty ae (6) 


ein dem Wirkungsfaktor des Héhenleitwerkes entsprechend gebildeter Wirkungsfaktor des Seiten- 
leitwerkes ist!. Wir nehmen dabei naherungsweise an, daB keine Abschirmung des Seitenleit- 
werkes erfolgt; setzen also qs/q ~ 1. 

_Aus dem bereits Gesagten kénnen wir bei der iiblichen Lage des Seitenleitwerkes fiir den 
Wirkungsfaktor des Seitenleitwerkes folgendes aussagen: 


Tiefdecker : Co oaT (stabilisierend), 
Schulterdecker: C5 < 1 (instabilisierend), 
Mitteldecker: ¢s = 1 (neutral). 


Gegenstand der vorliegenden Untersuchung ist es nun, die Wirkungsfaktoren einmal theoretisch 
zu berechnen und zum anderen experimentell zu bestimmen. 


A Caan, 


cael 


\\ 


40 —08 —O6 —O4 -2 7] G2 g4 G6 G8 10 
voraere Fligelhaltre hintere Flugelnaitre 


Abb. 2. Auftriebsverteilung eines schiebenden Flugzeuges. Rechteckfliigel: A = 5, ee 270 Ce = 1,05 


I. Fligel allein, 8 = 0°. II. Fliigel allein, 8 = 20°. Ill. Fligel allein mit 3° V-Stellung, 6 = 20°. 
IV. Fligel + Rumpf (Schulterdecker), B = 20°. 


3. Die theoretische Berechnung des Wirkungsfaktors. a) Die Formel zur Berechnung 
des brtlichen Wirkungsfaktors. Ein theoretischer Ansatz zur Abschatzung des Wirkungs- 
faktors des Leitwerkes ist bereits in der zitierten Arbeit von H. Schlichting und W. Frenz? ent- 
halten. Beschrankt man sich auf Punkte, die sich unmittelbar ober- oder unterhalb der Wirbel- 
schicht befinden, so findet man fiir die induzierte Seitengeschwindigkeit in bekannter Weise 


ee (7) 


Hierin bedeutet ['(y) die Zirkulationsverteilung langs Spannweite. Das obere positive Vorzeichen 
gilt fiir Punkte oberhalb der Wirbelschicht, das untere negative Vorzeichen entsprechend fir 
Punkte unterhalb der Wirbelschicht. Fiir die Berechnung des induzierten Seitenwindes am Ort 
des Seitenleitwerkes kommt es also auf den Gradienten der Zirkulationsverteilung d/‘/dy in der 
Umgebung der Fliigelmitte an. Wie wir aus Abb. 2 ersehen, gibt besonders der Rumpf den iiber- 
wiegenden Beitrag zum Gradienten d/'/dy und damit zum induzierten Seitenwind. 
Fir die Zirkulation fiihren wir den dimensionslosen Wert 

a (8) 
ein. Die Zirkulation kénnen wir aufteilen in einen Anteil, der von der Geradanblasung herriihrt, 
y,, und in einen solchen, der infolge Schiebens zusatzlich entsteht, By,. 


Setzen wir 
Y¥=YetPYe (9) 


a /4 00 iA: 
1 Fiir das Héhenleitwerk ist der Wirkungsfaktor durch(y = a Ga) a == = definiert. 


/ 


2 Siehe FuBnote 10 von S, 1. 
iE 
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und 


Uf) = 


ies y) 
$s 


Te ere 8 10) ¥ 
abe (19) 0 
dann kénnen wir fiir (6) auch schreiben 
dy, \| 

a mae if! 
Cs aa 1 + dy . ( ) i 


Wie bereits gesagt wurde, gilt diese Beziehung fiir Punkte in Hohe der Wirbelschicht. Da die — 
induzierten Geschwindigkeiten in der Wirbelschicht im allgemeinen eine starke Spitze aufweisen, 
ist zu erwarten, da diese Abschatzung zur Berechnung des Wirkungsfaktors des Leitwerkes auBer- 
halb der Wirbelschicht nicht ausreicht, sondern stets zu groBe Werte liefert. Man muf bei der 
genaueren Berechnung auch auf die Verhiltnisse in gréBeren senkrechten Abstanden von der 
Wirbelflache eingehen. Es kommt also darauf an, das gesamte Strémungsfeld am Ort des Seiten- 
leitwerkes zu bestimmen. Fiir unsere theoretischen Uberlegungen nehmen wir an, daB der Fliigel 
ungepfeilt ist. 

Nach der Prandilschen Tragfliigeltheorie ersetzen wir den Tragfliigel durch die gebundenen 
Wirbel am Ort des Fliigels (1/4-Linie) und die freien Wirbel, die sich bei Vernachlassigung des 


a symmetrische Verfeilung b antimetrische Verreilung 
b=25 


/ 
/ 
/ 
| 
/ 
fy 
} 
CLIYULL / 
/ Ge | 


Geschwintigkelts- —»|<— 
Komponenten Dy: ~=\=—-— 


Abb. 3. Das Wirbelsystem des schiebenden Abb. 4. Zum Einflu8 des schraggestellten freien Wirbelbandes auf die induzierte 
liigels. Seitengeschwindigkeit. 


Aufrollvorganges in Anstrémungsrichtung geradlinig nach hinten bis ins Unendliche erstrecken 
(Abb.3). Durch die Schragstellung der freien Wirbel hinter dem Fliigel gestaltet sich die mathe- 
matische Behandlung des Problems etwas umstiandlich. Durch eine qualitative Abschatzung 
dieses Einflusses kann man aber zeigen, da die Schragstellung bei unseren Betrachtungen nur 
untergeordnete Bedeutung hat, so da wir sie fiir unsere Rechnungen vernachlassigen kénnen. 

Unsere Uberlegungen zur Abschatzung des induzierten Seitenwindes durch die Schragstellung 
der freien Wirbel fithren wir so durch, daB wir bei vorgegebener Auftriebsverteilung einmal die 
freien Wirbel nach links und zum anderen nach rechts nach hinten abgehen lassen und uns je- 
weils die Richtung und Gréfe der induzierten Seitengeschwindigkeit qualitativ tiberlegen (Abb.4). 

Die Auftriebsverteilung spalten wir in einen symmetrischen und einen antimetrischen Anteil 
auf, die wir jeweils getrennt voneinander behandeln. Das Wirbelsystem, das zu der symmetrischen 
Verteilung gehért, ersetzen wir niherungsweise durch einen schraggestellten Hufeisenwirbel der 
Spannweitenerstreckung b = 2 s (Abb. 4a). Das entsprechende Wirbelsystem der antimetrischen 
Verteilung denken wir uns durch zwei schraggestellte Hufeisenwirbel ersetzt, deren Breite je- 
weils nur die Halbspannweite b/2 = s ist (Abb. 4b). Auf der einen Seite ist der gebundene Wirbel 
positiv und auf der anderen Seite negativ drehend. Entsprechend ist auch der Drehsinn der freien 
Wirbel. Insbesondere geht in der Fligelmitte ein starker Witbel ab, der aus der Summe der 
beiden gebundenen Wirbel gebildet wird. Wie aus Abb.4 hervorgeht, liefern die gebundenen 
Wirbel keinen Beitrag zum induzierten Seitenwind. Das Seitenleitwerk (Aufpunkt) mége sich 
soweit hinter dem Fliigel befinden, da® die Lingen der freien Wirbel, gleichgiiltig wie ihre 
Richtung hinter dem Fliigel ist, naherungsweise gleich grof sind. 


r 
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Betrachten wir zunichst die symmetrische Anordnung. Gehen die freien Wirbel nach rechts 
ab, dann erzeugt der linke Wirbel eine gréBere Seitengeschwindigkeit als der rechte Wirbel, so 
dafs eine Geschwindigkeitskomponente nach rechts vorhanden ist, wenn wir annehmen, daf sich 
die Lage des Leitwerkes oberhalb der Wirbelflache befindet. Gehen dagegen die freien Wirbel 
nach links ab, dann verbleibt eine Geschwindigkeitskomponente nach links (Abb. 4a). Daraus 
folgt, daB bei symmetrischer Auftriebsverteilung der Einflu8 der Schragstellung der freien Wirbel 
auf den induzierten Seitenwind proportional dem Winkel der Schragstellung (Schiebewinkel) ist. 
Wie wir nun aber aus (7) ersehen, ist die induzierte Seitengeschwindigkeit proportional dem 
Gradienten dI"/dy. Dieser ist aber bei symmetrischer Auftriebsverteilung in Fliigelmitte nur ge- 
ring, so daf} auch der Einflu8, der durch die Schragstellung der freien Wirbel entsteht, klein sein 
wird und bei unseren weiteren Uberlegungen vernachlassigt werden kann. 

Anders liegen die Verhaltnisse bei der antimetrischen Anordnung. Gehen die freien Wirbel 


nach rechts ab, so erzeugen der linke und rechte Wirbel nach links gerichtete Geschwindigkeits- 


komponenten. Der starke in Fliigelmitte abgehende Wirbel erzeugt aber infolge seines wesent- 


lich kleineren Abstandes zum Leitwerk eine sehr viel gréBere Seitengeschwindigkeit nach rechts, 


so daf insgesamt eine Komponente nach rechts tibrig bleibt. Gehen die freien Wirbel nach links 
ab, so ergibt sich bei ungeanderter Auftriebsverteilung aus entsprechenden Uberlegungen die- 
selbe nach rechts gerichtete Seitengeschwindigkeit (Abb. 4b). Daraus folgt nun das wichtige 
Ergebnis, daB bei antimetrischer Auftriebsverteilung die Schragstellung der 
freien Wirbel die GréBe und Richtung der induzierten Seitengeschwindigkeit 
in erster Naherung nicht beeinfluSt. Wir kénnen also naherungsweise die Schrag- 
stellung des freien Wirbelfeldes bei der durch den Rumpf erzeugten antimetrischen Auftriebs- 
verteilung vernachlassigen. Diese Tatsache bringt natiirlich wesentliche Vereinfachungen fir 
die Berechnung des induzierten Geschwindigkeitsfeldes hinter dem Fliigel 1. 

Wir gehen noch einen Schritt weiter in der Vereinfachung, indem wir von der Tatsache Ge- 
brauch machen, da® sich das Seitenleitwerk verhaltnismaBig weit hinter dem Fliigel befindet, 
so dafS wir im mathematischen Sinne x — — oo gehen lassen kénnen. Die Anwendung des Biot- 
Savartschen Gesetzes liefert dann 


s eee ; ) 
DO aC reece rer | 

Fe (12) 
ato x=aa | POVig=ytap: | 


Diese Ausdriicke gelten unter der Annahme, daf die Wirbelschichtin der x, y-Ebeneliegt, so daB z der 
Abstand des Aufpunktes von der Witbelschichtist. Das hier benutzte Koordinatensystem (x, y, 2) ist 
identisch mit dem in Abb. 14angegebenen experimentellen Achsensystem (%,, ¥.,%)- Will man die 
Absenkung der Wirbelschicht beriicksichtigen, deren Lage wir mit 2, = z,(x, y) bezeichnen wollen, so 
ist in (12) z durch (z —z,) zu ersetzen. Fithren wir jetzt noch die dimensionslosen Koordinaten 


Raw EEE 2 aad a pele 13 
g et: le i= ign? ‘3 om C1 ae ( ) 
sowie nach (8) y = /'/bV ein, dann wird fiir die induzierte Seitengeschwindigkeit aus (12) 
1 
2s ee ee a en) n! a dn’ . 14. 
iyakes Raley) Sire ie ) Tq aye oF al 


ail 

1 Die induzierten Geschwindigkeitskomponenten fiir das in Abb. 3 dargestellte Wirbelsystem wurden 
vom erstgenannten Verfasser auf Grund des Biot-Savartschen Gesetzes berechnet. Die Ausdriicke, die man 
dabei findet, sind aber sehr uniibersichtlich, so daB auf die Wiedergabe hier verzichtet wird, zumal wir ge- 
zeigt haben, da die Schragstellung der freien Wirbel keinen entscheidenden EinfluS hat. Greift man auf 
eine neuere Arbeit vom zweitgenannten Verfasser (Das Geschwindigkeitspotential der tragenden Flache 
bei inkompressibler Strémung (erscheint demniachst in der Z. angew. Math. Mech.) zurtick, so kann man 
den allgemeinen Ausdruck fiir das Geschwindigkeitspotential folgendermafen finden: 


ddueti rm “i 1 x cos B —(y —y’) sinB dy’ 
2 ead ais eonheaG a yeahhe |. -\etooyrte |} 


Hieraus berechnen sich die Geschwindigkeitskomponenten zu 
vx = OD/0x., vy = 0@/dy , vz = OD/0z, 
was fir B =0 und x = — o sofort auf (12) fihrt. 
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Ist jetzt die Auftriebsverteilung langs Spannweite bekannt, so liefert die Auswertung des | 
Integrals (14) den gesuchten induzierten Seitenwind in Abhangigkeit von der Hoch- und Seiten- }- 
lage des Aufpunktes (Seitenleitwerkes). Im Fall, daB ¢ = ¢, ist, d. h. in Héhe der Wirbelschicht, ] 
ist das Integral iiber eine singulare Stelle zu bilden. Eine genaue Grenzwertbildung liefert dann 
den Wert der Gleichung (7), der sich in dimensionsloser Foim folgenderma8en schreibt: 


55 (—c0, 9, 0) = J (a). (15) J 


Fir die Zirkulationsverteilung hatten wir nach (9) bereits geschrieben y = y, + Py,. Hierin be- ij 
deutet y, die Verteilung bei Geradanblasung. Das Glied Pyp erfaBt die zusatzlich beim Schieben | 
erzeugte Verteilung. Sie umfaBt die folgenden beiden Anteile: | 

1. Die Zirkulationsinderung des Fliigels allein, herrithrend von der Herabsetzung der wirk- | 
samen Anstrémungsgeschwindigkeit, von der Schragstellung des freien Wirbelbandes und von 
der seitlichen Umstrémung der Fligelenden. 


2. Die Zirkulationsinderung herrithrend von der seitlichen Rumpfumstrémung. | 

Den ersten Anteil, den man fiir einen ungepfeilten Fliigel z. B. nach Weissinger + berechnen 
kann (vgl. Abb.2), wollen wir, wie bereits erwahnt wurde, vernachlassigen. Auch der Einflu8 
einer V-Stellung, der ebenfalls in Abb. 2 gezeigt ist, kann unberiicksichtigt bleiben. Der zweite 
Anteil, der RumpfeinfluB8, ist dagegen besonders stark, so da® wir mit Py, fortan nur die zusatz- | 
liche Auftriebsverteilung durch den Rumpf bezeichnen wollen. / 


Setzen wir jetzt (9) in (14) ein, fithren die Differentiation nach dem Schiebewinkel f aus und | 
setzen dieses Ei gebnis in (6) ein, dann finden wir den értlichen Wirkungsfaktor in einem bestimm- | 
ten Punkt des Seitenleitwerkes zu | 


1 
2K 1) (Neg 


Osan (i Sabai Lae (Cee eiGE ad (16) | 
ey 
Den Grenzwert ¢ = ¢, hatten wir bereits in (11) zu | 
= fody 
Csi (0) = 1 = Spe (16a) 


bestimmt. 


Die Auswertung der Formel (16) setzt nun die Kenntnis der durch den Rumpf erzeugten Auf- 
triebsverteilung voraus. AuBerdem muf die Lage der Wirbelschicht €, bekannt sein. In einer 
friiheren Arbeit wurden vom erstgenannten Verfasser? bereits die beim Schieben zusatzlich auf- 
tretenden Auftriebsverteilungen berechnet. Dabei wurden Fliigelform, Fliigelhochlage, Seitenver- 
haltnis, Rumpfform und Verhaltnis von Fligelbreite zu Rumpfbreite geandert. Der Rumpf- 
einflu8 auf den Abwind hinter dem Fliigel macht sich dadurch bemerkbar, daf eine antimetrische | 
Verformung der Wirbelflache eintritt. Da wir der Einfachheit halber mit einer ebenen Wirbel- | 
flache rechnen wollen, wird sich deren mittlere Héhenlage durch diesen EinfluB8 nicht andern. 

Zu dem unsymmetrischen Kinflu8B des Rumpfes kommt noch ein symmetrischer hinzu, der 
davon herriithrt, daB der Rumpf bei Anstellung eine Anstrémungskomponente von unten er- | 
fahrt. Dadurch ergeben sich beidseitig des Rumpfes zusitzliche Anstellwinkelanderungen, die dann | 
eine entsprechende Auftriebsanderung zur Folge haben. Dieser Einflu8 wurde in einer Arbeit | 
von I’. Riegels und W. Liess* fiir verschiedene Rumpfformen berechnet. Durch diese sym- | 
metrischen Verteilungen, die itberdies stark von der Hochlage des Fliigels zum Rumpf abhingen, | 
wird die Lage der Wirbelschicht am Ort des Leitwerkes entscheidend mitbestimmt. 


b) Die numerische Berechnung des Wirkungsfaktors fiir einigeRumpf-Fligel- | 
Anordnungen. Wir haben zunichst das Feld der induzierten Geschwindigkeiten fiir einen sym- 
metrisch angestrémten Rechteckfliigel vom Seitenverhaltnis 1 = 5 berechnet. In Abb. 5 
haben wir die Kurven dy/a = konst. bzw. 6s/a = konst. einmal in der (n,€—€,)-Ebene und | 
zum anderen die Werte dy/a bzw. ds/a iiber 7 fiir verschiedene Hochlagen (¢—C,) aufgetragen. — 
Aus den Geschwindigkeitskomponenten v, und v, wurde auch das Stromlinienbild der induzicrten 
Stiémurg ermittelt. Der Fall des schiecbenden Fligels wurde nicht durchgerechnet. Aus dem | 
Vergleich der Auftriebsverteilungen des schiebenden und nicht schiebenden Fliigels, vel. Abb.2, | 


J. Weissinger, Jahrbuch 1940 der Deutschen Luftfahrtforschung, S. I 138. 
* W. Jacobs, Jahrbuch 1941 der Deutschen Luftfahrtforschung, S. I 170. 
3 Siehe FuBnote 4 von S. 1. 
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kann man schlieBen, da durch das Schieben lediglich in den uns hier nicht interessierenden 


Randgebieten Abweichungen auftreten werden. Diese Untersuchungen am Fliigel allein dienten 


lediglich der Orientierung. 
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Wie wir schon mehrfach betont haben, liefert die beim Schieben durch den Rumpf erzeugte 
antimetrische Auftriebsverteilung den entscheidenden Einflu8 bei der Berechnung des induzierten 
Seitenwindes. Wir greifen jetzt eine spezielle Fliigel-Rumpf-Anordnung heraus, und zwar einen 
Schulterdecker, bestehend aus einem Rechteckfliigel vom Seitenverhaltnis A =5 und einem 
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Rumpf von kreisférmigem Querschnitt mit der Fligelhochlage z/R = —1(R = Halbmesser des 
Rumpfquerschnittes). Das Verhaltnis von Fliigelspannweite zu Rumpfdurchmesser betragt 
b/D = 7,5. Wir betrachten jetzt nur den zusitzlich durch die Queranstrémung des Rumpfes er- 
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Abb. 6, Zusitzliche Auftriebsverteilung und zusitzliches Stromungsfeld eines Schulterdeckers beim Schieben. &é =—oo, 
Fligel: Rechteck, A = 5, one = 2a. Rumpf: Kreisquerschnitt, b/D = 7,5, a/R =—1l. 


zeugten KinfluB. Abb.6 stellt die der Abb.5 entsprechende Auftragung dar. Fiir die uns besonders 
interessierenden induzierten Seitengeschwindigkeiten entnchmen wir der Abbildung folgendes: 
Die Kurven 65/8 = konst. verlaufen zur (x, z)-Ebene symmetrisch. Oberhalb der Wirbelschicht, 
(C—¢,) <0, ergeben sich im Bereich der Flugzeugsymmetrieebene negative Werte, untethalb, 
(C—C,)>0, entsprechend positive Werte. Innerhalb der Wirbelschicht sind die Werte am g16Bten,. 
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Aus der Auftragung ist der starke EinfluB des Strémungsfeldes des Fliigels auf die Lage des 
Seitenleitwerkes ersichtlich, wobei jedoch auch noch die Lage der Wirbelschicht C, bekannt 
sein muf. Ausgehend von (6) haben wir die értliche Seitenleitwerks-Wirksamkeit 


Csortl = tt —( 


3g ) 
e. ortl 


gebildet und in Abb.7 dieKurven €¢ 5, = konst. fiir den Schulterdecker aufgetragen. Fiir ein 
zentrales Seitenleitwerk iiblicher Hochlage tritt z. B. eine betrachtliche Abminderung der Seiten- 
leitwerks-Wirksamkeit ein. Kurz iiber der Wirbelschicht ist sie sogar null. Die ungiinstigste Lage 
des Seitenleitwerkes ist im vorliegenden Fall die, wo das Leitwerk gerade auf der Wirbelschicht 
aufsitzt. Mit wachsendem Abstand von der Wirbelschicht, z. B. beim Hohersetzen des Leit- 
werkes, wird der Rumpfeinflu8 geringer ; die Leitwerkswirksamkeit nahert sich dem Wert 1. Ver- 
schieben wir das Leitwerk aus der Mitte heraus, gehen wir also zum doppelkieligen Leitwerk iiber, 
so steigt die Wirksamkeit des Leitwerkes an und itherschreitet sogar den Wert 1. Dieses bedeutet, 
da8 der .Rumpf beim Schieben auf das Seitenleitwerk stabilisierend wirkt. Befindet sich das 
Leitwerk unterhalb der Wirbelschicht, so kehren sich die Verhaltnisse um. 
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Abb.7. Kurven konstanten theoretischen 6rtlichen Wirkungsfaktors 
des Seitenleitwerkes (eg grt] = konst.) hinter einem 


Schulterdecker, &=—o, 


Fligel: Rechteck, A = 5, ee 
Rumpf: Kreisquerschnitt, 6/D = 7,5 


ee 
z,/R = —1. 
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Abb. 8. Kurven konstanten theoretischen értlichen Wirkungs- 
faktors des Seitenleitwerkes (CS 6rtl = konst.) hinter einem 


Schulterdecker. &=—wo. 
Fligel: Rechteck, A = 5, oe = 2a. 
Rumpf; Kreisquerschnitt, b/D = 7,5, z,/R = — 0,5. 


Bei einem Tiefdecker sind die Verhaltnisse gerade umgekehrt wie bei einem Schulterdecker, 
wenn man die Lage des Seitenleitwerkes jeweils von der Wirbelschicht aus mift. Da aber fir 
die Lage der Wirbelschicht, die ja im wesentlichen vom Anstellwinkel des Fliigels abhangt, die 
Umkehrbarkeit zwischen Schulter- und Tiefdecker nicht gilt, kann man die Wirkungsfaktoren 
eines Tiefdeckers nicht ohne weiteres entgegengesetzt gleich denen eines entsprechenden Schulter- 
deckers setzen. Bei gleichem Fligel und Rumpf haben wir noch fiir die Hochlage z/R = —0,5 
die ortlichen Wirkungsfaktoren am Ort des Seitenleitwerkes berechnet (Abb.8). Die Hinfliisse 
sind kleiner als die bei z/R —=-—1. In einem anderen Fall (Abb.9) haben wir noch die Rumpf- 
breite bei sonst gleichen Verhaltnissen wie in Abb. 7 (a/R = —1) gedndert; und zwar haben 
wir gewahlt b/D —5. Mit gréBerer Rumpfdicke nimmt auch die Anderung des Wirkungsfaktors zu. 

Die Lage des Seiten- und Héhenleitwerkes im Seitenwind- bzw. Abwindfeld ist in Abb.10, 
11 und 12 fiir einen Schulter-, Mittel- und Tiefdecker bei einkieligem Seitenleitwerk dargestellt. 
Fiir den Mitteldecker ist nur das Abwindfeld gezeichnet worden, da in diesem Fall die induzierte 
Seitengeschwindigkeit verschwindet. Die Rechnung wurde durchgefihrt fiir einen Rumpf mit 
kreisférmigem Querschnitt und einem Rechteckfliigel mit symmetrischem Profil. 

Rechtecks 1h 5s cz, =. 2, ep = 0°; 
Kreisquerschnitt, b/D = 7,5, 
Héhenleitwerk: by/b = 0,33, Seitenleitwerk: bs/b = 0,175 ; 


Die Daten sind: Fliigel : 
Rumpf: 
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er bedeutet die Schrankung des Fliigels zur Rumpfachse, bs die Héhe des Seitenleitwerkes und | 
bj, die Spannweite des Héhenleitwerkes. Die Abhangigkeit vom c,-Wert wurde durch die 
Wahl von c, = 0; 0,5 und 1,0 untersucht. Das Koordinatensystem ist flugzeugfest. Der Null- | 
punkt liegt auf der Rumpfachse in 1/4 der Fliigeltiefe. Die Kurven Cy —konst. und Cs = konst. 
sind gerechnet worden fiir £ = —oco, wihrend die Absenkung der Wirhbelflache sich auf die 
Lage des Leitwerkes bei & = —1,2 bezieht. Die theoretische Lage der Wirbelschicht bei ver- 
schiedenen Schrinkungen des Fliigels zur Rumpfachse ist fir den Schulter- und Tiefdecker 
in Abb. 13 dargestellt. Fiir ein bestimmtes Leitwerk ist tiber seine Hohe bzw. Breite zu inte- 
grieren, um die mittlere Wirksamkeit zu erhalten. Fir den Wirkungsfaktor des Seitenleitwerkes 
te wird der folgende Mittelwert gebildet: 
~b7 


ate Biss erin ace { (5), dz. (11) 
10 O95, O98 10 op bs y OB Jort 
S 
4 bone ba Fiir das gewahlte Flugzeug sind diese 
és Mittel-Werte in den Abbildungen mit- 
105 98 105 angegeben. Eine etwas genauere Formel 
08 fiir den mittleren Wirkungsfaktor wiirde 
wt 95 41 sich dadurch ergeben, das man die Tiefen- | 
12 y te verteilung des Seitenleitwerkes berick- 
0 sichtigt. An Stelle von (17) mite man 
98 08 dann schreiben 
09 NX 15 09 Le ; 1: (00s a ” 
0 aes ft 03 3 Fs |G y a: ee 
gu Y In der nachfolgenden Arbeit des erst- 
genannten Verfassers! wird ein bequemes 
10 105 105, 70 Rechenverfahren fiir eine serienmabige Be- | 
Re ps mal 7 rechnung des Wirkungsfaktors des Seiten- 
52 Be Ca uy ae ye 9° leitwerkes angegeben. 
Abb. 9. Kurven konstanten theoretischen 6rtlichen Wirkungsfaktors 
des Seitenleitwerkes (fg g-4j = konst.) hinter einem 4. Die experimentelle Bestimmung des 
Pee ae oe ge Es Wirkungsfaktors. a) Die MeBypethoden 


Fliigel: Rechteck, A = 5, ie, —tovare i _ 
Rumpf: Kreisquerschnitt, b/D = 5, 2,/R = —1. und die Auswertung. Zur Bestimmung 


des Wirkungsfaktors eines Seitenleitwerkes 
aus Windkanalmessungen fiir eine vorgegebene Fliigel-Rumpf-Anordnung stehen folgende zwei 
Methoden zur Verfiigung: 

1. Sechskomponentenmessung mit und ohne Leitwerk. Aus einer Messung des 
Giermomentes fiir das Gesamtmodell (Fliigel + Rumpf + Leitwerk) wnd fiir das Modell ohne 
Leitwerk (Fliigel -- Rumpf) erhalt man durch Differenzbildung den Beitrag des Seitenleitwerkes 
zum Giermoment. Aus dem Verhiltnis dieser Differenz zu dem Weit des Leitwerks allein er- 
gibt sich dann der Wirkungsfaktor des Seitenleitwerkes. Unter Vernachlassigung einer Staudruck- 
anderung am Leitwerk stellt die Abweichung des Wirkungsfaktors vom Wert 1 den Einflu8 des 
induzierten Seitenwindes dar. Wir schreiben 

» (Cx )riRyL— (Cy) FLR 

Conia (ez : (19) 

ew )r 

Der Index N hei Cs bedeutet, daB die Giermomentenmessungen zur Bestimmung von Cs benutzt 
worden sind. cy = dcy/df gibt den Giermomentenanstieg mit dem Schiebewinkel an. Um den 
Beitrag des Leitwerks allein zum Giermomentenanstieg zu bestimmen, nimmt man zweck- 
maBigerweise die Differenz aus einer Messung von Rumpf + Leitwerk und Rumpf allein, wo- 
durch man in erster Naherung die Abschirmung durch den Rumpf miterfaBt hat: 


(ew )n = (Cy )rix— (ew)a - (20) 
Diese Beziehung setzen wir in (19) ein und erhalten 


(Cy) ry RyL —(¢y)ryR 


Csy = (CAN or : (21) 


1 W. Jacobs, Ing.-Arch. 21 (1953) S. 23. 
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Analog hierzu kénnen wir den Wirkungsfaktor auch aus den Seitenkraftsmessungen finden: 


(c,) F4R4L— (rir 
(Cry 7 = (ie 


aa (2) 
Der Index y bei Cs deutet die Herleitung aus der Seitenkraft an. Es bedeutet cy = dc,/06 den 


kurven £5 =const. 
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Abb. 10. Lage des Seiten- und Héhenleitwerkes im Seiten- und Abwindfeld eines Schulterdeckers 
‘ ? fiir die drei c,-Werte 0; 0,5 und 1,0. 
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} 


Seitenkraftsanstieg mit dem Schiebewinkel. Die Werte fiir Cs, und Cs, miissen in erster Naherung 


gleich sein. 

2. Punktweises Ausmessen des Stroémungsfeldes durch Staudruck- und Rich- 
tungsmessungen. Diese Methode ist sehr mithsam. Sie liefert natiirlich einen besseren Ver- 
gleich mit der Theorie, da die Werte in den einzelnen Punkten direkt verglichen werden kénnen, 
wahrend bei der Kraftmessung nur Integralwerte tiber das ganze Leitwerk erhalten werden. 
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Abb, 11. Lage des Seiten- und Héhenleitwerkes im Seiten- und Abwindfeld eines Tiefdeckers 
fiir die drei c,-Werte 0; 0,5 und 1,0. 
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Von uns wurden beide Methoden an einem 
Schulter-, Mittel- und Tiefdecker durchgefiihrt. 
Es wurde die Abhangigkeit des Wirkungsfaktors 
von der Hochlage des Fliigels zam Rumpf und von 
der Lage des Seitenleitwerkes untersucht. Auch 
die Einstellung (Schrankung) des Fliigels zum 
Rumpf wurde verandert. 

b) Beschreibung der Modelle. Fir die 
Messungen wurden die in Abb.15 dargestellten Mo- 
delle benutzt. Der Fliigel war ein gerader Recht- 
eckfliigel ohne V-Stellung und ohne Verwindung 
vom Profil NACA 23012. Er besa8 Endkappen, 
deren Abrundungshalbmesser gleich der értlichen 
Profildicke war. Der Fliigel war in Metall-Gips-Bau- 
weise hergestellt. Die Abmessungen des Fliigels sind: 


kurven Cy=const. 


Spannweite: De OO mats 
Fliigelflache : P= 0.1146 m2, 
Fliigeltiefe : i015 0m. 


Seitenverhaltnis: A = 5,18. 


Der Rumpf ist aus dem in fritheren Messungen 
verwendeten Rotationsellipsoid vom Achsenverhalt- 
nis 1: 7 entstanden, indem der hintere Teil eine 
Zuspitzung erhielt, um méglichst glatten AbfluB 
und damit geringe Stérungen hinter dem Rumpf zu 
haben. Er war aus vergiitetem Holz hergestellt. 
Durch auswechselbare LEinsatzstiicke mit ent- 
sprechenden Fliigeldurchbriichen konnten die ge- 
wiinschten Hochlagen des Fligels erzielt werden. 
Das Rumpfende war mit einer auswechselbaren 


Metallspitze versehen, so daB eine Haltestange fiir 
das in der Hohe verschiebbare Seitenleitwerk nach | G65 
der auf Abb.15 ersichtlichen Art eingesetzt werden ee Be t— ga 
konnte. Die Haltestange war gekrépft, um auch i 05: 
gréBere Tieflagen des Leitwerkes zu erreichen. it 
Die Abmessungen des Rumpfes sind: ——- Ga aN 
Lange: a = 0,650 m , _ gy 
GréBter Durchmesser: D = 0,104 m. 2 Fi Nae ; oe Vbe 
Als Leitwerk wurde ein einkieliges Seitenleit- x G65) Fliige/ 
werk ohne Hékenleitwerk vom Profil NACA 0010 070 


benutzt. Es war in der Héhe verschieblich an- Se AP OY 


gebracht. Die Abmessungen des Seitenleitwerkes Abb, 12. Lage des Héhenleitwerkes im Abwindfeld eines 


sind: Mitteldeckers fiir die drei c,-Werte 0; 0,5 und 1,0. 
Schulterdecker lietdecker 
. o 
&=0° a: Ep=0 


Abb, 13. Theoretische Lage der Wirbelschicht bei verschiedenen Schrankungen des Fligels zum Rumpf fiir Schulter- 
und _ Tiefdecker. 
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Hohe (Spannweite): bs = 0,113 m , Flache: Fs = 0,00918 m?, 

Seitenverhaltnis: As = 1,40, Hebelarm: rs = 0,50m. 

Als Fligel-Rumpf-Anordnungen wurden Schulter-, Mittel- und Tiefdecker gewahlt. 

Die Hoch- bzw. Tieflagen fiir Schulter- und Tiefdecker wurden so bestimmt, da Ober- bzw. 

Unterseite des Fliigels in der Mitte gerade mit der Rumpfoberfléche abschlossen. Als charakte- 

ristische Gr6éBe wird der Ab-| 

is oe stand z, der 1/4-Linie der Fliigel-_ 
| 


M ve ebene von der Rumpfachse an- 
gesehen. Es ergaben sich fol- 
gende Hochlagen: 


be Schulterdecker : 

Sap ers 2 z/D! = —0,154.4 
Mitteldecker : = + 0,050, 

= Tiefdecker: = + 0,854. 


D’ bedeutet den Durchmesser 
des Rumpfes am Ort der 1/4- 
Linie des Fliigels. 

Die Riicklage des Fliigels 
(e = Abstand der 1/4-Linie von | 


“Soe sa OM 


dimungsebene : 


Ne x, y, 2 flugzeugfestes Achsensystem, 
Xo» Yeo #, experimentelles Achsensystem. 


1 der Rumpfnase) war in allen | 
‘ | Fallen gleich und betrug 
Abb. 14. Achsensysteme. ela = 0,35 ; 


Die Schrankungen er zwischen Fliigel und Rumpf (positiv im Sinne erhéhten Anstellwinkels) 
waren 
Schulterdecker: ep = 0°, 3°, 6°, 
Mitteldecker: Ole 
Tiefdecker : 0° Ges 
Der Fligel war mit einer 
Fliigelwurzelverkleidung 
versehen, die in Abb. 15 
y angedeutet ist. Neben 
WRN “ey scene h or aaa: it 
(+) Mitteldecker mae ce 
A — Tetecer und ohne Seitenleitwerk 
wurde auch noch der 
Rumpf allein mit und 
ohne Seitenleitwerk ge- 
messen. 
Zur Richtungsmessung 
diente ein Vierfingergerat 
Abb. 15. Modellibersicht. iitblicher Bauart. 
Mabe jn wm c) Durchfithrung 
und Umfang der 
Messungen. Die Mes- 
sungen wurden durch- 
gefiihrt im 1,2-m-Wind- 
kanal des Aerodynami- | 
schen Institutes der Technischen Hochschule Braunschweig bei einer Windgeschwindigkeit von 
40 m/sec. Dieses entspricht auf die Fliigeltiefe bezogen einer Reynoldsschen Zahl Re = 4,2 - 10°. 
1. Sechskomponentenmessungen. Fir die Sechskomponentenmessungen wurde das 
experimentelle Achsensystem, Abb. 14, verwendet. Der Bezugspunkt lag bei den Fliigel-Rumpf- 
Anordnungen in der Mitte der Fliigelvorderkante, bei den Messungen am Rumpf allein wurde 
er auf die Rumpfachse mit derselben Riicklage verlegt. Sadmtliche Messungen wurden mit 
und ohne Seitenleitwerk durchgefiihrt. Im letzten Fall wurde nur das Leitwerk abgenommen; 
die Haltestange also am Modell belassen, um ausschlieSlich den Beitrag des Seitenleitwerkes zur 
Richtungsstabilitat ermitteln zu kénnen, Fir den Schulterdecker und fiir den Tiefdecker wurden | 
jeweils drei Hochlagen des Seitenleitwerkes gewahlt. Als charakteristische GréBe wird hierbei | 
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der Abstand des tiefsten Leitwerkspunktes von der Rumpfachse angesehen. Wir bezeichnen 
diesen Abstand mit 2s. Auferdem wurde noch die Schrankung des Fliigels geandert. Die Drei- 
komponentenmessungen wurden fiir eine bestimmte Fliigel-Rumpf-Anordnung nur einmal durch- 
gefiihrt, da die Hochlage des Seitenleitwerkes kaum von Einflu8 sein diufte. Die Sechskompo- 
nentenmessungen wurden fiir den Schiebewinkelbereich f = +30° im allgemeinen fiir drei An- 
‘stellwinkel durchgefithrt. Nur fiir den Mitteldecker beschrankten wir uns auf einen Anstell- 
winkel und auf eine Héhenlage des Leitwerkes, weil hier die erwarteten Hinfliisse sehr klein sind. 
Bei den Messungen am Rumpf allein mit und ohne Seitenleitwerk zeigte sich, daB die Abschir- 
mung des Seitenleitwerkes durch den Rumpf ziemlich stark vom Anstellwinkel abhingt. Es muB- 
ten daher hierfiir auch dieselben Anstellwinkel gewahlt werden wie bei den Gesamtmodellen. 
Kine Ubersicht der durchgefiihiten Sechskomponentenmessungen ist in Tabelle 1 enthalten. 
2. Richtungsmessungen. 

Bei den Richtungsmessungen Tabelle 1. Ubersicht der Drei- und Sechskomponentenmessungen. 
merden das Leitwerk sowie die) © 


} 


Haltestange entfernt, und der Anordnung | Leitwerk Schrinkung ep oe me ae 

Meee sider pansenden spitze) 2 
ee : | 

versehen, um méglichst geringe sehulterdecker . .| einkielig | 0°, 3°, 6° | —15, +18, +42 

Stérungen infolge AbreiBens der Tiefdecker ; 5 Ne. 6 —15, + 5,-+25 

Strémung zu haben. Es wurde Mitteldecker .. . * [20° +5 


jeweils das gesamte Richtungs- 
feld, also sowohl der Abwind wie auch der Seitenwind gemessen. Die Messung geschah so, daB in 
einem bestimmten Abstand bei einem festen Schiebewinkel eine Ebene x = konst. ausgemessen 
wurde und zwar in verschiedenen Schichten z = konst. Es wurden die Schiebewinkel f = 0°, 
5°, 10° und 15° gewahlt. Der Koordinatenursprung wurde hierbei stets in das Rumpfende gelegt. 

Das bei der Messung benutzte Vierfingergeraét wurde vorher im ungestérten Strahl geeicht. 
Mit geniigender Genauigkeit ergab sich fiir den Seitenwind dieselbe Eichkurve wie fiir den Ab- 
wind. Diese Methode war in unserem Fall einfacher als die Nullmethode, wobei das Gerat so- 
lange gedreht wird, bis die zwei gegeniiberliegenden Pitot-Réhrchen den gleichen Druck anzeigen. 
Bei Messungen mit Vierfingergeraten, deren gegeneinander geschaltete Réhrchen einen gewissen 
Abstand voneinander haben, ist besonders darauf zu achten, daB in Gebieten mit starken 
Geschwindigkeitsgradienten die Messungen gefalscht werden, weil die Réhrchen dann mit ver- 
schiedenen Geschwindigkeiten angeblasen werden. Dieses tritt ein beim Durchgang durch die 
Nachlaufdelle von Rumpf und Fligel. Da diese Dellen bei gesunder Strémung sehr schmal sind, 
kann man die etwas auferhalb der Delle gemessenen Werte in die Delle hinein extrapolieren. 
Dieses wurde bei den Auswertungen beachtet. 

Eine Ubersicht der durchgefithiten Richtungsmessungen geht aus Tabelle 2 hervor. 


Neben den Messungen dicht hinter ge 
Tabelle 2. Ubersicht der Richtungsmessungen. 


dem Rumpfende (x = — 10mm) wurden 

auch einige Stichmessungen fiir den Sehriinkung | a | Eapendiiiater 

Me8Bquerschnitt bei «x = —120 mm or: €R a 2 ee | (ny 

durchgefiihrt. <a 
d) MeBergebnisse und Ver- Schulterdecker | Oe Ge — 1°, 5,7° —10 

gleich mit der Theorie. Als Bezugs- Tiefdecker 0° Seb astyl © —10 


gréBen sind bei der Bildung der dimen- 
sionslosen Beiwerte die GréBen des Fliigels benutzt worden. Es gilt also fiir die Seitenkraft 
= Y/Fq, fiir das Rollmoment cy, = L/F'sq und fiir das Giermoment cy = N/F sq. 
1. Dreikomponentenmessungen. Abb.16enthialt die Polaren und Momentenkurven der 
verschiedenen Anordnungen. Die Widerstande der verschiedenen Modelle unterscheiden sich bis 
auf den Bereich hoherer c,-Werte fast gar nicht voneinander. Trotz guter Fliigelwurzelverklei- 
dung reiBt beim Tiefdecker die Strémung eher ab als beim Mittel- und Schulterdecker. Dieses ist 
aber fiir unsere Messungen belanglos, da wir in dem gesunden Bereich der anliegenden Strémung 
bleiben. 
2.Sechskomponentenmessungen. Auf die gesamten Messungen wollen wir im einzelnen 

nicht eingehen. Es wurde jeweils fiir eine bestimmte Fligel-Rumpf-Anordnung eine Messung mit 
und ohne Seitenleitwerk gemacht. Hinsichtlich des Auftriebes, des Widerstandes und des Kipp- 
momentes in Abhangigkeit vom Schiebewinkel wurden die schon aus friiheren Messungen ! 


Cy 


1 Vgl. FuBnoten 5, 10 und 11 von S.1 sowie E. Moller, Jahrbuch 1942 der Deutschen Luftfahrt- 
forschung, S. I 336. 
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bekannten Interferenzen beobachtet. Es spielt hierbei die Hochlage des Fliigels zum Rumpf eine) 
wichtige Rolle. 
Cu Gy 
a 


06 0 Schulterdecker Ep=0° 
es ” 6 2 


° Mitteldecker 0° 
° Miefdecker os 


G2 


? 
f 
B | | | na 0 O20 
& 


Viefdecker. 


Ac z 
er 
ee on ee 
Ofe La 
5s) Abb.17. Rollmomenten-Anstieg c= be, [a8 
Z Zps= -15 TM. in Abhangigkeit von « fiir die Modelle 
008. ZZ ae =+5mm mit und ohne Seitenleitwerk. 
fog 
SE —-— = 6 mm 
4 SS Dinh 
| “guy — = +42 mm 
——-- =ohneSlW 


a sen 
ao = ee act geen eaaeege Oo - 1A 


Das Schieberoll moment einer Fliigel-Rumpf-Anordnung ist stark abhangig von der Hoch- 
lage des Fliigels zum Rumpf (Abb.17). Beim Schulterdecker ist dieses zusitzliche Rollmoment 
positiv, beim Tiefdecker negativ und beim Mitteldecker null. Es 1a4Bt sich gut vorausberechnen, 


ae 
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wie frithere Arbeiten zeigen 1. Beim Gesamtmodell mit Seitenleitwerk wird das Rollmoment 


natiirlich noch durch das Seitenleitwerk beeinflu®t. Es kommt hierbei 


auf die Hochlage des 


0) Hefdecker -——--— =ohneSlLW 
a 
% Schulter- x 


/ 
hey 
*. me] Oo]. 7 T T = 
eget —o— Zyg=— mM 
=a = Po 
=| O40 0 ——_|_ _- = + 4mm 
a | Fef- —4+— =+é5mm 
‘4 decker —r»— =+42mm 


decker SS 


¥ = ) Schutter - 
020. See atl ia decker 


ells 
6 & 70° 10° 
Abb. 18. Seitenkraft-Anstieg <p = @Cy/0p 
in Abhingigkeit von « fiir die Modelle 
Fi ca aan mit und ohne Seitenleitwerk. 
Ye 
ie 
| EL. —— 
2 se 
jl meee ff 
70 aaa ia 
& 
a 
Sse. pon ae ae ee 
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= 
—o— Zy,=— 15mm 


ae UY 
—— =+8mm 
—i—- |S s = + mM 
—1—- ==+ 42mm 


= ohne SLW 


Abb. 19. Giermomenten-Anstieg cK = 0C n/p 


bchulter- in Abhangigkeit von « fir die Modelle 
echer mit und ohne Seitenleitwerk. 
atl AZ 
10° 


Seitenleitwerkes an, so da® sich die Einfliisse hier ttberdecken. Mit wachesnder Hochlage des 
Seitenleitwerkes waichst das Rollmoment ebenfalls an, was aus Abb.17 zu entnehmen ist. 


1 Vel. FuBnote 2 von S.6 sowie G. Braun und H. Scharn, Jahrbuch 1942 der Deutschen Luftfahrt- 


forschung, S. I 246. 


2 
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Die Schiebe-Seitenkraft und das Schiebe-Giermo ment wollen wir genauer betrach- 
ten, da hierauf der induzierte Seitenwind von groBem EinfluB ist. In Abb.18 sind die Stabilitats- | 
beiwerte der Seitenkraft (c;) und in Abb.19 die des Giermomentes (cn) fiir verschiedene Schran-— 
kungen und verschiedene Hochlagen des Seitenleitwerkes ttber dem Anstellwinkel aufgetragen. 
AuBerdem sind die Ergebnisse ohne Seitenleitwerk mit darin enthalten. Die Stabilitatsbeiwerte 
sind auf das flugzeugfeste Koordinatensystem umgerechnet worden. Der Koordinatenursprung 
liegt auf der Rumpfachse und hat die Riicklage der Fliigelvorderkante. Nach unseren fritheren 
Uberlegungen folgt, daB die Schiebeseitenkraft des Tiefdeckers iiber der des Schulterdeckers liegt, 
und daB eine Abhiangigkeit von der Hochlage des Seitenleitwerkes und von der Schrankung 
zwischen Fliigel und Rumpf besteht. Entsprechendes gilt fiir das Schiebegiermoment. Unsere 
Uberlegungen werden durch die Messungen bestitigt. Beim Tiefdecker wirkt der induzierte 
Seitenwind stabilisierend und beim Schulterdecker instabilisierend; beim Mitteldecker ist kein 
besonderer EinfluB festzustellen. s a 

In Abb. 20 und 21 sind die Wirkungsfaktoren Cs und €s, nach den Beziehungen (21) 
und (22) iiber der Hochlage des Seitenleitwerkes (293) und dem Anstellwinkel « fiir die ver- 
schiedenen Fliigel-Rumpf-Anordnungen aufgetragen. Betrachten wir zunachst den Wirkungs- 


™ 
Mtleldecker = 


= -_——e poe 


Schulter- 
deck 


ea eas Zos 
cee Mii, Dele RA lias 2 ete MES od -20 -0 a a a 2} 
Abb. 20. Wirkungsfaktoren des Seitenleitwerkes lon Abb. 21. Wirkungsfaktoren des Seitenleitwerkes sy 
gebildet aus dem Giermoment in Abhingigkeit von der gebildet aus der Seitenkraft in Abhangigkeit von der 
Hochlage des Seitenleitwerkes. Hochlage des Seitenleitwerkes. 


faktor Csy in Abhangigkeit von der Hochlage 2s (Abb. 20). Bei der Messung wurden drei Hoch- | 
lagen des Seitenleitwerkes gewahlt. Diese Hochlagen wurden durch Vorversuche festgelegt, um 
den maximalen EinfluB des induzierten Seitenwindes zu erfassen. Wie die Kurven zeigen, tritt 


dieses Maximum auch deutlich hervor. Der Schulterdecker zeigt durchweg [sy-Werte kleiner als 
1, wie auch theoretisch zu erwarten ist. Fiir « = —1° ergeben sich bei einer Fligelschrankung 
én = 3° und 6° und einer Hochlage des Seitenleitwerkes von zs = +20 mm sogar Wirkungs- 
faktoren um 0,4, Mit wachsender Hoch- bzw. Tieflage steigen sie wieder an. Weiter ist bei allen 
Schulterdecker-Anordnungen ein Ansteigen der Wirkungsfaktoren mit dem Anstellwinkel zu 
beobachten. Dieses ist auf folgende Weise zu erklaren: Nehmen wir an, daB bei « = —1° die 
Wirbelflache das Leitwerk unten berithrt bzw. etwas schneidet, so wird mit wachsendem Anstell- 
winkel die Wirbelschicht in das Leitwerk hineinwandern. Da zu beiden Seiten der Wirbelschicht 
entgegengesetzte Seitengeschwindigkeiten herrschen, so heben sich beim Eintauchen des Leit- 
werkes in die Wirbelschicht die Einfliisse z. T. auf; die Wirkungsfaktoren miissen also ansteigen. 
Der Wert 1 wiirde erreicht werden, wenn das Leitwerk durch die Wirbelschicht genau in der 
Mitte getroffen wiirde. Hierbei ist ein Leitwerk gleicher Tiefe und gleicher c, 5,y-Verteilung 
vorausgesetzt. 

Beim Tiefdecker mite bei entsprechender Lage des Leitwerkes zur Wirbelschicht genau das 
Entgegengesetzte festzustellen sein, und zwar miBte mit wachsendem Anstellwinkel der Wir- 
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kungsfaktor von einem Wert gréBer als 1 an abnahmen. Diese Abnahme ist bei den Messungen 
nicht aufgetreten. Dieses rithrt wohl daher, daB die Wisbelschicht bei « = —1° noch etwas Hee 
halb des Leitwerkes lag, so daB mit wachsendem Anstellwinkel zuniichst noch eine Zunahme ein- 
treten muf. Erst mit gréBeren Anstellwinkeln mite dann die Abnahme erfolgen. Von der 
Schrankung ist der Wirkungsfaktor natiirlich aus den gleichen Griinden ebenfalls abhangig. Als 


-maximaler Wert tsy ergibt sich beim Tiefdecker mit eg — 0° und « = 5,7° der Wert 1,56. Dieser 


° 
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Abb. 22. Kurven gleichen Abwindes (5;, = konst.) und gleichen Seitenwindes (6, = konst.) aus Richtungsmessungen 
hinter dem Schulterdecker bei verschiedenen Schiebewinkeln. ER 0°, « 1°, « = —10 mm. 


Wert liegt also ungefahr soviel tiber 1 wie der kleinste Wert des Schulterdeckers unter 1. Fiir 
den Mitteldecker liefert die Messung fast genau 1, wie nach der Theorie zu erwarten war. Es 
wurde der Mitteldecker mit der Schrankung ¢g = 0° nur bei einem Anstellwinkel vermessen, da 
gréBere Einfliisse vermutlich nicht auftreten werden. Ganz ahnliches Verhalten wie die Wirkungs- 
faktoren ‘acta zeigen die aus der Seitenkrafi abgeleiteten Werte Csy, Abb. 21. 

3. Richtungsmessungen. Im AnschluB an die Sechskomponentenmessungen wurden um- 
fangreiche Richtungsmessungen durchgefiihrt, deren Umfang aus Tabelle 2 hervorgeht. Es wurde 
dabei das gesamte Abwind- und Seitenwindfeld ausgemessen. Es wurden jeweils fiir die Schiebe- 

9* 
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winkel 6 = 5°, 10° und 15° die Kurven konstanten Abwindes (6) und konstanten Seitenwindes 


(ds) in dem flugzeugfesten Koordinatensystem bestimmt (Abb. 14). 


Wir verzichten darauf, hier samtliche MeBergebnisse wiederzugeben. Als Beispiel einer voll- 
standigen Messung haben wir in Abb.22 die Verhaltnisse fiir den Schulterdecker mit der Schran- | 
kung ep = 0° bei dem Anstellwinkel « = —1° dargestellt. Fiir die anderen Anordnungen sehen | 


die Kurven ahnlich aus. Die Kurven gleichen Ab- bzw. Aufwindes sind um 90° gegeniiber den 
Kurven gleichen Seitenwindes gedreht. 


Aus den Richtungsmessungen am Ort des Leitwerkes haben wir den 6rtlichen Wirkungs- | 
faktor nach (6) bestimmt. Die Ergebnisse fiir samtliche Anordnungen sind in Abb.23, 24 und 25 | 


aufgetragen. Die Kurven C35,, = konst. sind etwas aus der Flugzeugsymmetrieebene heraus 
aie versetzt, was wahrscheinlich die Folge 
100) davon ist, daB wir nur positive Schiebe- 
winkel gemessen haben und zur Be- 
stimmung von 0063/06 den Mittelwert 


5p| Se ~ 


~150 =700 
die kleiner als 1 sind, unterhalb solche, 


theoretisch mit dem Anstellwinkel zu 
erwartende Aufwartsbewegung des 
ganzen Geschwindigkeitsfeldes nur so 
gering ist, riihrt wohl daher, daf die 
Wirbelflache hinter dem Fligel stark 
vom Rumpf gefiihrt wird. Durch die 
Schrankung des Fliigels gegeniiber dem 
Rumpf um ér = 6°, Abb.24, erfahrt 
das gesamte Geschwindigkeitsfeld eine 
Absenkung nach unten. Fir den Tief- 
decker, Abb.25, kehren sich die Ver- 


etwa um. Im allgemeinen kann man 
feststellen, daB die im Abschnitt 3 theo- 
retisch abgeleiteten Beziehungen durch 
die Messungen gut bestatigt werden. 

Um den Wirkungsfaktor des ge- 
samten Seitenleitwerkes zu erhalten, 


700 5 : ° : 
Abb. 23. Ortlicher Wirkungsfaktor des Seitenleitwerkes ist nach (17) ein. entsprechender Mittel- 


aus Richtungsmessungen hinter dem Schulterdecker. ER=9°, x =—10mm. wert zu bilden. Der Vergleich der auf 


“% =—I1° = 0) und « = 5,7° es 0,5). p : : 
Oe) (Ca ) diese Weise berechneten Wirkungs- 


faktoren mit denen der Sechskomponentenmessung ist in Abb.26 durchgefithrt worden. Die Uber- 
einstimmung ist einigermafien gut. Sie wiirde besser werden, wenn man entsprechend (18) noch 
die Tiefenverteilung des Seitenleitwerkes beriicksichtigen wiirde. 


In Abb.27 haben wir noch aus den Ergebnissen der theoretischen Rechnungen des Ab- 
schnittes 3 die Wirkungsfaktoren fiir die gemessenen Modelle ermittelt. Die dargestellten Kurven 
sind mit denen aus Abb.20, 21 und 26 zu vergleichen. Die Ubereinstimmung ist wieder einiger- 
maBen gut. Nach den Messungen sind die durch den induzierten Seitenwind hervorgerufenen 
Stabilitatsanderungen starker als nach der Theorie. Die Abhingigkeit vom Anstellwinkel kommt 
nach den Messungen wesentlich schwacher heraus als nach der Theorie. 


5. Zusammenfassung. Bei einem schiebenden Flugzeug wird durch die Verdrangungsstrémung 
des Rumpfes eine antimetrische Auftriebsverteilung am Fliigel erzeugt, die insbesondere fiir den 
Schulter- und Tiefdecker recht betrachtlich ist. Als Folge des sich hieraus am Tragfliigel aus- 
bildenden Druckgradienten in Querrichtung entsteht ein induzierter Seitenwind, der analog dem 
Abwind am Héhenleitwerk die Wirksamkeit des Seitenleitwerkes und damit die Richtungs- 


zwischen B = 0° und f= 10° genommen | 
haben. Fiir den Schulterdecker, Abb.23, | 
ergeben sich in Ubereinstimmung mit | 
y den theoretischen Uberlegungen ober- | 
mm 70"  halb der Wirbelschicht C5 5,,-Werte, 


die gréBer als 1 sind. Der EinfluB des | 
Anstellwinkels ist nur gering. Daf die _ 


haltnisse gegenittber dem Schulterdecker | 
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ao = —I1° (cg + 0) und a = 5,7° (cg =~ 0,5). 
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Abb, 25. Ortlicher Wirkungsfaktor des Seitenleitwerkes aus Richtungsmessungen hinter dem Tiefdecker. 
€ 
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stabilitat stark beeinflussen kann. Im AnschluB an theoretische Untersuchungen wird tiber um-) 
fangreiche Sechskomponenten- und Richtungsmessungen an Modellen eines Schulter-, Mittel-» 
und Tiefdeckers zur Bestimmung des Wirkungsfaktors des Seitenleitwerkes berichtet. Unter-) 
sucht wurden die Abhangigkeit des Wirkungsfaktors des Seitenleitwerkes vom Anstellwinkel, 
von der Schrankung des Fligels zur Rumpfachse und von der Hochlage des Seitenleit-) 
werkes. Bei den Richtungsmessungen wurde sowohl der Abwind wie auch der Seitenwind ge-) 


Tief- 
Tieflecker decker 
y eo 
Ca 
Vs | Schulfer- Schulter- 
decker decker 


= (i Sechskomp—Messg. 

6° ———-Richtungsmessung 

Zos 
“20-10 0 0 20 30 40 50mm C0 i/o 0 0 20 JO YO 50 tM 
Abb. 26. Wirkungsfaktoren des Seitenleitwerkes te Abb. 27. Theoretische Wirkungsfaktoren des Seitenleitwerkes 


fiir Schulter- und Tiefdecker in Abhangigkeit von der 


Vergleich zwischen Sechskomponentenmessung (sn) nea Ra etein: ig 
ochlage des Seitenleitwerkes. 


und Richtungsmessung. ~« = —1°. 


messen. Kurven gleichen Abwindes bzw. gleichen Seitenwindes geben ein anschauliches Bild 
der gesamten Zusatzstrémung. Die Messungen werden mit theoretischen Rechnungen ver- 
glichen. Es zeigt sich hierbei, da} der nach der Theorie erwartete EinfluB in den Messungen noch 
starker hervortritt. Bei normaler Lage des Seitenleitwerkes bewirkt der zusatzliche induzierte 
Seitenwind beim Schulterdecker eine Erniedrigung (€¢s< 1), beim Tiefdecker eine Erhéhung 
(Cs >1) und beim Mitteldecker keine Anderung der Richtungsstabilitat. 


(Eingegangen am 25, April 1952.) 


Anschriften der Verfasser: Dr. habil. W. Jacobs, Stockholm-Bromma (Schweden), Stora Mossens Backe 6. 


Dr.-Ing. E. Truckenbrodt, Braunschweig, Technische Hochschule, 
Institut fiir Strémungsmechanik, 
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Ein vereinfachtes Verfahren 
zur Berechnung des induzierten Seitenwindes von Flugzeugen *. 


Von W. Jacobs. 


1. Einleitung. In der vorstehenden Arbeit! haben wir ausfiihrlich iiber theoretische und experi- 
mentelle Untersuchungen iiber den induzierten Seitenwind von Flugzeugen berichtet. Hierbei 
ergab sich in gewissen Fallen ein betrachtlicher Einflu8® des induzierten Seitenwindes auf die 
Wirksamkeit des Seitenleitwerkes, so daB er fiir die Richtungsstabilitat des Flugzeuges er- 
hebliche Bedeutung hat. Der Vergleich der Theorie mit Messungen hatte recht befriedigende 
Ubereinstimmung gezeigt. Da jedoch die theoretische Ermittlung des Feldes der induzierten 
Seitenwindgeschwindigkeit in dem ganzen Leitwerksbereich bisher ziemlich umstandlich ist — 
es mu fiir jeden einzelnen Punkt des Feldes ein Integral ausgewertet werden — erscheint 
es wiinschenswert, das Rechenverfahren fiir die Ermittlung des Feldes des induzierten Seiten- 
windes erheblich zu vereinfachen. 

Im folgenden soll deshalb ein Verfahren angegeben werden, das es gestattet, auf einfache 
Weise das Geschwindigkeitsfeld des induzierten Seitenwindes in der ganzen Querebene weit 
hinter dem Fliigel (x = — oo) zu berechnen. Wie man aus der Theorie des Abwindes weiS, 
ist das Geschwindigkeitsfeld am Ort des Leitwerkes (x ~ — 6/2) nicht wesentlich verschieden 
von demjenigen sehr weit hinter dem Fliigel. Dasselbe gilt auch fiir den induzierten Seitenwind, 
so da die nachfolgende Theorie, die fiir « = — oo durchgefiihrt wird, die Verhaltnisse am Ort 
des Leitwerkes gut wiedergeben diirfte. 


2. Bezeichnungen und BezugsgréBen. Das benutzte Achsensystem (x, y, 2) ist identisch mit 
dem ,,experimentellen Achsensystem:: (x,; y = y,3 2,). Der Index e ist wegen der bequemeren 
Schreibweise fortgelassen worden. Ferner bedeuten: 


b=2s die Fliigelspannweite, 

I die Fliigeltiefe, 

F die Fliigelflache, 

Eas =, = 2 Ta = mit der Halbspannweite dimensionslos gemachte Koordinaten, 
¢, die Hochlage der Wirbelschicht, 

z die Hochlage des Fliigels zum Rumpf, 

a den Anstellwinkel, 

B den Schiebewinkel, 

V die Fluggeschwindigkeit, 


Vz, Vy, Vv, die Komponenten der induzierten Geschwindigkeit, 


v 5 ° . . 0 
6s = — =e den induzierten Seitenwindwinkel, 


Cs =1— a den Wirkungsfaktor des Seitenleitwerkes, 
I'(y) die Zirkulationsverteilung, 

y(n) = aa die dimensionslose Zirkulationsverteilung, 
K,,...K, charakteristische Konstanten, 

J ee = das Seitenverhaltnis des Fliigels, 

D den Rumpfdurchmesser. 


* Die vorliegenden Untersuchungen wurden zum Teil bereits 1944 ausgefiihrt und in dem heute 
nicht mehr zuginglichen Bericht 45/5 des friheren Aerodynamischen Institutes der Technischen Hoch- 
schule Braunschweig (Leitung: Prof. Dr. H. Schlichting) mitgeteilt. In zwei Berichten des Flugtechnischen 
Institutes Stockholm-Ulvsunda 2»? wurden die Untersuchungen fortgefihrt. 

1 W. Jacobs u. E. Truckenbrodt, Ing.-Arch. 21 (1953) S. 1. ; - 

2 W. Jacobs, The influence of the induced sidewind on the efficiency of the vertical tail; a simplified 
method for calculation. Flygtekniska forséksanstalten, Medd. 35 (1950). ; N 

3 W. Jacobs, The induced sidewind behind swept wings at subsonic velocities. Flygtekniska forsdks- 
anstalten, Medd. 41 (1951). 
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3. Verfahren zur vereinfachten Berechnung des induzierten Seitenwindes. Auf die physi ie 
kalischen Grundlagen und die Entstehung des induzierten Seitenwindes soll hier nicht nahen)> 
eingegangen werden. Dies ist ausfithrlich in der erwahnten friiheren Arbeit geschehen. Es magx}y 
nur erwahnt werden, daB fiir gewisse Hochlagen des Fliigels zum Rumpf beim Schieben solchey 

zusatzlichen antisymmetrischen Auftriebsverteilungen erzeugt werden, die ihrerseits einen i 
betrachtlichen Seitenwind induzieren. Dieser induzierte Seitenwind ist entgegengesetzt gleich i inti / 
spiegelbildlich zur Wirbelflache gelegenen Punkten. Da das induzierte Feld an die Wirbelflachet) 
gebunden ist, kommt es bei der Beurteilung der Richtungsstabilitat auf die Lage des Seitenleit-}) 
werkes relativ zur Wirbelflache stark an. Der induzierte Seitenwind ist proportional dem Schiebe-})) 


winkel f, er ist unabhingig vom Anstellwinkel a. 


Fiir den Einflu® des induzierten Seitenwindes auf das Seitenleitwerk wurde analog zum| r 
Abwindfaktor des Héhenleitwerkes ein Wirkungsfaktor des Seitenleitwerkes eingefiihrt von)): 
der Grobe! 1) 


Q 


euity ae 


Cs = 1— 


Ferner wurde gezeigt, daB fiir groBe Abstande hinter dem ine der EinfluB der Schragstellung s 
des Wirbelbandes vernachlassigt werden kann. Wir rechnen also mit der Auftriebsverteilung \ 
des schiebenden Fliigels, nehmen aber ein zur Symmetrieebene paralleles Wirbelsystem an. fl 
Durch Zerlegung der Zirkulation in einen bei Geradanstrémung vorhandenen Anteil y,(7) und 
einen vom Schieben herrithrenden zusitzlichen Anteil? 8 yg(7) in der Form 


v(n)= y(n) +B ve(n) » (2) f 


hatte sich fiir den értlichen Wirkungsfaktor Oe 


' py 1 a2 =< : =" r | : 
son (yn. 6 —G) = 1+ i [( —n? + (E—G)E ya(n')dn’ . (3) | 


Die Auswertung dieses Integrales ist in jedem Einzelfall nur durch miihsame numerische 
oder graphische Integration méglich. Um diese miihsame Auswertung zu vermeiden, wollen 
wir die antisymmetrische Auftriebsverteilung yg(yj) durch eine Summe von Polynomen dar- 
stellen und fiir die letzteren die Auswertung des Integrales (3) ein fiir allemal ausfiihren. Die 
meisten praktisch vorkommenden antisymmetrischen Auftriebsverteilungen lassen sich durch 
die folgende Polynom-Summe darstellen: 


volt) = SK fil) (—1<n< +1). (4 


Die K, sind freie Konstanten, die fiir jede Auftriebsverteilung passend zu wahlen sind. Die 
Polynome f,(j) miissen in der Mitte und am Fligelende verschwinden, f,(0) = 0, f,(- 1) = 0. 
Es wurden fiir die Polynome f,,(7) die folgenden Ansatze gewahlt : 


Alm) = n( — In) 

San) = 7 |n| 2 —|n)* » 

f(n) = 7°(1 — |n))3 | 

ft) =n. —|n) | 

fon) = n|n\(l —|n)) » | 

Se(n) = 77(1 —|n)) « 
Diese Funktionen sind in Abb. 1 und 2 dargestellt. 


Fiir den ortlichen Wirkungsfaktor erhailt man dann nach Auswertung des Integrales in (3) 


tsa —G) =1— SKF (mo —L)- (6) 


1 Vel. Gl. (6) der in FuBnote 1 von S, 23 zitierten Arbeit, 
2 Ebenda Gl. (9). 
3 Ebenda Gl. (16). 


.B : > Hi i 
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Abb. 2. Die Naherungspolynome f,, f;, f, fiir antisymmetrische Auftriebsverteilungen, Gleichung (5). 


Wir bezeichnen diese Funktionen F,, als EinfluBfunktionen der Grundpolynome f,. Diese 
EinfluBfunktionen sind allgemein giiltig. Zur Ermittlung des induzierten Seitenwindes brauchen 
wir diese Funktionen F,, nur mit den entsprechenden Konstanten K, zu multiplizieren, die fiir 
eine bestimmte antisymmetrische Auftriebsverteilung durch das Naherungspolynom (4) gegeben 
werden. Im allgemeinen geniigt es, nur drei der Grundpolynome fiir die Approximation zu ver- 
wenden. Handelt es sich um antisymmetrische Auftriebsverteilungen mit einem ausgepragten 
Maximum im inneren Teil des Fliigels (Rumpfeinflu8), so kénnen die ersten drei Polynome 
fs fos fz herangezogen werden; liegt das Maximum weiter auferhalb (z. B. Querruderaus- 
schlag), so werden die Polynome f,, f;, f; in den meisten Fallen eine gute Naherung liefern. 
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Abb. 4. Kurven F, = konst. in der (7, ¢)-Ebene. 
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Abb. 6. Kurven F, = konst. in der (ny, )-Ebene. 
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c-¢, 
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Abb. 7. Kurven F, = konst. in der (n, ¢)-Ebene. 


Abb. 8. Kurven F, = konst. in der (n, ¢)-Ebene. 


XXI. Band 1953. Jacobs: Ein vereinfachtes Verfahren zur Berechnung des induzierten Seitenwindes usw. 29 
ag ee a ih 


Die Einflu8funktionen F, sind nun fir die gesamte Querebene —1 < y<-+1 und 
—1<¢€—£, < +1 berechnet worden. Die Ergebnisse sind in Kurvenform in Abb. 3 bis 8 
dargestellt, und zwar wurden die Kurven F,, = konst. (n=1,..., 6) in der [n, (€ —C. )]-Ebene 
gezeichnet. ' ; 


Abb. 9. Zusatzliche Auftriebsverteilung am Schulterdecker beim Schieben. Rumpf: Kreisquerschnitt b/D = 7,5; 


2 F “ D 
Fligel: Rechteck A = 5, Can = 2%; Fliigelhochlage z, an 1; ausgezogene Kurve: Rechnung nach Multhopp; 
gestrichelte Kurve: Naherungspolynom nach’ Gl. (4); SS — 1,88 f, + 2,99 f, — 2,21f,. 
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Abb, 10. Induziertes Seitenwindfeld (Cg 5,,] = konst.) hinter einem Schulterdecker beim Schieben, vgl. Abb, 9. 


Fiir eine vorgegebene antisymmetrische Auftriebsverteilung, fiir welche die Konstanten K,, des 
Naherungspolynoms ermittelt worden sind, lassen sich sodann mit Hilfe der EinfluBfunktionen 
F,, auf einfache Weise das ganze praktisch interessierende Seitenwindfeld und der drtliche 
Wirkungsfaktor des Seitenleitwerkes bestimmen. Da der induzierte Seitenwind sich stark mit 
dem Abstand von der Wirbelschicht dndert, ist bei der Bestimmung des Wirkungsfaktors des 
gesamten Seitenleitwerkes iiber die Héhe des Seitenleitwerkes zu integrieren. Bei Beriick- 
sichtigung der verdnderlichen Tiefe des Seitenleitwerkes gilt fiir den gesamten Wirkungsfaktor 
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Abb. 11. Auftriebsverteilung mit Anstellwinkelsprung in Fligelmitte (Fliigel mit V-Stellung). 
Fligel: Rechteck, A = 5, cZ . = 2m; ausgezogene Kurve: Rechnung nach Multhopp; 
gestrichelte Kurve: Naherungspolynom nach Gl. (4); y/o = 2,00 f, — 2,82f, + 3,12f,. 
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Abb. 12. Induziertes Seitenwindfeld (Kurven 6g/o = konst.) fiir Auftriebsverteilung 
mit Anstellwinkelsprung in Fliigelmitte (Fliigel mit V-Stellung); vgl. Abb. 11. 


4. Beispiele. An einigen Beispielen mége die Anwendung des beschriebenen Verfahrens 
erldutert werden: 


a) Fliigel-Rumpfanordnung. Abb. 9 zeigt die Zirkulationsverteilung, welche an einem 
schiebenden Schulterdecker durch den Rumpf erzeugt wird! *. Fiir die Anniherung der Auf- 
triebsverteilung wurden die drei ersten Polynome f,, fo, fg verwendet, die die tatsichliche Auf- 


1 H, Schlichting u. W. Frenz, Jahrbuch 1941 d. deutschen Luftfahrtforschung, S. I 300. 
2 W. Jacobs, Jahrbuch 1941 d. deutschen Luftfahrtforschung, S. I 165. 
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triebsverteilung sehr gut wiedergeben. Das zu dieser Verteilung gehérende induzierte Seiten- 


windfeld (¢; = konst.) ist in Abb. 10 dargestellt. 


4 b) Stiickweise antisymmetrische Anstellwinkelverteilungen. In Abb. 11 haben 
wir die Auftriebsverteilung dargestellt, die sich ergibt, wenn sich der Anstellwinkel in Fliigelmitte 
(7 = 0) sprunghaft um den Betrag von Ax = 2c andert. Dieser Fall liegt z. B. bei einem 
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Abb. 13. Auftriebsverteilung infolge Querruderausschlag. Fligel: Rechteck,A =5, 4 ., = 27; 


ausgezogene Kurve: Rechnung nach Multhopp; gestrichelte Kurve: Naherungspolynom nach Gl. (5); 
ylo = 0,432 f, — 2,27 f, + 5,75 fy. 
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Abb. 14. Induziertes Seitenwindfeld (Kurven 6¢/o = konst.) infolge Querruderausschlag, vgl. Abb. 13. 


schiebenden Fliigel mit V-Stellung vor. Die genaue Auftriebsverteilung ist mit der Naherungs- 
verteilung (Polynome f,, fs, fg) verglichen. Die Ubereinstimmung ist befriedigend. Die Kurven 
konstanten induzierten Seitenwindes ds/o = konst. fiir die Naherungskurve enthalt Abb. 12. 


In Abb. 13 ist die Auftriebsverteilung dargestellt, die sich bei einem antisymmetrischen 
Anstellwinkelsprung nur in der 4uBeren Fliigelhalfte ergibt. Dieser Fall entspricht etwa einem 
Querruderausschlag. Bei der Naherungsdarstellung wurden wieder nur die drei letzten Polynome 
fa. fs» fe verwendet. Die Ubereinstimmung ist befriedigend. Das induzierte Seitenwindfeld 


(ds/o = konst.) fiir die Naherungskurve zeigt Abb. 14. 
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5. Zusammenfassung. Es wird ein Verfahren angegeben, das es gestattet, das induzierte 
Seitenwindfeld hinter einem Flugzeug fiir antisymmetrische Auftriebsverteilungen auf einfache 
Weise zu berechnen und damit die Wirksamkeit des Seitenleitwerkes zu bestimmen. Dabei wird 
die antisymmetrische Zirkulationsverteilung ersetzt durch ein Naherungspolynom mit freien 
Konstanten, die sich leicht fiir eine vorgegebene Verteilung bestimmen lassen. Die EinfluB- 
funktionen der einzelnen Polynome auf den induzierten Seitenwind sind berechnet und in Kurven- 
form dargestellt worden. Die ganze Aufgabe der Berechnung des Seitenwindes fiir eine vor- 
gelegte antisymmetrische Auftriebsverteilung beschrankt sich jetzt im wesentlichen auf die 
Bestimmung der Konstanten des Naherungspolynoms. 


(Eingegangen am 28. Mai 1952.) 


Anschrift des Verfassers: Dr. habil. W. Jacobs, Stockholm-Bromma (Schweden), Stora Mossens Backe 6, 


XXI. Band 1953. Schmutz: Koeffizientenbedingungen zur Kontrolle des Dampfungs grades usw. 33 


Koeffizientenbedingungen zur Kontrolle des Dampfungsgrades bei 
Ausgleichvorgingen (verallgemeinerte Hurwitzbedingungen). 


Von O. Schmutz, 


Lf Ubersicht. Fiir Systeme linearer Differentialgleichungen mit konstanten Koeffizienten 
geben die komplexen Wurzeln p,—b, +iw,, bzw. die reellen Wurzeln Py =65, der ,,Stamm- 


gleichung’* den Verlauf der Einzelschwingungen C, ert) byw, C, e”” * eines Ausgleichs- 
_vorganges an, Hurwitz stellte ohne Ermittlung der Wurzeln die Koeffizientenbedingungen dafiir 
auf, daB8 nur abklingende Schwingungen mit 6, bzw. b,< 0 auftreten (Ziff. 3). Die Hurwitz- 


_bedingungen lassen sich bekanntlich in einfacher Weise auf den Fall erweitern, da® ein Maximal- 


x= CeP' die Stammgleichung (charakteristische Gleichung) 


dann klingen alle angestoBenen Eigenschwingungen ab. Beim Aufhéren 


wert der Abklingzeit vorgeschrieben ist (Ziff. 4). Dariiber hinaus werden nachstehend (Ziff. 6) 
mit Hilfe bekannter Kigenschaften der Polynome mit komplexen Koeffizienten Koeffizienten- 
bedingungen angegeben, die sicherstellen, dafB Dampfungsgrad bzw. logarithmisches Dekrement 
samtlicher Einzelschwingungen eine vorgeschriebene Grenze nicht unterschreiten (Ziff. 8). Als 
Grenzfall ergeben sich die Koeffizienzenbedingungen dafir, da® alle Wurzeln reell und negativ 
sind, das also der Ausgleichsvorgang rein aperiodisch erfolgt (Ziff. 9 und 12). SchlieBlich wird 
noch der Grenzfall behandelt, daB alle Wurzeln gleichen Betrag besitzen (Ziff. 13). 


2. Einleitung. Andert sich in einem linearen System eine StérgréBe um die Abweichung 
s= s(t) von einem Gleichgewichtszustand, dann andert sich eine von dieser Stérung abhangige 
Gri Be — in der Regelungstechnik die RegelgréBe — um die Abweichung x. Fiir die Abweichung x 
gilt dann folgende Differentialgleichung : 

dx d? x ? dns, ds, 


ea lg foes 1 ee ot seas 
Ui Ost ee SUR ar re "ae S0 OL a ae Sn 


(1) 
mit konstanten Beiwerten 
tne tr aon.) und, Ie Ser «ar 


Um sie zu lésen, sucht man zunachst n linear voneinander unabhangige Lésungen der verkiirzten 


Differentialgleichung 


dx | ax “A d" x ia 


Sots] mG, ' “2 de Fe On oe ) 


und erganzt ihre sich aus diesen linear aufbauende allgemeine Lisung, die n willkiirliche Kon- 
stanten besitzt, durch eine partikulare Lésung der vollstindigen Gleichung (1). Die n Konstanten 
erméglichen die Anpassung an die Anfangsbedingungen. Fiir (2) liefert zunachst der Ansatz 


dy+ 4, p+ a, p> +°*:+a, p"=0. (3) 


Diese Gleichung besitzt n Wurzeln fiir p von der Form p=6, + iw, 
bzw. p,=b,. In unserem Falle der reellen Koeffizienten a), a, usw. 
sind die komplexen Wurzeln paarweise konjugiert komplex. Auf den 
Sonderfall mehrfacher Wurzeln (mehrerer gleicher Wurzeln) soll zu- 
nichst nicht naher eingegangen werden. 

Sind alle reellen Anteile b, der Wurzeln negativ, liegen sie also in 
Abb. 1 in dem linken nicht schraffierten Teil der Gaufschen Zahlenebene, 


7 F 5 4 ‘ ; ] he Abb. 1. Stabilitiit, wenn jede 
der Stérung, bzw. bei einer von einem bestimmten Zeitpunkt ab kon- Warzel p, imnichtschraffierten 


stanten Stérgriéke s strebt die GréBe x einem konstanten Endwert zu. Gebiet der komplexen Zablen- 
* . ebene liegt. 
Das System ist stabil. 
3. Stabilitatsbedingung von Hurwitz. Ohne Aufsuchung der Wurzeln der Gleichung (3), was 
bei Stammgleichungen hdherer Ordnung nicht mehr direkt méglich ist, konnte Hurwitz die 


3 
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Koeffizientenbedingung?! fiir diese Stabilitat angeben: 


literate ae 


Gy gl, Gage, sere 


Ree Oath CPO OP Ost O) 080 CO" Gurr oe (Od 


Sind alle Beiwerte a, a), 4, °** @, positiv und die ,,Abschnittsdeterminanten“ der eingerahmten 
quadratischen Schemata gréfer als Null, dann ist Stabilitat vorhanden. 

Es treten n—1 Determinanten auf mit der Hauptdeterminante vom Grade n—l. Als 
Beispiel fordert bei n= 3 die Hauptdeterminante a, a,—a,a,>0. 


1 
Man kann in die verkiirzte Differentialgleichung auch x= C e!'1 einfiihren, also p = re setzen. 
Es wird dann die Stammgleichung erhalten 
Qn + On1 9+ 4n2 P+ °° +4, gq" 14+ a) q"= 0. (5) 


Besitzt eine Wurzel p, negativen Realteil, dann hat auch q, einen negativen Realteil. Der 
Determinantenbedingung (4) ist also die folgende gleichwertig?: 


| } 
I |an—| Gn—3| ns Ons eee 
An Gas Cae] ag=6 ee F, 
0 Urast cUprms|| | Upese | omob t 
0 Gp Coon eae 
alr iN Bie nets tees MSM chine Ree cans Re j 


4, Forderung einer maximalen Abklingzeit. In der Regelungstechnik kann man sich nicht 
damit begniigen, daB die angestoBenen Eigenschwingungen irgendwann abklingen, sondern sie 
miissen mit einer gewissen Sicherheit und geniigend rasch abklingen. Fir die Wurzel 
P,=), +1, ist — 1/b,= T, die Ab- 
klingzeitkonstante der zugehérigen 
Teilschwingung. Schreibt man einen 
Maximalwert Ty) = —1/b) vor, so 
bedeutet dies: alle |b,| miissen gréBer 
sein als |b |, d.h. alle p,-Werte 
miissenin der komplexen Ebene links 
des Punktes 0’ liegen, der durch den 
Zeiger by festgelegt ist (Abb. 2). 

Der vom Punkt 0’ nach dem Punkt 
/ P, gezogene Zeiger q,= p,—b,) muB 
eh one APES crore einen negativen Realteil besitzen. 


derTeilschwingungen unterschreiten Abb. 3. Teilschwingungen a, b, c Es ist also die bekannte Koordinaten- 
einen Grenzwert, wenn die Wurzeln 7 


‘ i : gleicher Abklingzeitkonstante T, aber 7 3 = ;, 
Py bzw. q, im nichtschraffierten verschiedenem Dampfungsgrad verschiebung A Pea =i by AE LA 
Gebiet liegen. bzw. verschiedenem w. nehmen, wobei man fiir die Stamm- 
gleichung erhalt 
Pee awe | 
@y+ 4\(q+ bo) + (q+ bo)? + + Gn(q-+ bo)” = 0. (7) 


Das Ausmultiplizieren kann bei gegebenen Zahlenwerten a, rasch nach dem Horner-Schema 
vorgenommen werden, und man erhilt 


Ag aA dang pe, aed BO (8) 


Fiir diese Koeffizienten A), A,...A, muB nun wieder die Hurwitzbedingung erfiillt sein. 


* W. Schmeidler, Vortrage iiber Determinanten und Matrizen mit Anwendungen in Physik und Tech- 
nik, S. 68ff. Berlin 1949. 


® A. Leonhard, Die selbsttatige Regelung S. 186, Berlin 1949, 
3 W. Schmeidler a. a. O., S. 71 


} 
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Fir den Grenzfall T,—=—1/b)= T, miissen in Abb. 3 die Teilschwingungen bei entsprechen- 
dem Ma8stab entweder mit der Exponentialkurve a zusammenfallen oder als abklingende 
Schwingungen zwischen der Kurve a und ihrer an der t-Achse gespiegelten Kuive (gestrichelt 
gezeichnet) liegen wie die Kurven b und c. Die Kurve c entspricht hohen w-Werten, wobei das 
Verhiltnis zweier aufeinander folgender Ausschlige nur wenig kleiner als 1 ist. Es + hier also 
nur eine geringe Sicherheit gegen dauernde Pendelungen vorhanden. 


5. Dimpfungsgrade ! der Einzelschwingungen und Mindest-Dimpfungsgrad fiir den Gesamt- 
vorgang. Fiir die Stammgleichung zweiten Grades einer Kinzelschwingung 


a + a p +a, p?= 0 (9) 
lauten ihre zwei Lisungen 
a ati Wie j 5 ot 
Pi2= sl Sn = | | a — 1) =b+iw=\|p,je*"*. (10) 
Hierbei sind die folgenden Abkiirzungen cingefithrt: 
= pal — ape 
eee Wiss EO pao l (11) 
eee) ses Dampfungsgrad . (12) 
2 Vayay 
ee yas 
tg Oe et b aS | iF Ie. (13) 


Die Richtungswinkel a,,, der Zeiger p, und p, sind nur von dem Dampfungsgrad D abhingig. 
Das Amplitudenverhaltnis fiir eine 


Halbschwingung betragt him 
ba Ef X mu 
a : — —— lee, 
pe 2 w : V1/D2—1 os a, tg B Purafftufe oo 
== == = <i} J 
, (14) Ce Wo 
mit 4 ke 
(Si ec, 0 oe aad 
teB=|/ m1 und cosf=D. (14a) je 
: > re ¢ Z 
Das logarithmische Dekrement L’ fiir die 2e7//e"® 
Halbperiode bzw. L fiir die ganze Periode 
ergeben sich damit zu Abb. 4. Die Winkel zur Abb. 5. Die Daimpfungsgrade der 
1 ai Kennzeichnung der Lage der Teilschwingungen iiberschreiten 
Lge = ss = konjugiert komplexen Wurzeln einen Grenzwert D = cos f, wenn 
(15) jug P B 
sy erie : p, und p,. alle Wurzeln innerhalb des nicht- 
|/ ss schraffierten Winkelraumes liegen. 
D2 


Der Dampfungsgrad Null entspricht der ungedimpften Dauerschwingung, der Dampfungs- 
grad 1 dem exponentiellen Abklingen und der Dimpfungsgrad 0,5 einer gut gedampften Schwin- 
gung mit 

Pp ae Vahey adie = AIBN. 1418'S /607,, a = 0103'S: a3 


Die folgende Amplitude ist also in diesem Fall nur rund 1/6 jeder vorhergehenden. 

In Abb. 4 ist die Lage der Wurzeln p, und p, dargestellt, wobei der Winkel fp =a —a, er- 
sichtlich ist. Dieser Winkel ist positiv bei den praktisch interessierenden Fallen eines negativen 
Realteiles. 

Mit dem Dampfungsgrad D ist nach (15) das logarithmische Dekrement, nach (14) das 
Amplitudenverhaltnis von zwei aufeinanderfolgenden Halbschwingungen und nach (13) der 
Richtungswinkel « der Wurzeln bzw. der Winkel f bestimmt. 

Fiir eine gute Regelung miissen nun alle Ausgleichsschwingungen Dampfungsgrade besitzen, 
die grifer als ein vorgegebener Mindestdimpfungsgrad sind. Alle Wurzeln miissen also innerhalb 
des in Abb. 5 nicht schraffierten Winkelraumes liegen, der durch einen der Winkel 


IU = 5 1G 
a SOS, ge eo Vas 


1 Regelungstechnik-VDI-Fachausschuf, 5. 48, Berlin 1944. 
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definiert werden kann}. Ob diese Bedingung erfiillt ist, 14Bt sich nun ohne Bestimmung der 
Wurzeln der Stammgleichung direkt durch Determinantenbedingungen fiir ihre Koeffizienten 


angeben. Die Hurwitzbedingungen bilden nur den Sonderfall mit « = > bzw. B = + bzw. y =0. 


6. Eigenschaften der Polynome mit komplexen Koeffizienten ?. Es werden fiir die Aufstellung 
der gewiinschten Koeffizientenbedingungen gewisse Eigenschaften der Polynome n-ten Grades 
mit konstanten aber beliebigen komplexen Koeffizienten Ay, A,,...A, benutzt. Die Wurzeln 
besitzen simtliche dann und nur dann positive Imaginarteile, liegen also in der oberen Halb- 
ebene, bzw. ,,links der Richtung der positiven reellen Achse“‘, wenn die folgenden Abschnitt- 
determinanten positiv sind. Werden hierbei die Abkiirzungen verwendet 


(hl) = 5 (Ay A) — Ay Ae) = | Aa | Ail sim (er — &) (16) 
A= Ape PE AS Ae (17) 
so gilt fiir n=2: 
Lo) | (02)) 9, (18) 
(02) a12)| 
fir n=3 
KG) (02) (03) 
ee (03)-+ (12)) (13)|>9- (19) 
(03) (13) (23) 
fiir n= 4 
op] 2), | 08) 
(02) (03) + (12) (04) + (13)) (14)|S 0, (20) 
(03)  (04)+ (13) (14) 4 (23)| (24) 
(04) (14) (24) (34) 
Fires 
| (on) (02) (03) (04) (05) 
| (02) (03) + (12) (04) + (13) (05)+ (14); (15) 
(03) (04)-+ (13) (05)+ (14)+ (23)] | (15)-+ (24)} (25)|> °° Ce 
(04) (05) + (14) (15) + (24) (25) + (34); (35) 
(05) (15) (25) (35) (45) 


Das Bildungsgesetz ist nun leicht zu iitbersehen, wenn beachtet wird, daB bei den in Summenform 
auftretenden Gliedern die Quersumme jedes Klammerausdruckes gleich gro8 ist. 


7. Ubergang zu den Bedingungen fiir vorgegebene Mindestdimpfung D bzw. fiir den ent- 
sprechenden Winkelgrenzwert 3. Sind die in Ziff. 6 aufgestellten n Koeffizientenbedingungen 
fiir das Polynom n-ter Ordnung erfiillt, liegen also simtliche Wurzeln links der reellen Achse, 
dann mu fiir alle Wurzeln 0<a<a@ sein. Zum Zwecke der angestrebten Verallgemeinerung 
wird eine Koordinatendrehung p= qe'? vorgenommen (0 <fB <a; B= konstant), also die 


1 A, Leonhard a. a. O., S. 198 ff. 
2 W. Schmeidler, a.a.O., S. 62. 
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Stammgleichung statt in Koordinaten der p-Ebene jetzt in Koordinaten der q-Ebene ausgedriickt : 
. 5) . : q 

Qy+ a, eF gta, eh P41... 41a, eat ghee 7 (22) 

AvP A, g+ 4,9 >: +A, = 0. (23) 


Der Punkt p=|ple® der p-Ebene geht iiber in den Punkt C= | p | ei( —?) der q-Ebene. Werden 
| die obigen Ungleichungen fiir die in q ausgedriickte Stammgleichung mit den kom lower Koeffi- 
zienten Ap, A,, A, erfiillt, dann liegen die Wurzeln in der Halbebene links der Nene der 
reellen Achse der q-Ebene. Dieser Richtung entspricht in der p-Ebene die Richtung aaa In 
der a interessierenden p-Ebene liegen die Wurzeln dann nur in der Halbebene links der Ricbeane 
=p . 

Beriicksichtigt man ferner, daB die Wurzeln in der p-Ebene nur konjugiert komplex auftreten 
-kénnen, so treten sie dann hier weder in der Halbebene rechts der Richtung «=f noch in der 
Halbebene links der konjugierten Richtung «= —f auf, so daB die Wurzeln nur in dem in Abb. 5 
nicht schraffierten Raum liegen kénnen. Sie kinnen also nur Richtungswinkel miiecken 
a=a—B und a=2-+f besitzen. 

In dem bisher zuganglichen Schrifttum ist nur der Sonderfall f =90° behandelt!. Hierbei 
ergeben sich aus den Ungleichungen (18) bis (21) Determinantenbedingungen, die zwar weniger 
itbersichtlich wie bei Hurwitz sind, aber auch weniger Auswertarbeit erfordern. Doch gibt es 
auch fiir die schnelle Auswertung der Hurwitzdeterminanten gecignete Methoden?. 


8. Koeffizientenbedingungen fiir eine Mindestdampfung (. Die Bezichungen (16) bis (21) sind 
jetzt sinngem48 auf die komplexen Koeffizienten A,—|A,|e'”® der Gleichung (23) anzuwenden. 
Gleichung (16) ergibt 


(kl) = |A,| |A)| sin(l —k) B, (24) 
fur — 2 
| aya, sin B. dy a, sin 2B 
este ees >0, (25) 
a,a,sin2B a, a, sin B 
fir n=3 
I 4 ay sin B dy a, sin 2B dy d3 sin 3B | 
a eee | | 
I dya,sin2B aya; sin 3B + a, a, sin B| a,a,sin 2B |> 9° (26) 
| dy a; sin 3B a, a, sin 2B ay ay sin B | 


Damit ist das Bildungsgesetz dieser Determinanten aus ihrer allgemeinen Form ersichtlich. Es 
ergeben sich hier n Determinanten mit einer Hauptdeterminante n-ten Grades. 

Liegen r Wurzelpaare an der Stabilitatsgrenze, dann sind die 
2r héchsten Determinanten der Stabilitatsbedingung (5) gleich 
Null’. Analog wird in Erweiterung dieses Satzes gefolgert: Liegen r 
Wurzelpaare an der Mindestdimpfungsgrenze, dann sind die 2r 
héchsten Determinanten der Mindestdampfungsbedingungen gleich 
Null. Umgekehrt kénnen unter Nullsetzung dieser 2r héchsten 
Determinanten die Koeffizienten so bestimmt werden, daf§ 2r 
Wurzelpaare gerade die geforderte Mindestdampfung besitzen. 


9. Koeffizientenbedingung fiir eine aperiodische Dampfung. Fir 
den Grenzfall 6 > 0 liegen, wie Abb. 6 andeutet, alle Wurzeln p 
auf der negativen reellen Achse, wobei im aperiodischen Grenz- Abb. 6. Die Dampfungsgrade der 


“ Teilschwingungen werden gleich oder 
fall auch Doppelwurzeln bzw. Mehrfachwurzeln auftreten kénnen. snc e gis 1 bzw. alle Wurzeln werden 


Es wird hierbei sinvB=vf, wobei mit den gemeinsamen Fak- negativ reell, wenn diese in dem 
eae nichtschraffierten Winkelraum mit 
toren B dividiert werden kann. Man hat also dura Renaiibescanigiocs0 lurvoas 


1 W. Schmeidler, a. a. O., 5. 66. 
2 W. Schmeidler, a. a. O., S. 68. 
3 H. Cremer, Z. angew. Math. Mech, 25/27 (1947) S. 160, 
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fiir t= 2 : 
| a, a, | 2 ty | 
>0, (27) 
2M) M% 4 Ay | 
tur == 5s 2 
4% a, | 2 My Ay | 3 dy As 
2)4, 34,43, a, ay 2, Gs > (25) 
3 dy a, 2 A, Gs My As 
fir n=4 ; 
| May . 2 Ay Ay | 3 dy Gs | 4 dy My 
Dee 3 Gy A, ay a, 4) 4,+ 20,43) 3a,a, Se (29) 
34,4, 4a,a,+2aa, 34, a,+ a, 43 | 2 Gy Ay 
| 4 ay cal aa ay 2 Ay My Bs Uy 


10. Eigenschaften eines Polynom mit komplexen Koeffizienten nach Bilharz!. Diese andere 
Darstellung soll angefiihrt werden, da sie ein ithersichtliches Schema ergibt. Fiir die rechnerische 
Auswertung erscheint es aber weniger geeignet. 

Die Wurzeln von (23) mit beliebigen komplexen Koeffizienten A,—= A, +iA) haben dann 
und nur dann alle einen negativen Realteil, wenn A, ~0 und die nachfolgenden Abschnitts- 
determinanten gerader Ordnungszahl positiv sind: 


=>) eee | A 


| | Ay Ay 2 3 4 5 6 —A, 
| 40 Ay fay AN Me Apes 
| Ay A Nee A i ae A, Aes Cah ee Se 
| ee AD iis Akt Wed A Ale 
0 0 ve Ae Ear Ad AG At 
0 0 A A A ed eed A 


11. Ubergang zur Dampfungsbedingung fiir einen vorgegebenen Grenzwert D der Dampfungs- 
grade unter Verwendung des Bilharzschen Satzes, Es wird wieder eine Koordinatendrehung 


P= |p |e = qei” vorgenommen (<> <= » vgl. Abb. 4) : 


Gy a, ebY g- a, ec 2? g?-+: a, ef"? gh = 0. (31) 
Dem Punkt p der p-Ebene entspricht ein Punkt g= |p| e(*—) der neuen q-Ebene. Wir wenden 
nun den Billharzschen Satz (30) auf das q-Polynom mit A,—<a, e’”” an und erhalten dann die 
Bedingungen dafitr, da alle q-Wurzeln nur in der Halbebene links der positiven imaginaren 


Richtung mit dem Richtungswinkel > =o«-—y in der q-Ebene liegen. Dies bedeutet, daB® die 
p-Wurzeln dann nicht in der in Abb. 5 schraffierten Halbebene rechts des Richtungswinkels 
a= +y liegen kénnen. Da nun wieder bei reellen Koeffizienten a, die p-Wurzeln nur reell 
oder konjugiert komplex auftreten kénnen, miissen alle diese Wurzeln in dem nicht schraffierten 


Winkelbereich 26 liegen mit 6 = ~ —y. Fir die Koeffizienten ist nun 


1 W. Schmeidler, a. a. O,, S. 70, 
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A, =a, "7 = A,+ iA; =a,cosy y +ia,siny y, (32) 
Vee. a 7 It 
A, = a, cos» y = a, cos y le — p) =a, (:0s » = cos y B+ sin y + sin »®) (33) 
ee . 3 . 4 Bs fy 5 
| Ar= ssiny y= a, siny(F—B) = a,(siny F cos f — cos’ sin vf), (34) 
Wir wenden nun die Billharzsche Bedingung (30) auf diese Werte an: 
| ay a, cos B | a,cos2B ascos3B| a,cos4f ----- 
0 a, sin B a, sin 2 B az;sin3B | asin4B ----- 
Orata. a, cos B a,cos2B | a;cos3B ----- 
0 0 a, sin B a,sin2B| a,sin3B ----- gi oe) 
0 A a, cos B a,cos2Bh -::-- 
0 0 a, sin B a,sin2B --+:: 


Die n Abschnittsdeterminanten gerader Ordnungszahl miissen alle gréfer als Null und auBerdem 
a S0 sein, dann liegen alle Wurzeln im Winkelbereich 2 £. 

Bedingung 35) besteht aus n Determinanten mit der Hauptdeterminanten der Ordnungszahl 
2 n. Das tibersichtliche Schema ist also durch eine doppelt so hohe Ordnungszahl gegeniiber den 
Bedingungen (25) bzw. (26) erkauft, erscheint also weniger giinstig fiir eine praktische Aus- 
wertung. 


12. Koeffizientenbedingungen fiir die aperiodische Dampfung aus (35). Die Bedingung, daB 
die Wurzeln der Stammgleichung alle reell und negativ sind erhalt man durch den Grenziibergang 
B—>0 aus der Ungleichung (35) 


= ) 
| a ay | ay (sy ay as Bote waht 
0 a, | 2a, 3a5, 4a, 5a,| ----- 
0 a% ay Gy | as Will FO SOO 
i < 
0 9 4 2a,) 3a, 4a) +> >0 mit ays 0. (36) 
0 0 sr) ay a, ae sS 
0 0 0 ay 2D As 33 As euetersts 
BR ones as Canes he ahs, nee ote a tae Celene que os j 


Zur Kontrolle wurden fiir n=3 aus beiden Ungleichungen (27) und (36) die drei im Ergebnis 
iibereinstimmenden Bedingungen fiir den aperiodischen Fall ausgerechnet : 


De=a,a,>0, ' 
D, = ay (3 dy a, 4g — 4a) a2 + a?a,) >0, i (37) 
Dz = dy 43(18 dy a, a, a, + a? a2 — 4a3a, — 4a, a} — 27 a§a3) >0. J 


Da alle Koeffizienten a, bereits nach Hurwitz fiir 0 < Bp <> positjv sein miissen, geniigt es, 


da®B auBerdem nur noch die Klammerausdriicke der Determnianten D, und D;>0 sind. Die 
Formel kann am Grenzfall der zusammenfallenden Wurzeln leicht nachkontrolliert werden, 
wobei alle Determinanten (36) von der Ordnungszahl gréBer als 1 zu Null werden. Das Beispiel 
mit drei Wurzeln p; = p. = p3 = — 1, namlich 


a) +a, p+ a, p?+a,p?=1+3p+3p? + p= (1-+p)r=—0 (38) 


zeigt, daB die Koeffizienten dieser Gleichung die Determinanten D, und D, von Bedingung (37) 
zu Null werden. 
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13. Sonderfall, da8& alle Wurzeln gleichen Betrag und gleiche Dimpfung besitzen. Wir bringen | 


die Stammgleichung (3) auf die Normalform, indem wir sie zunachst mit a) dividieren. 


17 OPO = 
iaritarar ss iat ios. (39) 
Mittels der Substitution 
On vo n pa f dy 2 
peal g's NL a (40) | 
wird schlieBlich erhalten 
Me Res ayia | 
ead! rece (ya Pr Sige ye) +g =o | (41 | 


It gt gt tog ge=O. 


Der erste und letzte Koeffizient ist hier gleich eins geworden. Fir die neuen Koeffizienten ergibt — 


sich also 


eee ae mit Kei 2, 30. ee, n—1. (42)8 


a an 


Die neue Verdnderliche q ist nur ein positives reelles Vielfaches von p. Damit sind auch die | 


q- Wurzeln ein reelles Vielfaches der p-Wurzeln. 
Die Normierung des Polynoms zweiten Grades (6) ergibt 


1+2Dq+@=1+2qceosh+ qg= 0 (43) | 


mit den Wurzeln vom Betrage 1: 


=u] rlieeen t(a+ 
Pea p(i iV pe ae ee (44) 


Die allgemeinen Gleichungen fiir Wurzeln gleichen Betrages und gleicher Dampfung haben dann 
die Form 


(1+2Dq+ q)r? = 0 fiir gerades n , (45) 


n— 1 


2 


(1+ 2D q+ q’) 


Bei ungeradem n muB8 eine reelle q-Wurzel vorhanden sein, und zwar ergibt der Wurzelwert 
q = —1 auch fiir die Polynome (46) symmetrische Koeffizienten 


(1+ q)=0 fiir ungerades n . (46) 


C=C, “nit v= 1, 273... 9¢n— Diz. 
Es ergeben sich damit folgende Beispiele von normierten Stammgleichungen: 
IM=s= (5 Ges 2. 1 ae 0; 


} 
n= 3, Iga gee 10, | 
Mien Sead) oe Gees t (47) 
n=5, I1+q+2¢@+2¢+q+q=0, | 
n=6, I1+3¢+3¢4 ¢=0; | 
D=0,5, n=2, 1+q+q=0, a) 
n= 3, 1+2q+2¢+q= 0, | 
n= 4, 14+29q+3¢+2q+ q=0, ( (48) 
n=0, 1+39--59¢ + 59+ 3g Gg = 0; | 


DS by ane ye gi 0), } 
n= 2, 1+2q+q@=0, | 

n=3,1+3q+3¢@+q@=0, (49) 
J 


n= 4, 1+ 4q+ 6¢+ 4q°+ gt = 0. 


Sind bei einer Regelaufgabe n—1 Konstanten C},¢,... wahlbar, dann wird man die Wahl 
nach Festlegung des erforderlichen Mindestdimpfungsgrades D z.B. 0,5 so treffen, da® die 
normierte Stammgleichung mit der entsprechenden Gleichung (48) identisch ist. 
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14. Koeffizientenbedingung, da8 siimtliche Wurzeln auferhalb eines Kreises liegen’. Auf 
iese bekannte Méglichkeit soll nur hingewiesen werden. 


15. Beurteilung der Ergebnisse. Die Nachpriifbarkeit des Dampfungsgrades an Hand der 
oeffizientenbedingungen stellt eine wiinschenswerte theoretische Erkenntnis dar. Praktisch 
ewahbren sie wohl zu wenig Einblick iiber den Verlauf und die BeeinfluBbarkeit des Ausgleichs- 
organges. Die langwierige zahlenmafige Auswertung fiir gegebene Koeffizienten gibt zunachst 
ur Auskunft, ob ein gewiinschter Dampfungsgrad eingehalten wird. Ist ein Teil der Koeffi- 
ienten noch wahlbar, so kénnen sie so bestimmt werden, da8 ein Teil der komplexen Wurzeln 
erade die Mindestdampfungsgrenze erreichen. Der tatsachliche Verlauf des Ausgleichsvorganges 
angt nicht etwa nur von den Wurzeln der Stammgleichung ab, sondern es wird die GréBe der 
eilschwingungen auch von den Koeffizienten der rechten Seite von (1) mitbestimmt. Diese 
ndern sich bei Variation der Koeffizienten der Stammgleichung meist ebenfalls. 


(Eingegangen am 16. Juni 1952.) 


Anschrift des Verfassers: Dipl.-Ing. Oskar Schmutz, Karlsruhe, Technische Hochschule, 
Elektrotechnisches Institut. 


1 W. Schmeidler, a. a.O., S. 71; ferner ein Vortragsmanuskript vom 2. 4. 1948 (in unseren Handen) 
‘on H. Tischner (Tiibingen) ,,Einfiihrende Betrachtungen iiber Stabilitat von Regelkreisen. 
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Analogie zwischen Wasserschwall und Verdichtungssto8. 
Von E, Becker. 


1, Einleitung. Schon seit langerer Zeit ist bekannt, daB die Gleichungen fiir die Str6mung; 
flachen Wassers iiber einer horizontalen, ebenen Sohle in erster Naherung den Gleichungen fiir die: 
zweidimensionale Bewegung eines kompressiblen Gases gleichen. Die theoretische Fundierung; 
dieser sogenannten hydraulischen Analogie verdankt man vor allem Preiswerk*. Es stellte sich. 
allerdings bald heraus, daB die Wasserstrémung eine ganze Reihe von Abweichungen gegen-- 
ber der Gasstrémung zeigt, die einen quantitativen experimentellen Vergleich sehr erschweren.. 
Ziel dieser Arbeit ist eine Diskussion der Abweichungen, die sich besonders beim Wassersprung,| 
dem Analogon des VerdichtungsstoBes, bemerkbar machen. Unter Wassersprung ist hierbei all-» 
gemein eine unstetige Erhéhung des Wasserspiegels zu verstehen. Wandert diese Unstetigkeit: 
relativ zur ruhenden Bodenflache in Wassergebiete hinein, die sich in hezug auf die Bodenflache) 
in Ruhe befinden, so soll insbesondere von einem Wasserschwall die Rede sein. Unter den Wasser-- 
schwallen werden solche mit stetig bleibender Wasseroberflache und solche mit einer Wasserwalze 
(=turbulente Durchmischungszone) unterschieden. 

Ihrer Natur nach lassen sich die Abweichungen zwischen VerdichtungsstoB und Wassersprung | 
in zwei Gruppen einteilen: Die erste Gruppe, auf die schon Preiswerk aufmerksam macht, wird | 
durch das verschiedene thermodynamische Verhalten von Verdichtungssto8 und Wassersprung ; 
verursacht. Die zweite Gruppe ist darauf zuriickzufithren, da®B die Wasserstrémung dreidimen- 
sional ist gegeniitber der zweidimensionalen Gasstrémung. Das heifbt: es treten in ihr Vertikalll 
geschwindigkeiten auf, die kein Analogon in der Gasstrémung haben. AuBerdem hat das Wasser 
zwei Grenzflachen: Bodenflache und Oberflache. An jener wirkt die Bodenreibung ein, an dieser 
die Oberflachenspannung. Als Ergebnis der vorliegenden Untersuchung zeigt sich, daB das Zu- 
sammenwirken von Vertikalgeschwindigkeit und Oberflachenspannung die Ursache fiir zwei bei 
den Experimenten an Schwallen (in obigem Sinne) manchmal beobachteten Erscheinungen ist : 

1. Wenn ein in ruhendes Wasser hineinlaufender Schwall nicht zu stark ist, bildet die Wasser- 
oberfliiche u. U. keine scharfe Unstetigkeit, sondern sie bleibt am Ort des Schwalls stetig, wenn 
auch sehr steil. In diesem Fall ist dem Schwall ein gut ausgebildetes System von Kapillarwellen 
vorgelagert, deren Phasengeschwindigkeit mit der Schwallgeschwindigkeit iibereinstimmt. 

2. Unter der Voraussetzung, daB der Schwall in einer solchen stetigen Erhéhung des Wasser- 
spiegels besteht, beobachtet man bei Erfiilltsein gewisser zusatzlicher Bedingungen (die in dieser 
Arbeit naher formuliert werden) schon vor dem Ort des eigentlichen Schwalls eine allmahliche. 
stetige Erhéhung des Wasserspiegels. Die Kapillarwellen itiberlagern sich im allgemeinen dieser 
stetigen Erhéhung, ohne sie wesentlich zu beeinflussen. Bei Wassertiefen, die kleiner als 0,4.cm sind, | 
beobachtet man in keinem Fall einen stetigen Anstieg der eben beschriebenen Art, sondern nur 
Kapillarwellen. Wird der Schwall zu stark, so verschwinden stetige Erhéhung und Kapillarwellen | 
und es bleibt am Ort des Schwalls eine heftige turbulente Durchmischungszone des Wassers. 


, 


2. Verlustfreie Strémung. Alle Vorgange werden durch Wahl eines geeigneten Koordinaten- | 
systems stationar gemacht. Zunachst miissen drei grundlegende Vernachlassigungen bei den Be- | 
wegungsgleichungen des Wassers formuliert werden, die notwendig sind, um diese auf eine den 
Gleichungen des Gases analoge Form zu bringen: 

I: Die Vertikalgeschwindigkeit und deren Ableitungen werden gleich Null gesetzt. 
II: Innere Reibung und Bodenreibung bleiben unberiicksichtigt. 

III: Die Oberflachenspannung wird vernachlassigt. 

AuSerdem mu8 die Adiabatenkonstante (Verhaltnis der spezifischen Warmen K = p[Cy) des 
Gases, dessen Strémung mit der des flachen Wassers verglichen wird, gleich 2 gesetzt werden. 

Es entsprechen sich dann folgende GréBen (Abb. 1): 


1. Enthalpie t | 
Temperatur O ¢ & Wasserspiegelhohe € ittber dem Boden, 
Gasdichte 0 


_| E, Preiswerk, Anwendung gasdynamischer Methoden auf Wasserstromungen mit freier Oberflache 
Diss, Ziirich 1938, 
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. Gasdruck P a£e (g = Schwerebeschleunigung) , 


. Schallgeschwindigkeit }/x Plo A Grundwellengeschwindigkeit ¢ = Vgc. 


Die wichtigsten sich entsprechenden Gleichungen der verlustfreien Strémung sind folgende 
links stehen jeweils die Gleichungen des Gases, rechts die der Flachwasserstrémung) : 


d(£ @ 
d 
WdW +-~ = 05 WdW + (F ual (d.h. WdW + ¢ de =0), 
W® = 2 g (iy — i) ; W? — 2g(f)—6), 
BO) a 80) ay. acu) | aco) 
Ce a eee er Te er Toca 


W bedeutet hierbei den Betrag der Geschwindigkeit mit den Komponenten u und v (parallel zur 
3odenflache beim Wasser). Der Index 0 bezieht sich auf den Ruhezustand (z. B. py = Kesseldruck, 
9 = Stautiefe). 


Abb. 1. Koordinatenbezeichnung. Abb. 2. Wassersprung (Verdichtungsstof) 
an einer einspringenden Ecke. 


3. Str6mung mit Verlusten. Beim sogenannten schragen Wassersprung bzw. VerdichtungsstoB, 
wie er etwa durch die Strémung in Umgebung einer konkaven Ecke oder durch Anstrémung einer 
ingestellten Platte realisiert ist, gelten folgende Beziehungen (Abb. 2, der Index | kennzeichnet 
>r6Ben vor dem Wassersprung, der Index 2 solche hinter dem Wassersprung) : 


Kontinuitatsgleichung : 


0, Wy = o Wy; 6, Wy = bo Wy; (1) 
Impulssatz parallel zur StoBfront: 
Oz Wy ty a 0, We W,, 5 ey We Wa == ce Wr. W,, 3 (2) 
Impulssatz senkrecht zur StoBfront: 
Wi +pi=0"i,trs GW t$G=4m,+2£9. (3) 


Bis hierher hat man bei obigen Vernachlassigungen vollkommene Ubereinstimmung zwischen 
len gasdynamischen und den Flachwassergleichungen. Daher entsprechen sich auch alle Be- 
iehungen, die allein unter Benutzung von (1) bis (3) abgeleitet werden. Die wichtigsten davon 
verden nun tabellarisch zusammengestellt 1. Mit den Abkiirzungen 


Ao = 02 — 01> AV Vg ae Boal, aes Ge, 
f §& g rs 8 Re 
Ap = P,— Py; Ls le 5 ope OT 
rgibt sich fiir die Geschwindigkeitskomponenten Wt parallel und Wn senkrecht zum Wasser- 
prung 


W,,= W,, » W,,= W,,, 
LOM, 

yw — 9 4P By otoulens & 

ee ae, Se oey 


1 Zur Herleitung der gasdynamischen Formeln vgl. R. Sauer, Theoretische Einfiihrung in die Gasdynamik. 


erlin 1943; die der Flachwasserformeln verlauft analog. 
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2 Gi ALY pe Git 
ie 0, AQ’ ed 
1 qs 
P 
EO Ao ? W.., Vinge 2 
fiir die Schall- bzw. Grundwellengeschwindigkeit c 
a(&@ 
ieee Bene iso” 
— dg’ dé 
Ferner ist | 
Ap a= Oy ae Wn A554) = 6, Woy AM 
peal thee Cp wn) Tega 
2 ee Gp te (= ae 
” Ap= mi = sin? a, AS =) hoe Wi S sino. 


Die Bedeutung der Winkel « und *: ist der Abb. 2 zu entnehmen. 

Durch (4) und (5) ist der Wassersprung bei bekantem €, und ? vollkommen bestimmt. Dies isi 
anders beim Verdichtungssto8 im Gas. Hier hat man drei unbekannte Gréfen hinter dem StoB 
OX, 09, Py» Man braucht also fiir das Gas noch eine dritte Gleichung zur eindeutigen Festlegung des 
StoBes. Daim Wasser selbstverstandlich die triviale Identitat 


1 9 
got 


® : 1) =a 
; Toe. ea 
eon % 86a °2 


gilt, erhielte man dann vollige Ubereinstimmung zwischen Verdichtungssto8 und Wassersprung. 
wenn die entsprechende (hier nicht triviale) Beziehung auch beim Gas gilte, also die Adiabaten 


gleichung mit K = 2: 


Die Annahme der Giiltigkeit der Adiabatengleichung durch den VerdichtungsstoB hindurch ied 
aber physikalisch nicht zulassig. An die Stelle der Adiabatengleichung tritt die Energiegleichung 


B Le eaten oe Us 4} 

ay = Cree t. + 2g ° (6}) 
(Die entsprechende Gleichung fiir Wasser, also 

eRe as stalg Saeteee ae 

Cy 9 = Ge Ig 


kann nicht gelten, da sie im Widerspruch zu (1) bis (3) steht.) 

Alle Ergebnisse, die man unter Benutzung von (6) fiir den Verdichtungssto8 herleitet, unter- 
scheiden sich von den analogen fiir Wasser (erste Gruppe von Abweichungen). Da die Kenntnis 
der Abweichung von Bedeutung fiir die quantitative Auswertung der Experimente ist, werden die 
beiden wichtigsten voneinander abweichenden Beziehungen angegeben: 


CL oe Lee: Lae ee hay 
ENA ?) a pT AC et eo = 
Ap Berea (pee a ORL) Wea wale ee 
Op eae AC Ga se Sp a 
t91 = l993 
Oo1 = Op 5 Cyn = Gee Che 
02 01 4 Ges 
Por __ Po2 : 
Qo1 2o2 | 


Fir spatere Anwendung soll die letzte Gleichung umgeformt werden: Je Zeit- und Langeneinheit 
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n Richtung der Schwallfront) durchquert den Wassersprung die Wassermasse ! 


iese verliert hierbei die Energie 


so gilt 


ap a Bee (7) 


Energiedissipation je Zeit- und Langeneinheit des Wassersprungs). Wie sich die Abweichung 
wischen VerdichtungsstoB und Wassersprung in den StoSpolaren fiir beide Falle auBert, kann 
aan der oben zitierten Arbeit von Preiswerk entnehmen. 


4, Einflu8 von Bodenreibung und Oberflichenspannung. Wir beschrinken uns nun auf den 
snkrecht in ruhendes Wasser hincinlaufenden Schwall. Schrage Schwalle erhalt man (ebenso wie 
shrage VerdichtungsstiBe im Gas) wegen der Beziehung W,, = W,, durch einfache Uberlagerung 
iner Strémung parallel zur Schwallfront. : 

a) Bodenreibung. Zur Einfithrung der Bodenreibung schreibt man den Impulssatz fiir eine 
uaderférmige, von der Bodenflache bis zur Wasseroberflache reichende Fliissigkeitssiule von der 
Jinge dx in Fortschreitungsrichtung des Schwalls und der Breite 1 senkrecht dazu an: 


sefertn tlie ed hd 
0 

dier bedeutet p wie im folgenden den statischen Druck im Wasser, t die Bodenschubspannung 
Jie Kontinuitatsgleichung lautet 


Nach Ejnfiithrung der Mittelwerte 


c x 
Sy = = 
ia [ude und enya | 2 fara 
rhalt man unter der Annahme, daB man u == k U mit konstantem k setzen darf (dies gilt nur dann 
treng, wenn sich das Geschwindigkeitsprofil u(z) in Abhangigkeit von x affin verandert), und bei 
Jenutzung der Kontinuitatsgleichung in der Form 


) anaes 


ie Gleichung 
dé sates Ja | 
dx (1 3) me sea (8) 


ut der Abkiirzung a? = W?C?/g (cm*). Da ¢, die Tiefe des 
ngestorten, ruhenden Wassers bedeutet, hat W, die Be- Abb. 3. Wasserschwall bei Beriicksichtigung 
eutung der Schwallgeschwindigkeit. PPE ROR EME 

Macht man die selbstverstandliche Annahme, daB die Bodenschubspannung t der Relativge- 
shwindigkeit zwischen Wasser und Bodenflache entgegengerichtet ist, und mit dieser gegen Null 
eht, so ergeben sich aus (8) Lésungen von einem in Abb, 3 skizzierten Typ. Da die Lésung bei 
=a eine vertikale Tangente bekommt, und da nach Definition von a die Beziehung ¢, <a gilt 
1s ist namlich 6, < W2/g), folgt, da der Schwall auch bei (alleiniger) Beriicksichtigung der 
sodenreibung eine Unstetigkeit ausbilden muB. Diese mu dann in der aus Abb. 3 ersichtlichen 
Veise verlaufen, wie aus Anwendung des Impulssatzes auf das schraffierte Kontrollgebiet folgt. 
Vahlt man namlich dieses Gebiet geniigend schmal, so fallt die Bodenreibung ganz aus der Im- 
ulsbilanz heraus, die damit dasselbe Ergebnis wie im reibungsfreien Fall liefert. Man kann das 
irgebnis kurz so zusammenfassen: Die Bodenreibung allein kann das unstetige Schwallprofil 
icht abrunden und gerade deshalb andert sich an der Impulsbilanz und damit an den in Ab- 
shnitt 3 gegebenen Beziehungen nichts. Die Anderung gegenitber den reibungsfreien Verhalt- 
issen besteht in einem schwachen, stetigen Anstieg des Wasserspiegels hinter dem Schwall. 


4 o bedeutet von hier ab die Dichte des Wassers (auf die Gleichungen des Gases wird nicht mehr 
ariickgegriffen werden). 
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Allerdings muB beriicksichtigt werden, daB andere Einfliisse die Ausbildung einer unstetigen, 
Schwallfront verhindern kénnen, so da man fiir die Impulsbilanz das Kontrollgebiet nicht mehr 
unendlich schmal wahlen darf und somit auch die Bodenreibung beriicksichtigen mu. (Die Wasser+ 
hohen ¢, und €, sind bei einem Schwall mit stetigem Profil natiirlich weit genug vor und hinte 
dem Schwall zu nehmen.) Man kann sich aber davon iiberzeugen, daB die Bodenreibung so klein. 
ist, da® sie auch in diesen Fallen vernachlissigt werden kann. Um einen Uberblick iiber die Gréfe} 
der Bodenschubspannung und damit auch der Spiegelneigung hinter dem Schwall zu bekommen, 
wurde nach dem Pohlhausenverfahren die Grenzschicht berechnet, wie sie sich hinter einem in) 
ruhendes Wasser hineinlaufenden Schwall ausbildet. Fiihrt man die Schwallstarke k—(¢,/C, ein, 
so ergibt sich 


bei laminarer Grenzschicht: 


ie 1 3k + 2 
ae eh 
oW’ Re | 15 


(1 — k)t.15 (: i =e mit Re =— 


bei turbulenter Grenzschicht: 
see aan 0,023 
OWT VT Re 


(« = Entfernung von der Schwallfront; y = 0,01 cm?/s = kinematische Zahigkeit des Wassers).| 
Diese Formeln gelten eigentlich nur fiir Schwalle mit unstetigem Beginn, kénnen jedoch nahe-| 
rungsweise auch auf Schwalle mit stetigem Profil itibertragen werden. In allen praktisch inter-| 
essierenden Fallen kann man die mit diesen Schubspannungen aus (8) berechneten Spiegel- 
neigungen hinter dem Schwall vernachlassigen, der Wasserspiegel hinter dem Schwall ist prak- 
tisch waagerecht. 

b)Oberflachenspannung. Auch die Oberflachenspannung hat keinen wesentlichen 
Einflu8 auf die in Abschnitt 3 gegebenen Beziehungen. Berechnet man namlich die von der 
Oberflichenspannung T herrithrende horizontale Kraftkomponente auf das in Abb. 4 skizzierte| 
Kontrollgebiet, so ergibt sich 


Ce dé 
GUR, 2S —= TP 2 
yi+¢72° 
also 
a 1 wy 
= | GUN {iS | 
| aa yi 0? 


Da der Wasserspiegel vor dem Schwall waagerecht ist und hinter dem Schwall in guter An- 
naherung ebenfalls, wird K,—0. Die Oberflachenspannung geht daher nicht in die Impuls- 
bilanz ein. Es bedarf keiner besonderen Erérterung, daB auch die Vertikalgeschwindigkeit nicht 
in die Impulsbilanz eingeht. Vertikalgeschwindigkeit und Oberflachenspannung beeinflussen } 
jedoch stark die Form des Schwallprofils. 


5. Einflu8 der Vertikalgeschwindigkeit, Einzelwelle. Es bleiben nun noch die beiden wichtig- 
sten, eingangs aufgeworfenen Fragen zu beantworten iibrig, namlich: 

1. Wieso bleibt die Wasseroberflache in manchen Fallen zusammenhangend und warum zer- } 
reiBt sie manchmal ? 

2. Wie kann man die Entstehung einer stetigen Erhéhung vor dem Schwall in tiefem Wasser} 
und ihr Fehlen in flachem Wasser verstehen ? 

In diesem Abschnitt wird zunachst vorausgesetzt, daB die Wasseroberflache zusammenhangend 
bleibt, und unter dieser Voraussetzung Frage 2 beantwortet, waihrend die Bedingungen fiir die} 
Giltigkeit dieser Voraussetzung erst im nichsten Abschnitt erdrtert werden. — Die experimen- 
tellen Ergebnisse (Abschnitt 8) sowie friihere Untersuchungen der sogenannten Einzelwelle 
(= Welle, die aus einem einzigen Wellenberg besteht) haben nahegelegt, die Ursache der all- 
mihlichen, stetigen Erhéhung vor dem Schwall in der Vertikalgeschwindigkeit zu suchen und zu 
der Vorstellung gefiithrt, da diese Erhéhung mit guter Naherung dem ansteigenden Teil einer| 
Kinzelwelle bis zu deren Wendepunkt entspricht (vgl. Abb. 5). DaB die Wasseroberflache beim 
Schwall nicht wie bei der Einzelwelle nach Erreichen eines Maximums wieder absinkt, hat seine } 
Ursache darin, da der Schwall im Gegensatz zur Einzelwelle ein Energie dissipierender Vorgang 
ist. Man darf daher (und dem entspricht die soeben ausgesprochene Annahme) nur voraussetzen, 
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laf in einiger Entfernung vor dem Schwall die Verhiltnisse in bezug auf die Wirkung der Vertikal- 
eschwindigkeit denen im Vordergebiet einer Einzelwelle entsprechen, die mit derselben Geschwin- 
ligkeit wie der Schwall wandert. 

Eine Untersuchung der Eigenschaften einer Einzelwelle unter der Wirkung der Oberflichen- 
spannung des Wassers findet sich in einer Veréffentlichung von Harrison!. Da die Theorie dieser 
Welle von ausschlaggebender Bedeutung fiir die erwahnten Beobachtungstatsachen ist,-soll sie 
aier in einer von der Harrisonschen Darstellung abweichenden Form noch einmal gegeben werden. 
Ausgangspunkt sind die Eulerschen Gleichungen. Da nach einer weiter oben getroffenen Verein- 
barung nur zweidimensionale Vorginge (d.h. ebene Wellen) betrachtet werden, entfallt die 
Abhangigkeit von der y-Koordinate. (Uber die Bezeichnung der Koordinaten vgl. Abb. 1: 
Ww ist die Komponente der Geschwindigkeit in vertikaler Richtung.) Man hat 


du du 1 Op 
u -W — = — ; 
Ox ie 0 Ox’ (9a) 
Ow Ow 1 op 
uu + wi —_— ——— .--g; 9b 
Ox Es @ Oz 2 (2b) 
ae ae 0 Ow 

Kontinuitat: ioe Gi Rep ae (10a) 

Ox Oz 
; it i) Ow 

Wirbelfreiheit : id gl a5 =0, (10b) 

) z % ° 

Randbedingungen: w(z=0) =0, p(z=l)=—TC". (11a, b) 

Hierbei ist 1+ (’2?= 1 gesetzt.) 

! 

LLL LALLA hdd deddils 2 
Abb. 4. Einwirkung der Oberflichenspannung. Abb. 5. AnschluB® des Einzelwellengebietes an den Schwall. 
Man setzt wu und w als Potenzreihen in folgender Form an: 
“ge by u, (x) Zn : (12a) 
n=—0 


== > Wy (x) 2”. (12b) 
yl 


Setzt man nun in (10a) und (10b) ein, so ergibt sich durch Koeffizientenvergleich 


a OF 
dwn—1 
nu, = - De =a 9) 
(realy 2584050) 
dun—1 
nw, =— oy Ps 


ieraus kénnen die u, und w, durch die Ableitungen von uy ausgedriickt werden. Kinsetzen der 
u, und w, in die Reihen (12) ergibt 


2" dy z* d*ty Wee 
wm — ge age (18) 
dug Lid?uy ts set 13b 
“ihaqiden awd oo 


was iibrigens folgende Schreibweise zulabt: 


u—iw= {ett#dd=) uo(x), 


1 W. J. Harrison, Proc. Lond. Math. Soc. 7 (1909), S. 107. 


48 Becker: Analogie zwischen Wasserschwall und VerdichiungsstoB. Ingenieur-Archiw} 


wobei {...} einen auf u9(x) anzuwendenden Differentialoperator bedeutet. Mit (13a) und (13b) 
geht man in (9a) und (9b) ein. Beriicksichtigt man nur die Glieder bis 2°, so erhalt man dies 
Gleichungen von Abschnitt 2. Beriicksichtigung von Gliedern bis 2° liefert / 


du, , [1 d?ujduy dup pea eR. l4a\ 
a . dg ds. dx ON ae @ Ox’ ( a} 
dii,\ 2 dtugtly aot deta ries A. | 
(a) —atelm meee | 


Integration von (14b) itber z von z = 0 bis z =€ liefert bei Beachtung von (11b) 
1 [/du,\2 d@u I 

D) eae ee) Le reeeO) ss ss 47 ae % 

¢ x |( ae) aio] = gle + peo] —86 
Man nimmt nun an, da® |\¢ —C,| <¢,, wobei C, die Spiegelhéhe des ungestérten Wassers be- 
deutet, und erhalt durch Differentiation nach x 
C2 [d? uy d* uy duy ofl f Ponayy 
% E 0 dat dx 18s ae ' 


Dies setzt man in (14a) fiir 2 =0 ein und befriedigt auBerdem die Kontinuitatsgleichung » 
naherungsweise durch den Ansatz 
1 
Ug = (amen 
s =¢ —f,= Erhebung iber den ungestérten Wasserspiegel). Setzt man dann in d?u,/dx? und 
du,/dx? noch C,/(¢,+ s) ~ (1 — s/C,), so ergibt sich fiir s die Differentialgleichung 
a a® s/s’ 5! b2 3’ 
=l]4 . 15) 
G+) Raper Gr] 29 ( 
Orr 


Hierbei ist 6 = \/ av = 0,38 cm fiir Wasser. 


Bezeichnet man die Scheitelhéhe einer Einzelwelle mit ¢,, so mu sich diese Einzelwelle aus; 
(15) als diejenige Lésung ergeben, die folgende Randbedingungen befriedigt: 


er Ss!) fur? 25 ==.) (16a) 
s = 0 fiir sot (16b) 
Durch die Substitution s’?/2 = q(s), also s’’= dq/(s)/ds, s’’’/s’ = d? g(s)/ds? geht (15) in eine 


lineare Gleichung fiir y iiber. Fiir my kann man ansetzen: 


wobei es geniigt, von der Reihe die beiden ersten Glieder zu beriicksichtigen. Durch Einsetzen 
in die Differentialgleichung (15) fiir m findet man 


eral a? — C3 phen 1 a 
C2 a — BC,’ <a Ci a —b*6,” 


Aus (16b) folgt dann fiir den Parameter a 


as 


fe 4 
a= Tae F 
Nach der Definition von a ergibt sich hieraus die Wellengeschwindigkeit (mit der Vernach- 
lassigung |€,—¢,| <¢,) zu 
Way stir 
Das heiBt: die Geschwindigkeit der Einzelwelle ist unabhangig von der Oberflaichenspannung 


gleich der Grundwellengeschwindigkeit in ihrem Scheitelpunkt. Die Oberflichenspannung hat 
nur Einflu8 auf das Profil der Welle. Dieses wird durch die Gleichung 


Gherenle = (17) 


gegeben. Da nur reelles ds/dx in Frage kommt, muf die rechte Seite gréBer als Null sein. 
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FallI:f? < b?; in diesem Falle muB (C,/0,— 1) —s/C, < 0 sein, oder 6, < ¢, d. h.: wenn die 
Wassertiefe kleiner ist als h = 0,38 cm, ist nur ein Einzelwellental méglich; Auf dieselbe Art 
erkennt man: 

Fall Il: Wenn die Wassertiefe gréfer ist als 0,38 cm, ist nur ein Kinzelwellenberg méglich. 

Es wird hervorgehoben, daB ohne Beriicksichtigung der Vertikalgeschwindigkeit die Er- 
scheinung der Einzelwelle nicht zu verstehen ist. Dies gilt zum mindesten fiir Kinzelwellenberge, 
also in Wassertiefen > 0,38 cm. Fiir flaches Wasser haben weitergehende Untersuchungen nahe- 
gelegt, daB sich Einzelwellentiler auch unter alleiniger Wirkung der Oberflachenspannung ver- 
stehen lassen, wie iiberhaupt in Wasser so geringer Tiefe die Wirkungen der Vertikalgeschwindig- 
keit immer mehr gegeniiber denen der Oberflachenspannung vernachlassigt werden kénnen. 

Die Integration von (17) ergibt 


Smax 


s= 
¥ yf2 Smax (x = Xo)” 
2 es 
(% = Abszisse des Wellenscheitelpunktes). Diese Gleichung stimmt mit dem Ergebnis von 
Boussinesq und Rayleigh iiberein, die allerdings die Oberflachenspannung vernachlissigen und 
somit die Vertikalgeschwindigkeit allein beriicksichtigen, abgesehen von einem Faktor y3/2 im 
Argument des Go). Diese Diskrepanz liegt in der Verschiedenheit der Naherungsmethoden be- 
griindet. — Der Wendepunkt der Einzelwelle liegt immer in halber Scheitelhéhe. In der Nahe 
der kritischen Wassertiefe b sind die Einzelwellen sehr steil. Physikalische Uberlegungen an 
Hand der Dispersionskurven (Abb. 6) legen fiir die kritische Wassertiefe den Wert y3/2 b= 
0,5 cm nahe, wie er sich auch bei Benutzung der Rayleighschen Rechenmethode ergeben hatte. 
An dieser Stelle sei darauf hin- 
gewiesen, daf} man aus (15) bei ge- 


. - . 
| eigneter Abanderung der Randbedin- 
gungen (16) Lésungen erhilt, die 
ges . a5 
yee Wellen von endjicher Linge ent- 
iy! 
| 280 
“4, 
Yr aaa Ba 
22,3 
GSS 
G8 
ht 
Lo 
Z 
10, = brs eilt De eS one 
! a f e 3 Z é a coe Abb. 7. Profil von Wasserwellen fiir konstante 
b. 6. Dis sionsk fir Wasserwellen. Wellenlinge bei verschiedenen Wassertiefen 
pes Sanang Sioa phaktior®: (Voraussetzung: Wellenlange > 1,73 cm). 


| sprechen. Um einen Uberblick iitber den Giiltigkeitsbereich dieser Lésungen zu erhalten, wurde 
(17) durch Spezialisierung auf (im Grenzfall) unendlich Kleine Amplituden linearisiert. (Die 
Rechnungen werden hier iibergangen.) Man erhalt dann eine Beziehung zwischen Wellenlange A, 
| Wellengeschwindigkeit W und Wassertiefe C,, die mit der exakten, aus potentialtheoretischen 
‘Uberlegungen bekannten Dispersionsformel verglichen werden kann. Diese exakte Formel ist 


b? A Gs 
Wes s| ae za Ugens, 


die Naherungsformel 
Am? b2 4-27? 
4m? C2 + 272° 
In Abb. 6 sind beide Formeln graphisch dargestellt. Die Naherung ist fiir groBe Wassertiefen und 
kleine Wellenlangen verhaltnismaBig gut. Wellen mit dW/dd > 0 nennt man gewohbnlich Schwere- 
A 


W2 = gly 
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| 


wellen, mit dW/dA <0 Kapillarwellen. Es ist wichtig, hier zu bemerken, daB Oberflachenwellen | 
auf Wasser sich in ganz anderer Weise ausbreiten als Wellen in einem Gas, denn im Gas wandern 
die Wellen unabhangig von ihrer Wellenlange mit Schallgeschwindigkeit, wahrend bei Wasser- | 
wellen die Geschwindigkeit in der angegebenen Weise von der Wellenlinge abhingt. Dies hat 
u.a. zur Folge, daB man bei einer Flachwasserstrémung nicht in genau analoger Weise wie bei. 
der Gasstrémung durch Winkelmessung an einem von einer punktférmigen Stéiung ausgehenden | 
Wellensystem die Machsche Zahl der Strémung bestimmen kann, wie dies verschiedentlich ge- | 
schehen ist. 

Das Profil von Schwere- und Kapillarwellen unterscheidet sich in charakteristischer Weise,. 
wie Abb. 7 zeigt, wo das Wellenprofil fiir konstantbleibende Wellenlange in abnehmender Wasser- 
tiefe skizziert ist. Man erkennt, da® Schwerewellen steile Berge und flache, breite Taler auf- | 
weisen, Kapillarwellen umgekehrt. Wird die Wellenlange gréBer, so werden bei den Schwere: | 


wellen die Taler breiter 
und flacher, bis im) 
Grenzfall unendlich | 
groBer Wellenlange nur | 
noch ein Einzelwellen- | 
berg iibrig bleibt. Bei 
den Kapillarwellen bleibt | 
im gleichen Grenzfall | 
das Einzelwellental. - 
Hierzu ist  allerdings | 
Voraussetzung, da die 
Wassertiefe, kleiner als 
etwa 0,4cm ist, denn 
nur dann gilt unabhan- | 
gig von der Wellenlange 
dW/di< 0, so daB auch | 
lange Wellen Kapillar- 
wellen sind. Es ist nun 
klar, warum sich in fla- 
chem Wasser (< 0,38cm) 
vor einem Schwall auch 
cae: Pe 22’ bei stetig bleibender 

G BIE OBOE. RT ORES se ee Wasser bi fliche keine 
allmahliche, stetige Er- 
héhung ausbilden kann, 
denn nach obigen Ergebnissen ware in solchen Wassertiefen allenfalls eine stetige Absenkung 
zu erwarten. Daf eine solche nicht beobachtet wird, liegt daran, dafS§ ein Einzelwellental 
langsamer wandert als es der Grundwellengeschwindigkeit des ungestérten Wassers ent- 
spricht, ein Schwall dagegen schneller. Auferdem kann man im Fall des Einzelwellentals | 
keine verniinftige Annahme mehr machen iiber den Anschlu8 des Einzelwellengebietes an den 
Schwall. 

Geht man von der zu Beginn des Abschnitts gemachten Annahme iiber den AnschluB des | 
Kinzelwellengebietes an den Schwall aus, so ergibt sich fiir die maximale allmahliche Erhebung © 
vor dem Schwall (Wassertiefe > 0,38 cm) 


eck Smax a ide tine 
a mleti 1). (18a) 


Vernachlassigt man die Oberflachenspannung, so berechnet man, daB diese stetige Erhéhung in 
einer Entfernung 


If Maximale stenige Erhohung 
vor der Schwalitront 


Lange der sterigen Lrhohung 
vor der Schwolitront 


Abb. 8. Maximale stetige Erhéhung 1/2 s,,4, vor der Schwallfront und Linge Ax der stetigen 
Erhéhung bis zum Abklingen auf 1/10 als Funktion der Wassersprunghéhe ¢,— ¢,. 


Ax a ar (18b) 


auf 1/10 ihres urspriinglichen Wertes abgeklungen ist. In Wassertiefen die kleiner als etwa 
0,5 cm sind, mu8 man die Oberflachenspannung beriicksichtigen. Die Formeln (18a) und (18b) 
_ verlieren dann ihren Sinn. In Abb. 8 sind 1/2 s,,,, und 4x als Funktionen der Wasserschwallhéhe 
mit C, als Parameter aufgetragen. 
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6. Kapillarwellen. Es bleibt nun noch die zu Beginn des vorigen Abschnittes zuerst auf- 
geworfenen Frage zu lésen, die in folgender Weise prazisiert werden kann: Wie kann im 
Schwall Energie dissipiert werden ohne Bildung einer Wasserwalze (= turbulente Durch- 
mischungszone) ? Die Beantwortung ergibt sich aus der Tatsache, da® vor allen beobachteten 
Wasserschwallen ein System von Kapillarwellen mit zur Schwallfront paralleler Wellenfront 
herlauft. Diese Kapillarwellen iiberlagern sich erfahrungsgema8 der allmihlichen stetigen Er- 
héhung (falls diese vorhanden) ohne sie wesentlich zu stéren. Die Phasengeschwindigkeit der 
Wellen ist der Schwallgeschwindigkeit gleich und somit immer gri®er als die Grundwellen- 
geschwindigkeit. Ein Blick auf Abb. 6 ergibt dann, daB die der Schwallfiont vorgelagerten 
Wellen sich stets im Gebiet dW/dd < 0 befinden, daB also ihre Gruppengeschwindigkeit gréBer 
ist als ihre Phasengeschwindigkeit. Hieraus folgt, daB der Kopf 
der Wellengruppe schneller wandert als die Phase und somit dem 
Wasserschwall vorauseilt, wie es Abb. 9 zu verschiedenen Zeiten 
zeigt. Dabei sind Reibungseinfliisse zunadchst auBer Betracht ge- 
lassen. Dem Wellensystem muB stindig Energie zugefiihrt werden, 
da es sich mit der Zeit iiber eine immer gréBer werdende Wasser- 
flache ausbreitet. Diese Energie wird vom Schwall geliefert. 
Bei starker werdendem Schwall miissen die Amplituden der 
Kapillarwellen immer gréBer werden, wenn das Wellensystem die 
ganze vom Schwall kommende Energie (oder doch den gréften 
Teil) noch aufnehmen soll. ErfahrungsgemaB® sind diese Ampli- 
tuden aber beschrankt, und man sieht daher ein, daB von einer 
bestimmten Schwallstéike an die Energie auf andere Weise ver- 
zehrt werden mu, namlich in einer Wasserwalze. 

Auf folgende Weise kann man ausgehend von diesen Vor- 
stellungen die Kapillarwellenamplitude sy) berechnen: Man setzt 


die Kapillarwellen als harmonische Schwingungen an (nur fiir ABbeO adettlichensbildane 
: : eerie ee des Kapillarwellensystems vor 
sehr kleine Amplituden streng giiltig): seep ese aly 


s=s, cos(kx—wt). 
Bei Beschrankung auf Tiefwasserwellen (fiir € > 0,4cem, A < 0,6 cm sicher zulassig) gilt fiir die 


Gruppengeschwindigkeit? W=3 /2 W. Die Energie unseres Wellensystems wird fiir eine Flache 
von der Lange | in x-Richtung und der Breite 1 senkrecht dazu ? 


1 
E50 kW? s2. 
Je Langeneinheit in y-Richtung wird daher in der Zeiteinheit die Energie 
1 ee 1 
a= Fe kW? s2(W — W) ira kW? s 


abtransportiert. Diesen Ausdruck vergleicht man nun mit (7) und erhalt 


& (A¢)8 : (19) 


Hierbei wurde ¢ < 1 eingefithrt, um der Tatsache Rechnung zu tragen, daf ein Teil der Energie 
ja wohl schon im eigentlichen Schwall vernichtet wird. Experimente haben ergeben, 
0,6 <e <0,9; die aus (19) berechneten Amplituden stimmen daun gut mit der Beobachtung 
iberein. 

Diese Ubereinstimmung bleibt auch dann erhalten, wenn man die innere Reibung beriick- 
sichtigt, die sich natiirlich darin duBert, daB die Wellenamplitude mit groBer werdender Ent- 
fernung vom Schwall abklingt. Ein bestimmter Teil der vom Schwall in das Wellensystem 


| hineingelieferten Energie wird namlich von der ersten Welle in Warme umgewandelt und nur 


der Rest gelangt zur zweiten, in der sich der Vorgang wiederholt. Unter dieser Annahme ergibt 
sich ein Dampfungsgesetz von der Form 


le 


a 
» 


So = So @ 


1 A, Sommerfeld, Vorlesungen iiber theoretische Physik. Bd. II, S. 180. Wiesbaden 1947. 
2 H,. Lamb, hore: der Hydrodynamik. S. 707 2, dt. Auflage. Leipzig und Berlin 1931. 


4* 
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Die Gréfbe B kann auf folgende Weise berechnet werden. Man nimmt an, dafs schon so lange Zeit | 


seit Entstehung des Schwalls vergangen ist, da sich der Kopf des Wellensystems in einem Ge- 
biet befindet, in dem die Amplituden praktisch auf Null abgeklungen sind. Dann wird die ganze 
in das Wellensystem hineingeschickte Energie dort in Warme verwandelt. Nach dem oben 


zitierten Lehrbuch von Lamb wird von einem Kapillarwellensystem von der Lange dx in Fort- | 
schreitungsrichtung und der Breite 1 senkrecht dazu folgende Energie je Zeiteinheit dissipiert: | 


i a dx =21,18 WS dx. 


Die gesamte dissipierte Energie wird dann 


es) 2% 
—_— 3 
i 2H WAS f dee 7 whe gS ep. 
2 
0) 


8 Durch Vergleich mit (7) erhalt 


man 


Wa 


Fortsohreitungsrichtung a 162?» 


(u = Zahigkeit, 
y =kinematische Zahigkeit). 
Die so berechneten f stimmen 


voll ausgebildere 
Linze/welle 


iiberein, wenn man die Ober- 


SC ——Messung —= etwas kleiner annimmt als dem 
\’ ——hechn nachGl17 Tabellenwert entspricht. Diese 
. \<-— fecha. nachhayleigh Verkleinerung 1aBt sich leicht | 
SS auf Verunreinigung der Wasser- 
Pe oberflache zurickfiihren +. 
Einzelwele hure nach 7. Experimentelle Ergeb- | 
ER ely nisse. Zur Bestatigung der 
=a, 0 are, 2 ome diskutierten Vorstellungen 
werden die Profile einer groBen 
Abb. 10. Vergleich zweier beobachteter Einzelwellenprofile mit der Rechnung Talli von Schwallen, Einzela 


von Rayleigh und der Rechnung nach Gl. (17). 
wellen und Wellensystemen bei 


Wassertiefen von 0,2 bis 1 em vermessen. Im allgemeinen ergab sich gute Ubereinstimmung mit 
der Theorie. Die Schwalle und Wellen wurden entweder dadurch erzeugt, daB der Wasserspiegel 
plétzlich durch aus einem Tank ausflieBendes Wasser erhéht wurde, oder indem eine flache Leiste 
oder eine gegen die Bewegungsrichtung angestellte Platte durch das Wasser geschleppt wurde. Man 
erhalt ein genaues Bild des Wellenprofils, indem man parallel zur Wasseroberflache einen diinnen, 
von einer starken Lichtquelle angestrahlten Draht spannt und diesen und sein Spiegelbild in der 
durch die Wellen deformierten Wasseroberflaiche photographiert. An Stelle des Drahtes kann man 


auch den genau gradlinigen unteren Rand einer gut reflektierenden Flache benutzen. Gemessen — 


wurden prinzipiell nur Langen. Wenn zum Vergleich mit der Theorie die Kenntnis von Ge- 
schwindigkeiten nétig war, so wurden diese aus Héhendifferenzen berechnet. 

In Abb. 10 ist das gemessene Profil einer Einzelwelle mit dem Ergebnis von Rayleigh und 
Boussinesq (nach Korrektur durch die Oberflachenspannung) verglichen. Auferdem ist das 
gema (17) berechnete Wellenprofil eingezeichnet. In der ersten Zeit nach der Entstehung einer 
Einzelwelle ist ihr ein deutlich sichtbares System von Kapillarwellen itberlagert. Diese stéren 
jedoch das Wellenprofil nicht wesentlich und verschwinden bald durch innere Reibung. (Im 
Gegensatz zum Schwall liefert eine Kinzelwelle ja keine Energie nach.) In Wassertiefen < 0,3 em 
wurden in Ubereinstimmung mit der Theorie keine Einzelwellenberge beobachtet. Da® auch 
keine gut ausgebildeten Einzelwellentiler beobachtet werden konnten, liegt daran, da® diese nie 
ganz von den unvermeidlichen Kapillarwellen endlicher Lange getrennt werden konnten und 
ebenso schnell wie diese durch Dampfung verschwinden. Abb. lle ist ein Beispiel fiir eine dem 
Kinzelwellental nahekommende Erscheinung. Die Abb. 11a bis lle zeigen eine Anzahl von Einzel- 


' G. Bakker, Kapillaritat und Oberflachenspannung. Hb. Exp. Phys. Bd. VI, S. 307. Leipzig 1928. 


mit den beobachteten gut | 


flachenspannung des Wassers | 
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wellen, die Abb. 12a bis 12k mehrere Schwalle bzw. Ubergangsformen zwischen Schwallen und 
Kinzelwellen. In den Fallen, in denen nach der Theorie eine allmahliche, stetige EKrhéhung zu 
erwarten war, ist diese gemaf den Gleichungen des Abschnitt 5 eingezeichnet. Bei allen Schwallen 
wurde fernerhin die vom Kapillarwellensystem dissipierte Energie berechnet. Wenn das Kapillar- 
wellensystem gut ausgebildet ist, werden inihm 
stets tiber 60% der gesamten Energie in Warme 
verwandelt. Man erkennt an den Bildern gut 


a 
die charakteristische Profilform der Kapillar- Pn for Mert i) aca 
wellen. Nach den Beobachtungen ist die Bil- | ides 
dung einer Wasserwalze beim Schwall dann zu é ee 
erwarten, wenn A/s)< 5 wird. Se ae 
_ 8. Zusammenfassung. 1. Die Analogie | 
zwischen Wassersprung und VerdichtungsstoB = {0 ~ i Pe 


ist schon im Rahmen der ,,gasdynamischen 
Naherung* der Flachwassergleichungen nicht 
vollkommen, da das thermodynamische Ver- |. GE 
halten in beiden Fallen verschieden ist (Ab- 
See ei aman e wearin er TRE STP) PSL hiner ete i a d 
2. Die in Abschnitt 3 gegebenen Beziehungen 
zwischen der Geschwindigkeit eines Schwalls 
und den Wasserhéhen vor und hinter dem [77777 7 
Schwall werden durch Bodenreibung, Ober- € 
flachenspannung und Vertikalgeschwindigkeit 
‘nicht oder doch nur sehr unwesentlich ge- 


7 hd 
f 
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i Il schiedener Starke Abb. 12. Beobachtete Schwalle verschiedener Starke und 
“ca aa De alison f.. Ubergangsformen zwischen Schwallen und Einzelwellen 


in verschiedenen Wassertiefen ¢,. 


andert (Abschnitt 4). Die Abweichungen gegenitber der Gasstrémung bestehen vielmehr darin, 
da die Profilform des Wasserschwalls im allgemeinen keine exakte Unstetigkeit aufweist. 

A. Sehr flaches Wasser: (<8. . 

3 Schwache Wasserschwalle: Die Oberflichenspannung verhindert ein Zerreiben 
der Wasseroberflache, und es kommt zur Bildung von Kapillarwellen. Diese Kapillarwellen iiber- 
nehmen die Energiedissipation des Schwalls. Eine stetige Frhéhung v or dem Schwall tritt nicht 
auf, da in so flachem Wasser allenfalls eine stetige Absenkung auftreten kénnte. Daf diese auch 
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nicht beobachtet wird, riihrt daher, da® der Schwall mit Ubergrundwellengeschwindigkeit | 
wandert, das Einzelwellental, das einer stetigen Absenkung entsprechen wiirde, aber mit 
Untergrundwellengeschwindigkeit (Abschnitt 5). 

4. Starke Wasserschwalle: Wenn die nach (19) berechneten Amplituden der _ 
Kapillarwellen zu gro® werden, bricht sich der Wasserschwall. Es kommt zur Ausbildung einer | 
Wasserwalze, die an Stelle der Kapillarwellen die Energiedissipation iitbernimmt. Infolge der 
értlichen starken Stérung der Wasseroberflache im Gebiet der Wasserwalze treten zwar auch | 
dann noch Kapillarwellen vor dem Schwall auf; ihre Amplituden sind aber sehr klein. 


BB. Tieferes Wasser: ¢, > B: 

5. Schwache Wasserschwalle: Es kommt wie in Fall 3 zur Bildung von Kapillar- | 
wellen vor dem Schwall, die die Dissipation der Energie tibernehmen. Die Ausbildung einer | 
Wasserwalze erfolgt unter diesen Umstanden nicht, die Wasseroberflache bleibt stetig. In diesem | 
Fall verursacht die Vertikalgeschwindigkeit eine stetige Exhéhung vor der Sto®front, deren | 
Gré8e mit guter Annaherung durch (18a) und (18b) gegeben ist. Die Kapillarwellen sind dieser 
Erhéhung iiberlagert. 

6. Starke Wasserschwalle: Sobald die Amplituden der Kapillarwellen gemaf (19) 
zu gro werden, setzt die Bildung von Wasserwalzen ein. Damit verschwinden wie in Fall 4 | 
die Kapillarwellen bis auf einen unbedeutenden Rest. Gleichzeitig verschwindet auch die stetige | 
Erhéhung: diese wird von der Wasserwalze iiberdeckt. 

Vom Standpunkt der Analogie mit der Gasdynamik ist die Bildung einer Wasserwalze 
wiinschenswert, da sie am ehesten einem VerdichtungsstoB entspricht und da die Wasser- | 
oberflache vor dem Schwall nahezu ungestért und damit eben bleibt. Man miifBte also im Ex- | 
periment die Wassertiefe und die Geschwindigkeit der durch das Wasser geschleppten Modelle 
so groB machen, dab die entstehenden Wasserschwalle eine Wasserwalze ausbilden. Dies hat den 
Nachteil, da die experimentelle Apparatur ziemlich gro wiirde. Doch darf man dessen un- 
geachtet feststellen, daB Wassertiefen von weniger als 0,5 cm schlecht geeignet sind zu einem 
genauen quantitativen Vergleich mit den gasdynamischen Formeln. Am giinstigsten sind wohl | 
Wassertiefen zwischen 0,5 cm und 2 cm. 

Diese Arbeit wurde 1951 im Max-Planck-Institut fiir Strémungsforschung zu Gottingen 
durchgefithrt. Herrn Prof. Betz sei an dieser Stelle fiir Anregung und stete Unterstiitzung der 
Arbeit gedankt. 


(Eingegangen am 24, Juni 1952.) 


Anschrift des Verfassers: E. Becker, Gottingen, Bottingerstr. 6/8, 


Max-Planck-Institut fiir Str6mungsforschung. 
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Uber die Berechnung der Hypoidgetriebe. 
(I. Mitteilung.) 
Von P, Matthieu. 


1, Einleitung. Wahrend die Theorie der Stirn- und Kegelradgetriebe langst ausgebaut und 


heute weitgehend abgeschlossen ist, ist diejenige der Hypoidgetriebe (auch Hyperboloidgetriebe 


_oder Schraubgetriebe genannt) noch verhaltnismafig wenig entwickelt. Es liegt das daran, daf 
die Theorie der Hypoidgetriebe trotz des einfachen mechanischen Sachverhalts auf betrachtliche 


| 


] 


mathematische Schwierigkeiten fiihrt. 

Der erste, der sich mit Hypoidverzahnungen befafte, war Th. Olivier 1. Er stellte zu Beginn 
des vergangenen Jahrhunderts gewisse mathematische Begriffsbildungen zur Betrachtung be- 
liebiger Zahnradgetriebe auf und fand, von diesen ausgehend, eine Reihe von Hypoidgetrieben 
mit Punktverzahnung sowie auch ein solches mit Linienverzahnung. Diese Untersuchungen 
Oliviers blieben lange Zeit die einzigen auf diesem Gebiete. Auch der Altmeister der Getriebe- 


~lehre in Deutschland, F. Reuleaux, befaBts sich nicht wesentlich mit Hypoidgetrieben. Nach 
einigen kiirzeren Untersuchungen verschiedenartiger Natur ? setzte eigentlich erst zu Beginn 


dieses Jahrhunderts eine systematische Weiterentwicklung der Theorie der Hypoidgetriebe ein. 
Der Eigenartigkeit des Problems entsprechend erstreckte sich diese vor allem in zwei Richtungen. 
Kinerseits lag es nahe, auch im jetzigen, raumlichen Falle nach einfachen speziellen Verzahnungs- 
formen zu fragen, nachdem man bei den ebenen Zahnriadern in Gestalt der Zykloiden-, Elvol- 
venten- und Triebstockverzahnung einfache geometrische Formen gefunden hatte, die mit- 
einander im Eingriff stehen. Aber dieser Weg erwies sich jetzt als bedeutend schwieriger. Es 
lassen sich zwar ohne Miihe die gleichen kinematischen Grundbegriffe wie im ebenen Fall 
(Achsoide, Hiillflachenbeziehungen usw.) einfiihren*®. Von diesen ausgehend lassen sich auch 
viele spezielle rdumliche Zahnformen finden, vor allem solche, die die Gestalt von Regelflachen 
haben 4. Jedoch erreicht diese Theorie niemals die Geschlossenheit des ebenen Falles. Auch ist 
sie fiir praktische Zwecke kaum verwendbar. 

Aus diesem Grunde hat sich, vor allem infolge des Aufkommens der Hypoidrader in der 
Praxis, neben dieser ersten Richtung cine zweite, mehr den Anwendungen entsprechende heraus- 
gebildet. Sie kiimmert sich im Gegensatz zur ersten weniger um theoretisch zwar interessante, 
praktisch jedoch kaum verwendbare Eingriffsformen. Vielmehr beschaftigt sie sich mit Zahnrad- 
flankenformen und -problemen (unter Umstanden auch nur in geniigender Approximation) wie sie 
in der Praxis auftreten. Hier sind vor allem die Arbeiten von E. Wildhaber ® und J. Capelle ® zu 
nennen. Diese zweite Richtung steht noch in ihren Anfangen, jedoch ist zu erwarten, daB sie 
der allgemeinen Tendenz der angewandten Mathematik entsprechend, noch zu bedeutsamen 
Resultaten fiihren wird. 

Die vorliegende Arbeit schlie&t sich zunachst an die erste dieser Richtungen an. Sie be- 
schaftigt sich mit dem urspriinglichen Problem, zu einer gegebenen Zahnflanke die mit ihr im 
Eingriff stehende Gegenflanke rechnerisch exakt zu bestimmen. Bekanntlich ist die Lésung 
dieser Aufgabe bei gegebenen Achsen und Rotationsgeschwindigkeiten stets eindeutig bestimmt. 
Bezogen auf das System der Gegenflanke fiihrt die gegebene Flanke eine Relativbewegung aus, 
deren Hiillflache die Gegenflanke ist. Von dieser Erzeugungsweise ausgehend, 1a8t sich die 


1 Er hat seine Resultate zusammengefaBt in dem Werk ,,Théorie géométrique des engrenages". 
Paris 1842. 

2 Diese sind angegeben in der Enzyklopadie der mathematischen Wissenschaften, Band IV (Me- 
chanik), Heft 3 (Kinematik, von A. Schénflies und M. Griibler), Nr. 27 (Zahnrader und andere Mechanis- 
men), S. 261 ff. ; sale 

3 Vel. z. B. das unter ? genannte Werk. Le 

2 Vel. die folgenden Arbeiten: M. Disteli, Z. Math, Phys. 53 (1904), S.51. M, Disteli, Z. Math. 
Phys. 60 (1911), S. 244. E. Stiibler, Z. angew. Math. Mech. 2 (1922), $.429. R. Crain, Werkstatts- 
technik 1 (1907), S. 88. dpe 

> a, Wildhaber, Basic Relationschip of Hypoid Gears. The Machinist (englische Ausgabe), 
Jahrgang 1946, Hefte vom 1. Juni, 8. Juni, 22. Juni, 28. Sept., 26. Okt., 23. Nov., 80. Nov., 14. Dez. 

6 J, Capelle, Théorie et calcul des engrenages hypoids. Edition Dunod. Paris 1949, 
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Gegenflanke nach bekannten Prinzipien leicht durch ein Gleichungssystem charakterisieren, wie} 
das z. B. H. Resal getan hat}. 

Offenbar wegen der zuniichst sehr kompliziert scheinenden Gestalt dieser Gleichungen sind 
aber bis jetzt nie Methoden zu ihrer Lésung entwickelt worden. Demgegeniiber soll in der vor-| 
liegenden Arbeit gezeigt werden, da sich das Gleichungssystem fiir die Gegenflanke in tiber- |i 
raschender Weise sogar in endlicher Form lésen laBt. Dadurch wird die Theorie der Hypoidrader |} 
auf eine neue Basis gestellt. Einerseits }/) 
kénnen die einfachen Gleichungen, die | 
zwischen einer Flanke und ihrer Gegen- | 
flanke bestehen, dazu benutzt werden, | 
um miglichst einfache Lisungen des | 
Verzahnungsproblems in analythischer 
Form zu finden. Andererseits kann man 
bei beliebig gegebener Flanke die Gegen- 
flanke stets punktweise numerisch be- 
rechnen, so daf man gleichzeitig auch 
eine einfache, fiir alle Falle giiltige, 
numerische Berechnungsmethode hat. 
AnschlieBend an diese exakte Lésung 
kénnen dann weiterhin auch noch die 
Abb. 1. Schematische Darstellung praktischen Methoden der erwahnten | 

Gos eae he zweiten Richtung zur Anwendung ge- | 
langen. 

In der vorliegenden Mitteilung wer- 
den die Grundgleichungen des Ver- 
zahnungsproblems aufgestellt und ge- 
list fiir den einfachsten und weitaus 
wichtigsten Fall von sich rechtwinklig 

kreuzenden Achsen. Der Fall beliebiger windschiefer Achsen sowie einige sich anschlieBend er- 
gebender Probleme werden in einer Reihe weiterer Mitteilungen behandelt werden. 


2. Aufstellung der Gleichungen. Wir legen das Getriebe auf die in Abb. 1 schematisch in 
Grund-, Auf- und Seitenri®B angegebene Art in einen ruhenden Raum bzw. in ein rechtwinkliges 
Koordinatensystem mit den Koordinaten x, y, z. Die Achse des Tellerrades fallt mit der x-Achse 


Abb. 2, Anschauliche Darstellung des Tellerradsystems. Abb. 3. Anschauliche Darstellung des Ritzelsystems. 


zusammen, diejenige des Ritzels verlauft im Abstand E parallel zur y-Achse. Die Ausdriicke 
»»Tellerrad’* und ,,Ritzel‘* sind nur der Bequemlichkeit wegen eingefiihrt und haben keine tiefere 
Bedeutung, insbesondere nicht in bezug auf die GréBenverhaltnisse. In iiblicher Weise werden 
alle GréBen, die sich auf das Ritzel beziehen, mit dem Index 1, diejenigen die sich auf das Teller- 
rad beziehen, mit dem Index 2 versehen. Die Winkelgeschwindigkeiten der beiden Rader seien 
mit w, und w, bezeichnet. w,, w, und E kénnen beliebig gegeben sein. 


* Comptes rendus de l’ac. des se 117 (1893), 
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Weiterhin fithren wir einen mit dem Tellerrad fest verbundenen Raum ein (der dann natur- 
gema gegeniiber dem ersten, ruhenden Raum beweglich ist) und in diesem die rechtwinkligen 
Koordinaten %,, y,, z, (Abb. 2), ebenso einen mit dem Ritzel fest verbundenen Raum mit den 
Koordinaten x,, y,, 2, (Abb. 3). Das ruhende Koordenatensystem und das des Tellerrades fallen 
zur Zeit t = 0 vollstaéndig zusammen, wahrend dasjenige des Ritzels ihnen gegeniiber um den 
Betrag E in der 2-Kichtung parallel verschoben ist. Nachher verdrehen sich die Systeme gegen- 
einander. Statt rechtwinkligen Koordinaten kann man im Ritzel- und im Tellerradraum auch 
Polarkoordinaten r,, 1, py bzw. ry, %, ~, einfiihren (Abb. 4u. 3). 

Der Ritzel-, der Tellerrad- und der ruhende Raum bzw. die in ihnen liegenden Koordinaten- 
systeme fiihren gegeneinander Relativbewegungen aus. Fiir unsere Zwecke wird vor allem die- 
jenige zwischen den beiden ersten Raumen von Bedeutung sein. Die Koordinaten der fiinf ein- 


22 


Abb. 4. Einfiihrung von Polarkoordinaten im Ritzelsystem. Abb. 5, Einfithrung von Polarkoordinaten im Tellerradsystem. 


gefiihrten Systeme lassen sich leicht aufeinander umrechnen. Zunichst bestehen zwischen den 
rechtwinkligen Koordinaten des ruhenden und des Ritzelraumes die Beziehungen 


x, = xcos ,t + (z — E): sin a,t, x = %, cos W,t — z, sin wt, | 
Ya oo N= Yivs g) 
% = — *xsin w,t + (2— E)cosw,t, «= E+ x sin wt +2, cos a,t 


und ebenso zwischen denjenigen des ruhenden und des Tellerradraumes die Gleichungen 


Xo ES a Xo: 
Yo = ycos w,t +2zsin w,t, ¥ = Y_ C08 Wat — 2 sin Wet, (2) 
Zo == — ¥ 81n Wt 4-2c08s Wt, 7 = Vo sitvW,t + % cos Wet. | 


Daraus folgen fiir die beiden bewegten rechtwinkligen Systeme die Umrechnungsformeln 
XM, = %, cos Wt — 2, sin w,t, | 
Yo = ¥1 COS Wot + x, sin Wt sin Wyt + % cos Mt sin wt + E sin wyt , (3) 
Zy = — ¥1, Sin Wet + X, sin W,t COs Wet -+ 2% COs W, t cos Wat +E cos w,t, | 
X%, = Y. Sin W,tsin Wet + x, cos W,t + 2, sin W,t cos Wt — Esin wt, | 
V1 = Yo COS Wt — %, SIN Wet, (4) 
Z, == Yo COS Wt Sin Wot — X sin W,t + % cos W,t cos Wt — Ecos w,t. | 


Hier kénnte man bereits die eine Zahnflanke, z. B. die des Tellerrades, als Flache in der 
Gestalt 


Jo (%2» Yo» %2) = 9 (5) 
beliebig annehmen und diese Gleichung gemaf (3) auf den Ritzelraum umrechnen, wodurch sich 
in diesem eine von t abhangige, also bewegte Flache in der Gestalt 


Si (1: Nv A 1) = 9 (6) 
ergibt. Die Gegenflanke ist dann die Hiillflache dieser bewegten Flache. Mathematisch bestimmt 
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sie sich nach bekannten Regeln dadurch, daB man (6) partiell nach t differentiiert, also die 
Gleichung 
of, 


Ot (1 Yr %1» t) = 9 (7) 


bildet, und t aus (6) und (7) eliminiert. 

Dieser prinzipiell einfache Weg fiihrt jedoch in rechtwinkligen Koordinaten auf grofe | 
rechnerische Schwierigkeiten. Wesentlich einfacher werden die Gleichungen, wenn man die 
oben eingefiihrten Polarkoordinaten verwendet. Als Umrechnungsformeln ergeben sich fiir die _ 
beiden Koordinatensysteme des Ritzelraumes gemaB Abb. 4 die Gleichungen 


es 1 9. 6 2 Pall 2 
x, =17,sinB,cosy,, 4 =le+yt+e2, 


Ren 
¥, = 7,008 B, , i, == are ya ; (8) | 


: . zy 
ZS 1 sin P, sing, —, = arete mae 
f 1 


Fiir den Tellerradraum dagegen hat man die Beziehungen 


Se Serres 2 yee 
NF COS Wy n= x2 + ¥g + 295 


; x 
ve == 7, S10 0, COS Cy, 0, care cop mas ; | (9) 
2 


My eat pect eee 
Z = T, SIN Vy SINP., Py = are bh. 


Daraus folgen weiterhin gemaf (1) und (2) zwischen den beiden Polarkoordinatensystemen einer- 

seits und den rechtwinkligen Koordinaten des ruhenden Raumes andererseits die folgenden 

Relationen, die sich auch nach Abb. 4 und 5 sofort aus der Anschauung gewinnen lassen: 
x =r,sin J, cos (y, + 0,1) = 1 cos D, , | 


y = 1, cos 0, = 7, sin Y, cos (Y,-+ Wet), (10) 


z= 71, sin 0, sin (gy, + w,t) + E = rsin J, sin (py + Wet). | 
Nach diesen Formeln lassen sich die beiden Polarkoordinatensysteme und die rechtwinkligen 
Koordinaten des ruhenden Systems leicht aufeinander umrechnen. Speziell ergeben sich fiir die 
beiden Polarkoordinatensysteme die Umrechnungsformeln: s 
rT, = {E +72 +2 Er, sin J, sin (y, + @,¢), 

Gare cee ie sin 0, cos (py + a ? 


= (11) 
PY, = — w,t + are tg fe - ae = °) j 

1, = /E? +r —2 Er, sin 8, sin (p, + a!) , 
9, = are cos ("207 Pacoe (oe Fos), (12) 


Y, = —,t + are tg Gas O,8in (Py + st) =) ' 


r, cos O, 


Wir wenden uns nun der oben im Prinzip bereits durchgefihrten Berechnung der Gegen- 
flanke als Hiillflache zu. Die eine Zahnflanke, z. B. die des Tellerrades, sei wieder gegeben. 
Wir nehmen sie jetzt in der Form 


P2 = Py (12, V2) (13) 
an’. Indem wir hier fiir p,, 0, und r, die Werte (11) einsetzen, erhalten wir die bewegte Tellerrad- 
* Der Fall, daB die Flanke in der allgemeinen Gestalt D* (Ay, Oy, Po) = 0 gegeben ist, laBt sich, wie man 


aus den folgenden Ausfiihrungen leicht erkennen wird, sofort auf den obigen Fall zuriickfiihren. Wir halten 
uns deshalb im folgenden an die Darstellung (13). 
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flanke im Ritzelsystem in der Gestalt 


ery Rk 
Coiaeree ( ts me ae + @,t) = D,(r,, 9.) , (14) 


wobei zur Abkiirzung die GréBen 


0, ae, 2 sin J, cos (p, + il (15) 
Up 
beibehalten wurden. Die partielle Ableitung von (14) nach t ergibt 
, rT, sin @, cos (p,+ «@,t) ) 
rT, cos B, _ 0®, Onn OD> 00, 
W, 4 14 E +7, sin 0, sin (p, + @,t) P + re ("a Ps) Ot ny 00, (To, Ya) a 
r, cos B, 
o@ Er,sin 8,¢ t 
= 04 F* (Fay By) SEES SOs a Ot (16) 
oR? (ran Os) aes . 
BS wo 00, _[r, sin 8 sin (p, + @,¢) ! Erj sin? 8, cos? (y, + 4) 
z ie Ws ¢ sin 8, cos (~; + @,¢) i 19 re 
re 


Die Hiillflache, also die gesuchte Zahnflanke des Ritzels erhalt man nun wieder durch Elimination 
von t aus (14) und (16) unter Beriicksichtigung der Abkiirzungen (15). Die Gleichungen (14) 
bis (16) stellen also das gesuchte Gleichungssystem fiir die Gegenflanke dar. 


3. Lésung der Gleichungen. Diese zunachst kompliziert scheinenden Gleichungen vereinfachen 
sich wesentlich, wenn man an Stelle von r,, ,, p, gemaB (10) wiederum die Veranderlichen r,, 
0, Po einfiihrt. (14) geht dann wieder iiber in (13), wahrend die Beziehungen (15) nicht mehr 
notig sind. An Stelle von (16) erhalt man schlieBlich gemaB (10) 


-- =2 + cotg %, cos (vy, + Wyt) = Ecos J, oe (tea) 
1 
: ‘ (17) 
r, sin O,sin(~,+ @.t)—E , E , 1 ao, 
I re aisle aan fe sin 3, 00, tear 


Diese Gleichung 148t sich noch in leicht ersichtlicher Weise vereinfachen, und man erhalt dann 
in Verbindung mit (13) im Tellerradsystem fiir die Gegenflanke die beiden Gleichungen 


Pe= Pp(r,, D2) » (18) 
a a4 + cotg J, cos (P, + Wyt) = Ecos 0, oe (12, Jp) | 
; , (19) 
‘ Jaf 0D 
+ |sin (Py + Wt) me sin 0, 30, (Fonts) | 


Von hier kann man natiirlich gemaB (11) bzw. (12) wiederum zum Ritzelsystem tibergehen. Wir 
schreiben diese Gleichungen nochmals an in der Gestalt 


r, = JE? +72 —2 Er, sin %, sin (p, + 8) , | 


r, sin B, cos (P, + Wot) 


cos oe = 


: (20) 
r, sin 9, sin (Pp, + wt) —E 
tg (M+ at) = ; 


r, cos 0, 
Aus (18) bis (20) 14Bt sich dann im Prinzip die Gegenflanke ausrechnen in der Form 

Py = P(r, 04). (21) 
Die Gleichungen (18) bis (21) stellen nun die Lésung des Verzahnungsproblems bei Hypoid- 
getrieben dar. Zwei Zahnflanken (18) und (21) stehen dann und nur dann miteinander im Linien- 


eingriff, wenn (18) und (21) sich vermége (19) und (20) auseinander herleiten lassen. Die Zeit t 
tritt dabei als Parameter auf in einer gleich noch naher zu besprechenden Weise. 
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Die Bedeutung der Grundgleichungen (18)—(21) liegt in verschiedener Richtung. Die 
wichtigsten Anwendungen sind die folgenden: 

a) Erstens scheint sich zwar im allgemeinen wegen der komplizierten Gestalt von (19) und 
(20) die zu einer gegebenen Tellerradflanke (18) gehérenden Gleichung (21) fiir die Ritzelflanke 
nicht in endlicher Form zu ergeben. Dagegen besteht durchaus die Méglichkeit, da das fir 


spezielle Flanken (18) der Fall ist, soda® man dann in diesen Fallen eine endliche Darstellung — 
fiir die Ritzelflanke in der Form (21) hat. Numerisch ]4Bt sich (21) natiilich stets leicht be- 


stimmen. Man vergleiche dazu Punkt c). 


b) Als wichtigste Anwendung der Grundformeln 1a48t sich in jedem Fall eine Darstellung | 
in geschlossener Form fiir die Ritzelflanke in Gestalt einer Parameterdarstellung angeben. Man 


bemerkt namlich zunichst, da®B gema®B (19) die Zeit t sich aus den iibrigen Veranderlichen 


dieser Gleichung im wesentlichen durch Lisung einer quadratischen Gleichung ergibt. In der | 


Tat kann man ja etwa sin (q+ We t) =u; cos (P.-+ We t) = V1 — u* setzen (vgl. dazu Abschnitt 4, 
S. 61). Zu gegebenen Werten von r, und #, kénnen demnach gemaf (18) und (19) die zu- 
gehirigen Werte von gy, und t in endlicher Form angegeben werden. Aus diesem folgen dann 


weiterhin gemaB (20) die zugehérigen Werte von r,, 3, @,. Diese letzteren GréBen ergeben sich | 


damit also in endlicher Form aus den Parametern #, und gy, womit die genannte Para- 
meterdarstellung der Ritzelflanke gefunden ist. 
c) Die hergeleitete Parameterdarstellung der Ritzelflanke eignet sich in hervorragender 


Weise auch zur numerischen Berechnung dieser letzteren, wenn die Tellerradflanke gegeben | 


ist. Man vergleiche dazu Abschnitt 4 (S. 61) dieser Arbeit. 

d) Bei gegebenen Funktionen (18) und (21) kann sofort verifiziert werden, ob die zugehoérigen 
Flanken miteinander im Eingriff stehen oder nicht. Man geht dazu von (21) aus, ersetzt darin 
r,, Oy, P, durch die Ausdriicke (12) bzw. (20). Weiterhin eliminiert man gemaf (19) t. Die tibrig 
gebliebene Beziehung zwischen ry, 2, P2 mu dann im Eingriffsfall mit (18) iibereinstimmen. 

e) Die exakt geltenden Gleichungen (18)—(21) kénnen approximativ durch einfachere 
Gleichungen ersetzt werden, wodurch sich Naherungsmethoden zur Berechnung von Hypoid- 
getrieben ergeben. 

Die geometrische Bedeutung dieser Entwicklungen ist die folgende: Jedes System der vier 
GréBen r,, Oy, P2, t, das (18) und (19) geniigt, stellt (zunachst im Tellerradsystem) einen Punkt 
der Gegenflanke dar, der dann gemaf (20) auch auf das Ritzelsystem umgerechnet werden 
kann. Gleichzeitig sind r,, 0., ~, aber auch die Koordinaten eines Punktes der urspriinglichen 
Flanke. Demnach ist t die Zeit, zu welcher beide zusammenfallen, also derjenige Zeitpunkt, in 
welchem der durch 15, 0, ~, gegebene Punkt der Tellerradflanke die Flanke des Ritzels berihrt. 
Jeder Punkt der ersteren wird wahrend der Bewegung zum Berithrungspunkt. Die in (19) 
und (20) auftretende Zeit t stellt den Zeitpunkt dar, in welchem dies stattfindet. Daraus ergibt 
sich noch unmittelbar der folgende Schlu8: Wenn man in (19) fiir t irgend einen Zahlwert ein- 


setzt, so stellt die dadurch hervorgehende Gleichung in Verbindung mit (18) die Berithrungs- | 


linie der beiden Zahnflanken zu dieser Zeit dar. 

Den beiden Lésungen der quadratischen Gleichung entsprechend gibt es sogar zwei Lé- 
sungen fiir u bzw. t. Welche der beiden Liésungen in einem konkreten Fall in Betracht kommt, 
erkennt man am besten aus der Anschauung gemafs Abb. 1, unter Beachtung der beiden letzten 


Formeln (10), nach welch letzteren die GréBen u und yl — u? im wesentlichen die Koordi- 
naten 2 und y des gegebenen Punktes im ruhenden System zur Zeit der Berithrung darstellen. 
Demnach miissen in praktischen Fallen im allgemeinen u und in jedem Fall der zugehdrige 


Wert von yl — u* positiv sein. Im folgenden sprechen wir, um die Betrachtungen méglichst ein- 
fach zu gestalten, stets nur von der einen, effektiv in Betracht kommenden Lésung und also 
auch nur von einem Zeitpunkt der Beriihrung. 

Diesen Entwicklungen gema8 erscheint die unter b) angegebene Parameterdarstellung 
in endlicher Form fir die Gegenflanke in einem sehr anschaulichen Lichte. Wie eben 
erwahnt, stellt der Parameter t den Zeitpunkt dar, zu welchem ein beliebig gegebener 
Punkt der Tellerradflanke mit den Koordinaten r,, Jy, g2. zum Berithrungspunkt wird. In 
diesem Moment fallt also der gegebene Punkt der Tellerradflanke mit einem Punkt der Ritzel- 
flanke zusammen, welch letzteren wir als den Gegenpunkt des gegebenen bezeichnen wollen. 
Indem wir schlieBlich diesen gemeinsamen Punkt gema8 den Gleichungen (20) auf das Ritzel- 
system umrechnen, erhalten wir die Koordinaten des Gegenpunktes in diesem letzteren. Die 
genannte Parameterdarstellung lauft also geometrisch gesprochen darauf hinaus, da8 man zu 
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einem gegebenen Punkt der Tellerradflanke zunachst den Zeitpunkt t bestimmt, zu welchem 
er die Gegenflanke im Gegenpunkt beriihrt, und dann gemaB (20) die Koordinaten des Gegen- 
punktes im Ritzelsystem berechnet. 

Erganzend sei noch bemerkt, daB die Grundgleichungen (18) bis (21) nicht nur in den bisher 
verwendeten Polarkoordinaten, sondern auch in Zylinderkoordinaten eine einfache Gestalt an- 
nehmen. Fiihren wir solche Koordinaten in den beiden bewegten Réiumen gema8 Abb. 6 und 
.. i so tritt an Stelle von (18) fiir die gegebene Zahnflanke des Tellerrades eine Gleichung 

er Gestalt 


Po = Pz (02, X2) « (22) 
Durch eine leichte Umrechnung geht ferner (19) itber in 
7) x One. 
iit Geos (Py + 28) = EF (05, x) | 
aw, a ie) 


; Be, 
Tin (Ps + wet) | 0352 (Oe 82) — on Get (Bos) | 


01; Pi: ¥1 = Zylinderkoordinaten 025 Por X, = Zylinderkoordinaten 
Abb. 6. Einfiihrung von Zylinderkoordinaten Abb. 7. Einfiihrung von Zylinderkoordinaten 
im Ritzelsystem. im Tellerradsystem. 


wahrend an Stelle von (20) die Umrechnungsformeln 


1 = 02 608 (Py + Wet) , 


0, = /x2 +[E — o sin (Y, + @,1)F, (24) 

O02 8in (P, + Wet) — E 

te(nt+o) =~ ip s 
treten. SchlieBlich erhalt man fir die Ritzelflanke entsprechend (21) eine Beziehung der Form 
py = Fy (0 1) + (25) 


Die Gleichungen (22) bis (25) sind das vollstandige Analogon zu (18) bis (21) und kénnen in genau 
derselben Weise verwendet werden. 


4, Bemerkungen zur numerischen Berechnung der Gegenflanke. Wie bereits erwahnt, sind 
die Grundgleichungen (18) bis (21) hervorragend geeignet fiir die Lisung der wichtigen Aufgabe, 
zu einer gegebenen Zahnflanke die Gegenflanke punktweise numerisch zu berechnen. Aus 
Griinden der Systematik wurde bei unseren Betrachtungen stets die Tellerradflanke als die ge- 
gebene, die Ritzelflanke als die gesuchte angenommen, aber selbstverstandlich ist das Verhaltnis 
ein absolut gegenseitiges. Die umgekehrten Formeln ergeben sich aus den oben abgeleiteten in 
leicht ersichtlicher Weise. 

Um aus einem Punkt P (r¥, 0%, py) der Tellerradflanke mit py = ©, (r3, 02) gemaB den Be- 
trachtungen des letzten Abschnitts den zugehérigen Gegenpunkt Q (r7, Jf, pi) der Ritzelflanke 
und die Zeit «* der Beriithrung zu berechnen, schreibt man Gleichung (19) am besten in der Form 

cos (P2 + Mgt) = sin (Pp, + Wt) tg Dy a (2, V2) | 
j (26) 


ap Hee vod ta Oe | 
+tg 0, a + E cos Oo Fr (a6 ah es sin ?, 38, (2 04)) | 
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Fiihrt man die Abkiirzungen ein 
u* = sin (gt + 0,0), 
0, 
A = tg 03 (r3, 99), (27) 


a0, 
= yz | Pe 2 OP.) on eels gx OP.) o* 
B= te 05 apatmal St Pherae) 2) 7 sin 0% 00, (3, 3) 9 


und quadriert (26), so ergibt sich nach leichter Rechnung zur Bestimmung von u* die Gleichung _ 


(42741) u¥®+2ABu* + BH 1=0. (28) | 

Daraus erhalt man u* in der Gestalt 
A AB A? — B24 1 
a Ltn Gt +o) 2S (29) 


woraus sich weiterhin cos (v3 + w,t*) und t* ergibt. Das Vorzeichen in (29) sowie dasjenige von 


cos (ps + w,t*) ergeben sich wie erwahnt leicht aus der Anschauung. Im allgemeinen miissen 
sin (y3-+ w, t*) und cos (y3-+ w, t*) positiv sein. ; 

Damit ist die Zeit t* berechnet, zu welcher P(r3, 03, p>) zum Beriithrungspunkt wird. Aus 
dieser und den Koordinaten von P berechnet sich dann der Gegenpunkt Q(ri, 07, yj) sofort 
gemadf (20). 

Fir die praktische Durchfiihrung der Berechnung von Q aus P halt man sich am besten an 
das folgende Rechenschema: 


: ; o® | 
1. Festlegung der Flanke g,=@,(r,, 0.) und der partiellen Ableitungen a (rz, 8) und | 


a®, 
00, (Tp, Do). 


2. Wahl eines Punktes P(r}, 03, y>) der Flanke mit y} = @, (r3, 0%). 
3. Berechnung von 


aD 
ce * IP 2 9x 
A= tg Ces 00, (r3, 2) ’ 
ce) o®D E 0®, 
B= tg 08/T* 4 Ecos OF 3p: (78, 98) — za sin OF age PE 88) 
he ys pe ey 
ee eed le ay: 528 
daraus sin(g} + wet*) = — C + jC? D und die Zeit t* der Berihrung. 


4. SchlieBlich Bestimmung des zu P gehérenden Gegenpunktes Q(ri, J}, yt) gemaB 
rt = jE + 73 —2 Er}? sin JF sin (yy + wyt*) , 


rx sin OF cos (py + @,t*) 


oa 
cos UF * > 
1 
P r% sin O* sin (p* + w, t* 
tg (v7 + @,t*) =- a ee ; , 
r* cos O* 


Auf diese Weise kénnen beliebig viel Punkte der Ritzelflanke auf einfachste Weise (speziell 
auch durch Verwendung von Rechenautomaten) berechnet werden. Ein ganz analoges Schema 
ergibt sich auch bei der in Abschnitt 3 erwahnten Verwendung von Zylinderkoordinaten. 


5. Zusammenfassung. In der vorliegenden Arbeit wird das Problem behandelt und ge- 
lést, zu einer gegebenen Zahnflanke eines Hypoidgetriebes die zugehérige Gegenflanke rech- 
nerisch exakt zu bestimmen. Die Lésung, die auf Grund des Enveloppen-Prinzips hergeleitet 
wird, wird dargestellt durch die Grundgleichungen (18)—(21). Diese kénnen in verschiedener 
Weise verwendet werden. Als wichtigste Anwendung ergibt sich aus ihnen die unter b) auf S. 60 
behandelte Parameterdarstellung in geschlossener Form fiir die Gegenflanke. Diese ist auch 
fiir die numerische Berechnung auferordentlich geeignet. Die geometrische Bedeutung dieser 
Parameterdarstellung ist die, da®B man zu einem beliebigen Punkt der gegebenen Flanke zu- 
nachst als Parameter die Zeit berechnet, zu welcher er die Gegenflanke im Gegenpunkt be- 
rithrt (im wesentlichen durch Auflésung einer quadratischen Gleichung), und daB man dann 
auf Grund der Kenntnis dieser Zeit nach den bekannten Transformationsformeln den Gegen- 
punkt selbst berechnet. Man erhalt auf diese Weise schluBendlich also den Gegenpunkt in Ab- 
hangigkeit vom gegebenen Punkt. 


(Eingegangen am 7. Juli 1952.) 
Anschrift des Verfassers: Professor Dr. P. Matthieu, Berlin-Charlottenburg, Jebensstr. 1. 
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Die Momentenflaiche der elastischen Platte oder Pilzdecke 
und die Bestimmung der Durchbiegung aus den Momenten. 


Von W. Miiller. 


1, Einleitung. In einer friiheren Mitteilung! ist die Durchbiegung fiir mehrere Belastungs- 
falle von rechteckigen Platten und durchgehenden Pilzdecken dargestellt worden. In Ergangung 
dazu wollen wir vor allem zwei bemerkenswerte Sonderfille herausgreifen, den Fall der frei 
aufliegenden Platte mit Punktbelastung (Einlastplatte) und die durchgehende Pilzdecke mit ver- 
schwindenden Stiitzflachen, um namentlich die Momentensumme oder die Niveaulinien der 
entsprechenden Flachen mit Hilfe der Thetafunktionen® bis zur zahlenmaBigen Auswertung zu 
ermitteln und schlieBlich an dem Sonderfall der vereinfachten Pilzdecke die Frage zu behandeln, 
wie man die Umkehrung der 4-Operation (4—!) in einem vorliegenden Fall auszufiithren, d. h. 
die Durchbiegung aus den Momenten unter Beachtung der Randbedingungen berechnen kann. 


2. Die Momentendarstellung fiir die Einlastplatte. Wir haben fiir die im Punkte x) =a é), 

Yo = NH durch eine Einzelkraft P belastete, frei aufliegende rechteckige Platte mit den Seiten a 
und b die Durchbiegung 
3 oo . . 

= — > ens : {Gin n2x(A — mp) (1+ naa Ctgnazd) SGinnay 


ad —nxzn8oj nan] —nx(A— mp) Goi naz (A— 7) Sitnan} (1) 


Ww 


gefunden, wo b/a =A das Seitenverhaltnis des Rechtecks bedeutet. Der auf der rechten Seite 
stehende Ausdruck ist eine biharmonische Funktion, die in allen Stellen auBerhalb des Last- 
punktes der Gleichung 


AdAw =0 (2) 
geniigt. Dabei ist, wie die Ableitung ergab, 7 < 7 vorausgesetzt. Fiir 7 > m erhalt man den 
Ausdruck fiir w durch Vertauschung von 7 und 7p. 

Fir die Momentensumme M ergibt sich dann 


mx my _2P Qsinnz&sinnaé Ginnan Sinn x(a — mh) 
M— I1+y oo n Ginna : (3) 


1 


Wenn man die komplexen Gréfen € =& +17, Co =& +1 und die entsprechenden kon- 


jugiert komplexen Koordinaten € =E— LN > C; = &)—in, einfihrt und den Realteil einer 
komplexen Gide mit Jie bezeichnet, so kann man fiir den Zahler unter dem Summenzeichen 
folgenden Ausdruck setzen: 

+ Re [eos na(¢ —Ly-+ id) + cos na(l +£,— it) —cosna(l—f,— if) —cosna(E+Cy+ iA). 
In der Theorie der elliptischen Funktionen? ist es iiblich, das Periodenverhaltnis t und die 
GréBe gq =e”? einzufithren. In unserem Falle des Periodenrechtecks wird 7 = tA =i b/a, also 
gq =e-7*, Dann erhalten wir 


M== >)" __ Me [eos na(€ —C,+ia)+cosna (€+C)—tA) 
mesg A Oe i) —cosna(l —C,—idA) —cosna(6+0)+iA)]. (4) 


Benutzt man nun die bekannte Reihenentwicklung fiir den Logarithmus der ?)-Funktion 


S 1 qn cos 2 
In 8,(w) = _>"In(1 — q?) —2 gq” cos 27n w 6) 


- 7 mb <= Gan) 


1 W. Miiller, Ing.-Arch. 20, (1952) 5. 278. eae 
2 Vgl. etwa G. H. Halphen, Traité des fonctions elliptiques et de leurs applications, Bd. I, Paris 
1886—91; ferner Fr. Tricomi, Elliptische Funktionen S. 156f, Leipzig 1948; Fr. Oberhettinger— W. Magnus, 
Anwendungen der elliptischen Funktionen auf Physik und Technik S. 13 und 35—37, Berlin 1948. 
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und weiter die Formel 
Do (~ ee 27) = fpig-Met'*" dy (w), (6) 
indem man fiir w die Argumente 4(¢ —,), 4(C-+40), $(€ —20)> $@ +00) und r= iA setat, 


so sieht man, daB nach Weghebung der in (6) auftretenden Faktoren von 0, die Momenten- 
funktion in der iibersichtlichen Form sich darstellt 


P| eet) 
w= 5 Re In 3 ez . (7) 
(=) a( 


fey) : 7) 


50 10, 


30r =f 
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60 
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Q7 
60 S 1—0 


50 7 


7 
45, 
Te OR UR OF 12 0=CO# 16 G0 G0 Ge G6 68 $0 32 gf 


Abb. 1. Darstellung der Funktion f (£7) = 21Im| 0, (¢/2)| + ky. 
a) fiir kleine Werte von 7 (< 2). b) fiir gréBere Werte von 7 (> 2). 


Da diese Funktion nur von den absoluten Werten der ,-Funktionen abhangt, so weist man 
leicht nach, daB die Unterscheidung der Gebiete 4 S no tiberfliissig wird, oder daB die Formel 
im ganzen Gebiet des Rechtecks Giiltigkeit hat. Die z. B. aus der Membrantheorie bekannte 
Formel (7) findet sich auch in der sehr umfassenden Arbeit von Fr. Télke, in der wohl zum 
ersten Mal die 9-Funktionen auf viele Einzelfalle aus der Plattentheorie angewendet werden}. 
Statt der ?-Funktionen kann man iibrigens auch die o-Funktionen mit veranderten Argumenten 
und schlieBlich auch die WeierstraBschen p-Funktionen einfiihren. Es wird dann 
Ls o(z — 2) o(z + %) Pp P (2) — p (2) 

M <3 8 Ga) ole al ae PGP) © (8) 
Diese Darstellung findet sich in einer verwandten hydrodynamischen Untersuchung von G Jaffe 
und W. Miiller?. 


1 F, Télke, Ing.-Arch. 5 (1934), S. 187. 


2 G. Jaffé, Ann. Physik (4), 61, (1920), 5. 173; W. Miiller, Z. angew. Math. Mech. 8 (1928), S. 447 
und 10 (1930), S. 227. 
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aa praktischen Berechnung der M-Fliche geht man am zweckmifigsten auf die 0,-Reihe 
zuriick : 


2n+1\? 


= n Bis) J 5 
0, (w) = 2 > (—1)"q sin (2 n+ 1)mw=2 g'/4(sin ww —g@sin 3w + gsin Sw —-+-.). (9) 
0 


Wenn man etwa 4 — 1,5 setzt, so wird g = 0,00898, und man sieht, daB man sich wegen der 


ausgezeichneten Konvergenz auf die beiden ersten Glieder der Reihe beschranken kann. Man 
hat dann mit w =u +iv 


0, (w) = 2 q'!4[sin x u Coj av +icosa u Sin x v— g?(sin 37 uCoj 3av+icos3au Cin3 x v)}. 
(10) 


Abb. 2. Niveaulinien der Funktion 21n| @, (¢/2)|in der Umgebung des singularen Punktes & = 0, 7 = 0. 


Setzt man 0, =, + iy,, so wird 
Q=2 ¢!4 (sin x u Gof xv — g? sin 37 wu Coj 3.2 v), 
y,=2q4(cosmu Cinav—qcos3nu Gin32). 


Daraus folgt nach entsprechender Umformung und Reduktion auf die Argumente mit n= 2 
|P, (w) |? = gy? + yp? = 4 q!? (Co 2 1 v — cos 2 wu) Games — 2 qCoj 2avcos2x u) foil (19 
Mit w —€/2 findet man also 


d, (zy ~ 4 g!!2 (Coj a4 — cos x €) bss —g—2¢qCoj x4 cosa é) ‘ (12) 


Wir haben zunachst die Funktion 


els) ual (2) 


GD qi4 


21n + ky = f (€n) 


tir eine Reihe von £-Werten in Abhangigkeit von 7 berechnet (Abb.1) und daraus die Niveau- 
inien der Funktionen 2 lg 3,(¢/2) veranschaulicht in der Umgebung des singularen Punktes €—0, 
5 
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in dem #, verschwindet, also In 0,(¢/2) negativ unendlich wird (Abb. 2). Im ganzen ist das 
Feld periodisch mit der Periode 2 a und gleichwertig mit der Potentialstrémung einer unend- 
lichen Reihe von aquidistanten Wirbeln in den Punkten x= 0,2a, 4a usw. 

Eine zweite Wirbelreihe liegt in den 
Punkten € =2n-+2iA. Aber in der 


zitat. Aus der Formel 
8, (w+ 2) =— qr} 4 8, (w) 


ergibt sich 
Ke In i(y — ia) == Inge | 


+ + Re In J, & | ‘/ 


Es tritt also noch auBer der Konstanten 
ein in 7 lineares Glied hinzu, das hydro- 
dynamisch einer Parallelstrémung ent- 
spricht. Das Feld in der Umgebung des 
singularen Punktes ¢=21/ haben wir 
in der Abb. 3 dargestellt. 


%—%N =b, n+ =b,, so werden 
die Quadrate der Absolutbetrage der 
vier in (7) auftretenden #,-Funktionen 


Abb. 3. Niveaulinien der Funktion 2 In|), (¢/2)| in der Umgebung des 
singularen Punktes § = 0,7 = 2d. 


“a 4 gil (Gof xd, —cos ma) — g —2 g? Co aby cosa). 


—¢ 
ba oe) 2 


0, (4%) 


5 Sh q'/2(Coj a b, — cos 7 ag) 5 — g—2 ¢Woj 7b, cos 7a) ; 


eee ND: 
,(F5% =~ 4 q'l? (Coj a b, — cos 7 ay) is — g— 2 q? oj a b, cos x a ‘ 


| > \\2 
6 (=e) 


| 
| 
a) 
| 
| 
| 


~ 4 q'!2(Coj a b, — cos 7 ay) (Ea g—2 q? Gof a b, cos 7 ay). 


Da $(¢ —,) einen Wirbelpunkt in ¢, gibt, so kann man sich das der Formel (7) entsprechende 
System der Niveaulinien dadurch entstanden denken, daB man dem im Lastpunkt des Grund- 


rechtecks gelegenen Wirbel die drei durch einmalige und zweimalige Spiegelung an den Seiten — 


— =0, 7 =0 entstandenen Wirbel hinzufiigt, wobei bei einmaliger Spiegelung der Drehsinn 
des Wirbels sich umkehrt (Abb. 4). 

Die Abb. 5 bezieht sich ferner auf die Héhenschichtlinien der Momentenflache fiir den 
Fall A =1,5,, &) =0,4, my) =0,5. Der Verlauf dieser den singularen Punkt € =C, umgebenden 
Linien, zu denen iibrigens bei freier Auflagerung auch die Grenzen des Rechtecks gehéren, wo 
M verschwindet, hat eine Ahnlichkeit mit den Héhenschichtlinien der w-Flache, unterscheiden 
sich aber von diesen Linien namentlich in der Umgebung des Lastpunktes, da fiir w hier der 
Grenzwert ! 


P 
W =p prior (15) 
zu nehmen ist, wenn unter r eines Nachbarpunktes vom Angriffspunkt, also | r | == tS —Cy| 
zu verstehen ist, wahernd M in der Nahe von € =0 durch 
IB 
M, = oi In Tr 


dargestellt werden kann. Man sieht, da8® w fir r— 0 verschwindet, wahrend M unendlich 
groB wird. Der Lastpunkt ist also fiir w keine eigentliche Singularitét; er hat aber fiir das 


Sei ferner 6 —£,=0, € +& =a,, | 


n-Richtung herrscht keine reine Periodi- 


* Vel. A. Pucher, Ing.-Arch, 12 (1941), S. 76f; ferner A. Pucher, EinfluBfelder elastischer Platten, , 


Wien 1951, 
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Moment den Charakter eines Wirbelpunktes, dem die M-Linien mit ins Unendliche wachsender 
Dichte sich asymptotisch annahern. 


3. Die Momentenflache fiir eine durchlaufende Pilzdecke mit Punktstiitzen. Wenn man den 
_Ursprung der Koordinaten in den Mittelpunkt des Periodenrechtecks 2a-2b legt, deren vier 
Eckpunkte mit den Punktstiitzen zusammenfallen und wenn man die fritheren Bezeichnungen 
einfithrt, so erhalt man die Durchbiegung fiir diesen Fall in Form 

be p b4 n\2|}2 pab (—1)"cosna& 
w =“oty en} = | Nat Oy a8 Girma Wl + 0d Cig nad) Cof nxn 
; —naynSinnan], (16) 

deren erster Fall Teil als Partikularlésung der inhomogenen Gleichung 


AAw = 


sich auf einen Plattenstreifen parallel der x-Achse von der Breite 2b bezieht, wahrend der zweite 
Teil als biharmonische Funktion der homogenen Gleichung geniigt. 


made fe Abb. 5. Hdhenlinien der Momentenflacbe 
Abb. 4. Zum Spiegelungsprinzip bei der Konstruktion der Einlastplatte. 
der Momentenflache einer Einlastplatte. (A ARS 0,4, "0 0,5) 


Daraus folgt dann fiir das Moment 


oc 


1 1 2pab (—1)"cosnaEGojnayH _ = 
M=pa(hi st) — 7 = n Ginna ee eware Be 
=, 


Mit Einfithrung von q haben wir ferner 


1 1 4pabo,. = (—1)" q" cosnae 18 
M=pa(tp 5 n?) a (lege) ; ( ) 


n= 


Die rechts stehende Summe kann man durch den Logarithmus der Funktion #, (¢/2) ausdriicken 
entsprechend der Formel 


a _ hes 1)"q" cos 2naw 19 
TD aC ae ‘eget aes (19) 

Dann wird ] 
M= parla, Z—2 Re ln B, (y)|=PebFn. (20) 


68 Miller: Die Momentenfliche der elastischen Platte oder Pilzdecke usw. Ingenieur-Archiv 


Setzen wir 0, (C/2) = +i ws, so kénnen wir schreiben 


jy 2 
F,, =In (93 + 93) +4 — 55 - 


Benutzen wir ferner die 3,-Reihe 


9,(£) =142 >'q"cosnal=1+42 qceosa6 +2 qtcos2aC4+2q@cos3aC, (21) 
1 


Bi 
aos s) so wird bei Vernachlassigung der 
20 92 héheren Potenzen von q 
3% 1+2qceosxéESojan, 
cs | y,~—2qsinné Sinan, 
16 also 
ga tye~1+4qceosaéE So) ay 
oe + 4 q? (cos? x € Co}? a 7 
of + sin? 7 & Gin? xn) (22) 
=1+4qcosxé oj ay 
10 | + 2 q? (oj 2an-+ cos2z76€). 
8 In dem angenommenen Fall J = 1,5 
kommt dann z. B. 
96 Fg gp + y2~ 1+ 0,0359 cosa ECoj xy 
7 y .0,00016 (Gof 2277 + cos 27). (23) 
ie Man kann auch die logarithmische 
pe 45 Reihe direkt benutzen und erhalt 
y damit bei Beriicksichtigung der vier 
: 06 ersten Potenzen von q 
G7 t 3 
-G2 GS Fo ee ee 
P % In 0, ( 5 ) q ul 
~o4 +2 q(1+ q?) cosa (24) 
oat —gPcos2nt +5 gq’ cos ae 
-08 woraus dann 
é=7\| \gs t eee 
= “eed 
-10F Me Ind, (5) = tacos 
2q(1 e 
Le + al + q’) cosa éCoj a7 (25) 
—@cos2rECoj 2a : 
-14 et - 4 <P cos 3x € Gof 307 
E sich ergibt. 
as. Wir haben zunachst, in Abb. 6 
CO OR GE Gi BE 0 ne Ee den Verlauf der Funktion 2 In ,(¢/2) 
Abb. 6. Darstellung der Funktion 2 In |, (¢/2)| in Abhingigkeit von § und y, £4 verschiedene Werte von ¢ in 


‘ht Abhangigkeit von 7 und daraus die 
Niveaulinien in der Umgebung des singulaéren Punktes¢ =1+i/ aufgezeichnet (Abb. 7), in 
dem entsprechend der Beziehung 


it Tt 1 sdk 
der Logarithmus negativ unendlich wird. Im Ganzen erhalten wird das Strémungsbild einer 
Reihe von gleichsinnigen Wirbelpaaren (€ = + (1 +2n-+ iA)), die symmetrische zur é- und 


n-Achse liegen. In Abb. 7ist nur der Teil 0 <<& <2,7> 0 dargestellt. Das vollstandige Feld 
ergibt sich durch periodische Vervielfachung in der positiven und negativen ¢-Richtung und 
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durch Spiegelung an der 7-Achse. Fiir gréBere Werte von n gilt die Beziehung 
é ; 4 
Re In (yp +12) =Ing'+27+ Re ln (5) ‘ (27) 
Ubrigens geben die Funktionen Do (C/2) und #, (¢/2) die gleichen Bilder wie die benutzten Funk- 


tionen, nur mit verdnderter Koordinatenlage, entsprechend den Umrechnungsformeln 


I (x + 3 = 9,(w), #,(w a e = 0,(w) (28) 


Aus den Werten von In #,(¢/2) sind dann die 
| 


entsprechenden Werte der Momentenfunktion F,, 
(vgl. Tabelle) berechnet und die Héhenschichtlinien 
des Momentes fiir den vierten Teil des Perioden- 


ee: rer 


4 \ 
4 a \ 
79033 G033N. Ms \ | 
7h Sali as7 ae | 
4 a) MR 2 \ | 
N Cf ‘N 
Na ‘ 5 
=7 ee 7 0 a 
Abb. 7. Niveaulinien der Funktion 2 In |, (¢/2)| in der Umgebung Abb. 8. Héhenschichtlinien der Momentenflache einer 
der singularen Stelle € = 1, 7 = 1,5. durchlaufenden Pilzdecke (A= 1,5). 


-echtecks aufgezeichnet (Abb. 8), Die Vervollstandigung des Bildes ist durch Spiegelung an 
Jen Rechteckseiten leicht herzustellen. Die Kurvenschaar unterscheidet sich wieder von der 
sntsprechenden, zur Durchbiegung w gehérenden und friiher dargestellten Schaar durch den 
singuldren Charakter des Punktes § =1,7 =A, dem eine ins Unendliche gehende Spitze der 
Momentenflache zugehért. Erst bei Beriicksichtigung der endlichen Stiitzflache wird das Mo- 
ment iiberall einen endlichen Wert annehmen. 

Fiir die Mitte des Rechtecks (£ =7 =0) wird F,, = 0,261. Die konstante Differenz fir 
fF, zwischen zwei benachbarten F,,-Linien betragt — 0,2/~ = — 0,0637. Das Vorzeichen 
wechselt ungefahr langs der Linie 5. 


Tabelle der Momentenfunktion Fm bei der Pilzdecke. 


n €=0,1 E=0,2 E€=0,3 E=0,4 E=050 E=0,6 E=0,7 E&=0,8 co 

0 0,261 0,2608 0,2567 0,2536 0,250 0,2465 0,2436 0,2407 0,2380 
0,2 0,2503 0,2477 0,2447 0,2410 0,2366 0,2325 0,2285 0,2254 0.2228 
0,4 0,2180 0,2140 0,2094 0,2036 0,1969 60,1900} 0,1836 0,1791 0,1746 
0,6 0.1675 0,1605 0,1528 0,1394 0,1310 0,1224 0,1076 0,0985 0,0901 
0,8 0,1019 0,0882 0,0774 0,0592 0,0382 0,0161 | —0,0057 | —0,0245 | —0,0406 
1,0 0,0366 0,0188 | —0,0032 | —0,0328 | —0,0702 | —0,1078 | —0,1550 | —0,1933 | —0,2315 
1,2 | —0,0212 | —0,0468 | —0,0791 | —0,1290 | —0,1857 | —0,2600 | —0,3471 | —0,4372 | —0,5450 
1,4 | —0,0551 | —0,0866 | —0,1268 | —0,1866 | —0,2690 | —0,3800 | —0,5330 | —0,7422 | —1,2740 
1,5 |—0,0598 | —0,0920 | —0,1337 | —0,1950 | —0,2793 | —0,3990 | —0,5631 | —0,8080 | —co 
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4, Die Ableitung der Biegefliiche aus der Momentenfliche. Im Allgemeinen wird es méglich 
sein, z. B. bei einer rechteckigen Platte mit Punktbelastung die Funktion M als Potential- 
funktion oder als Realteil einer komplexen Funktion W = M +i N mit Hilfe der bekannten 
Singularitiitenmethode und des Spiegelungsprinzips aufzustellen. Nicht so einfach gestaltet 
sich die Ableitung der im Allgemeinen biharmonischen Funktion w aus der Momentensumme. 
Wir wollen hier ein Verfahren angeben und auf den Fall der Pilzdecke anwenden, mit dem es 
durch eine doppelte Integration midglich sein wird, die Biegefliche aufzufinden. Dabei kénnen 
wir uns auf das zweite Glied M,in dem Ausdruck fiir M beschranken, weil der sich auf den Platten- 
streifen sich beziehende Teil M, ohne Schwierigkeit w, zu berechnen gestattet. 

Wir gehen davon aus, die beiden konjugiert komplexen GroéBen € = & +i 7 und —F =f i ” 
als unabhangige Variable einzufiihren. Damit kénnen wir zunachst M, die Form geben 


“Presb (—1)"(cosnaé + cosnzZ) 29 
Ya Tt n Ginna : CZ 


Um die Operation Aw, in den Variablen darzustellen, hat man 


du, Ou, , Ou, Qw, __ ; (Ow, du, 
a gee ar a (a z) 
daher 
Ow, A 0? w, 2 Pw, “2 Pw, Ow, we L. Pw, Pw, 
Gee eee Oe ae nce Ont” 7 PEE hs Weed h ool 
Also 1a48t sich die 4-Operation in folgender Weise darstellen: 
4 Pw, 
w, =—=——. 30 
Bae ae: (30) 
Daher wird 
4 N 0?w 
Wh == = IN Alie =—=— Sn 31 
2 ‘ 2 a2 ar ar ( ) 
und man sieht, daf es theoretisch méglich ist, umgekehrt w, als Doppelintegral 
2 —_ a 
=—in {| M,(¢&) dé dé (32) 
zu berechnen. Die Lisung ergibt sich in der allgemeinen Form? 
w, =CF (2) +CF (0) + AC) + flO) - (33) 


Die Schwierigkeit besteht nur darin, die auftretenden willkiirlichen Funktionen den Rand- 
bedingungen anzupassen. Wenn wir (26) in (29) einsetzen, erhalten wir zundchst fiir w, folgenden | 
Ausdruck: 
a a b (— 1)” ye : i i } 
2 IN 2u an) Gintnna te sinnal + nacsinnal + fi(C) +f2(6)] 
; (34) 
J 


a® b (2 1)FECe 
4N (an)? CinnawA * 


Die Funktionen f, und f, miissen so bestimmt werden, daB 0G/0& und dG/dn fir & —0 und 
1 und fiir 7 =0 und 4 verschwinden. Wir haben aber 


0G : 
ag = 2 nasinnz &Coj natn —2(nz)?ysinnzwé Cinnazy | 


+ 2(nz)? €cos nz & Co} nan +f; (C) SEE: P 
: 3 
I nacosnné Sin nay +2(nz)?é sinnaé Sin nay | “ 


+ 2 (nx)? y cos nx Col nxn + iffi(S) —f(6)) - 


Fir den Mittelpunkt des Periodenrechtecks wird also 


Filin) + fx(— in) = 0 ; Fi(S) —felo) = 9. 


} Vgl. auch W. Miiller, Einfiihrung in die Theorie der zahen Flissigkeiten, S$, 121f., 215, Leipzig 1932, 
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Daraus folgt aber, daB man f, =f,—=f setzen kann, und daB f eine gerade und f’ eine un- 
gerade Funktion annelle: Setzt man f= +iy, so folgt aus 
0G 
(FF), =2(nz)®cosnaCojnxn +f'(A+in) +f/(A—in) <0 
die Beziehung 
y (1, 7) =— (na) cos naz oj nan, (36) 
ebenso aus 


ac 
ale nzcosnz& SinnwdA(1+naAaCtg nad) +2(nx)2EsinnazéSin nwd 2 y'(EA) 


die Beziehung 
y (a) =na(1+naaCtg nad) cosnné GinnazA+ (na) Esinnazé Ginna. (37) 
Setzt man nun 
f(C) =—Deosna6—Enaésinnal, 
so findet man 
y (A, n)=— E(na)cosnaCojnan, 
yp (€,4)=(D— E) nacosnxé GinnazdA+ E(nz)[Esin na é Gin nad —Acosnaz ECoj nz A]. 
Durch Vergleich mit den Werten (33) und (34) wad dann 
B=), D—E— EnxdAGtgnzA =1+4 nxiCignad, | 


also (38) 
DS wai Cianel: | 
Es kommt mithin schlieSlich 
f(6) = —2(1+ na AGig nad)cosnal—nalsinnazl. (39) 


Damit erhalten wir nach (34) in*Ubereinstimmung mit der Ausgangsformel 


2 = 3b Ware howe Ea tc : F) 
m= | | M, dé ae 2 ea — [nax(€ —C) (sinnal —sinnzC) 


+2(1+nz2AStq nz) (cosnal + cos nxC)| (40) 


pab See) econ ey 
IN — (xn)? Ginn awd 


1+n2AStqnxd) Coj nay —nay Ginna]. 


Man kann nun die Endformel auch durch die einfachen und doppelten Integrale der ?,-Loga- 
rithmen darstellen. Geht man von der Formel aus 


(—1)"* 1G" cosn al IG 
2 n(l — 2”) == [hin ,(¢ +G,; (41) 
Bo Co Sn (1 — q?") , so erhalt man durch Integration zwischen den Grenzen 0 und ¢€ 
1 
~1 ; 
(—1)P"q"sinnaz g | EG 
cE == Tie See [mo($) +2 (42) 


0 
und durch eine weitere Integration 


n+1lon 


o0) =2, St = [fmol Jaz dt — 4. C,t?+€ 


(43) 


weil 


Set ee I 
C=2 >) n?(1 — q?”) SOS cians 


1? Miiller: Die Momentenflache der elastischen Platte oder Pilzdecke usw. Ingenieur-Archiy 


Beachtet man ferner, daB der von A abhiangige Faktor im zweiten Ausdruck von w, in der Form 


linaaGignna 1 42 1 Ky pea ee i 
Ginna ~~ Cinna dA\ Ginnazda} — 0A} Ginn awd 
a 1 
=(1+%495)) canae (i 


geschrieben werden kann, so erhalt man mit den neu eingefiihrten Funktionen J/, und @, fol- 
gende Darstellung fiir die Durchbiegung: 


a®b 3g e : 
1 = EET (C—O) UEC) —Ha@)] +2 (1-257) 104(0) + 0.) 7 


oder in der Form eines Realteils einer komplexen Funktion 


Ww, = ae Me EE) NEDA) (0) 
2 oa Nt 3 ( a1) ae 
Die an anderer Stelle durchzufiihrende Berechnung und Tabulierung der Integralfunktionen //, 
und @,, die meines Wissens in den bekannten Handbiichern der elliptischen Funktionen nicht 
vorkommen, macht bei der guten Konvergenz der zugehérigen Reihen keine groBe Schwierig- 
keiten. Entsprechende Funktionen lassen sich in Verbindung mit den anderen #-Funktionen 
aufstellen. 

Die Erweiterung dieser Betrachtung und die Behandlung anderer Belastungsfalle recht- 
eckiger Platten oder Decken nach den vorgetragenen Grundsatzen, insbesondere also die Zuriick- 
fiihrung der Durchbiegung auf endliche Ausdriicke mit Hilfe von charakteristischen Funktionen 
miBte den Gegenstand einer besonderen Untersuchung bilden. 


(46) 


(Eingegangen am 10. Juli 1952.) 
Anschrift des Verfassers: Prof. Dr. W. Miiller, Miinchen 9, Taubenstr. 41. 
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Fachwerk und Rahmenwerk. Ein syste miatiseher Grundri8 der Statik des ebenen ~ 
Tragwerkes. Von Dr.-Ing. Walter Fries., Mit 365 Bildern, X, 368 Seiten. 1953, 
Ganzleinen DM 42,— 


Inhaltstibersicht: Einleitung. Aufgab i i 
L rsi : 1 : gabe und Umfang der Statik. 1. Grundlagen: Das Tragwerk und seine 
Teile. Standfestigkeit. Kraftwirkungen. Die Waunueletoland: Hauptsiitze der Mechanik, Teretecck und Seileck. 
a) poe rragnectadas und Zerlegung von Kriften. b) Statisches Moment und Seileck. Spannungen und Form- 
a erungen des Stabes, a) Der gerade Stab, b) Der krumme Stab. Die Hauptgleichungen‘der Statik des ebenen 
hi ge Die Aufgabe der Statik. Lésungsverfahren, — |, Das einfach standteste Tragwerk: Allgemeines, 
gomcine Berechnung mit der Hauptgleichung. Allgemeine Berechnung mit den Gleichgewichtsbedingungen 
er Knoten. Die Stitzkrifte. Die inneren Krafte des Tragwerks. Stabvertauschungsverfahren. Der einfache 
‘Savant EinfluBlinien, Der einfache Trager unter einem verschiebbaren System von Einzellasten. Das Dreieck- 
‘ac werk. Mehrfache und mehrteilige Fachwerke. — II, Kinematische Untersuchung des einfach standfesten Trag- 
werkes: i. Zwangsliufige Ketten. Die zwangsliufige Kette. Der Polplan, Der Geschwindigkeits- und Ver- 
schiebungsplan. 2. Anwendungen. Berechnung der Stabkrifte aus dem Geschwindigkeitsplan. Kinematische 
Ermittlung der EinfluBlinien. Das kinematische Kennzeichen der Standfestigkeit. — JV. Die Forminderung des 
Tragwerkes: 1. Allgemeine Berechnung der Formanderung. 2. Besondere Verfahren. Fachwerketabziige. Die 
Biegungslinie als Momentenlinie (Seileck). Die Forminderung des geraden Stabes. Rahmenstabziige. — V. Das 
mehrfach standfeste Tragwerk. Kraftverfahren: Das Kraftverfahren. Lastengruppen. Das Festpunktverfahren. 
a) Unverschiebliche Knotenpunkte. b) Verschiebliche Knotenpunkte. Beispiele., Der eingespannte Stab! Der 
eigespannte Bogen. — VI. Das mehrfach standfeste Tragwerk, Formanderungsverfahren: Die Grundwerte und 
die Hauptsitze. Die Hauptgleichung. Stockwerkrahmen. Formanderungsgruppen. Naherungeverfahren. a) Un- 
yerschiebliche Knotenpunkte. a) Das Momentenausgleichsverfabren. $) Das Drehwinkelausgleichsverfahren, 
b) Verschiebliche Knotenpunkte. Beispiele, Nachtrag. 
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Grundlagen der Wasserwirtschaft und Gewdsserkunde. von Professor 
Dr.-Ing. Otto Streck. Mit 291 Abbildungen. X, 466 Seiten, 1953, Ganzl. DM 43,50 


é 4 


Inhaltsibersicht: #. Wesen und Zweck der Wasserwirtschaft: Das Wasserdargebot aus dem Wasser- 
kreislauf. Die Nutzungen des Wassers. Schadenverhiitung bei Wasser. Wesen und Aufgabe der Wasserwirtschaft. 
Praktische Beispiele fiir die Wasserbewirtschaftung. Die gewdsserkundlichen Aufgabengebiete. — II. Die atmo- 
spharischen Elemente des Wetters und das Wettergeschehen: Die atmosphirischen Elemente’ des Wetters. Das 


. Wetter, — III, Der Wasserhaushalt und seine Elemente: Niederschlag. Verdunstung., Versickerung. AbfluB. Be- 


ziehung zwischen Niederschlag, Abflu8 und Verdunstung. — IV. Allgemeine Eigenschaften oberirdisch flieBender 
Gewisser: Talbildung und FluBlaufgestaltung. Laufentwicklung. Bett, Lings- und Querprofil. Vermessung des 
FluBlaufes und ihre kartographische Verarbeitung. — V. Wasserstande und Abftiu8mengen oberirdisch flieBender 
Gewasser: Wasserstinde (Pegelstande), AbfluBmengen. — VI, Weitere wasserwirtschaftliche Verfahren: Wasser- 
wirtschaftliche Kennwerte, Leistungsplan eines Flu8kraftwerks (nicht speicherfihiges Laufwerk), Speicherung 
und Wasserwirtschaftsplane dafir. — VII. Hochwasser, Gewdsservereisungen, Wasserstandsnachrichtendienst: 
Hochwasser. Gewdsservereisung. Wasserstandsnachrichtendienst. — VIII. Die Schwerstoffe in den offenen Ge- 
wassern: Herkunft der Schwerstoffe. Der Schweb (Sinkstoff, Schwemmstoff). Das Geschiebe. — IX. Grundwassers 
Entstehung und Vorkommen des Grundwassers. Aufsuchen von Grundwasser, Messung and Berechnung. — 
X. Das Meer im Klstengebiet: Beschaffenheit des Meerwassers..Winde und Wellen. Die Gezeiten (Tiden). Meeres+ 
strémungen. — XI. Qualitative und biologische Gewdsserkunde (Wassergitewirtschaft): Wesen und Aufgabe. 
Qualitative Eigenschaften des Wassers. Biologischer Bereich der Gewassez. — Namen- und Sachverzecichnis. 
* : ; 
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Die Hartzerkleinerung. Maschinen, Theorie und Anwendung in den verschiedenen 
Zweigen der Verfabrenstechnik: Von Carl Mittag, Kéln. Unter Mitarbeit von Dr.-Ing. 
Hellmuth Weinrich, Kéln. Mit 190 Abbildungen. VI, 342 Seiten. 1953. 

Ganzleinen DM 40,50 


Inhaltstibersicht: J. Die Hartzerkleinerung und ihre wichtigsten Maschinen: Einleitung. — Geschicht- 
licher Uberblick. — Die wichtigsten Maschinen der Hartrerkleinerung. — StahlguB und VerschleiBwerkstoffe fir 
Zerkleinerungsmaschinen. — il. Theoretische Grundziige der physixalischen und technischen Zerkleinerung und . 
ihre Bedeutung fiir die Praxis: Einleitung. — Allgemeine Theorie der Hartzerkleinerung. — Theorie einzelner 
Maschinen. — Ausblicke fir die Weiterentwicklung der maschinellen Hartzerkleinerung und neue Entwicklungs- 
wege. — IM. Die Hartzerkleinerung in der Aufbereitungs- und Verfahrenstechnik: Eizleitung, — Hilfsmaschinen, 
— Die Hartzerkleinerung in Anlagen der Aufbereitungs- und Verfahrenstechnik. — Schriittum. — Sachverzeichnis. 


SPRINGER-VERLAG / BERLIN: GOTTINGEN- HEIDELBERG 


* 


Die Metallurgie der Ferrolegierungen. Bearbeitet von Ingenieur F. V. An- 
dreae, Dr.-Ing. W. Dautzenberg, Dipl.-Ing. A. Driller, Professor Dr. 


R. Durrer, Dr. G; Fiore, Dr. W. Freigang, Dipl.-Ing. W. Hilgers, 
Dipl.-Ing. K.-Kintzinger, Whlther Mareth, Dr.-Ing. O. Rosner, 
Professor Dr. E. Schwarzv. Bergkampf, P.Sernitz, Dr. 0. Smetana, 
Dr. phil. G. Volkert, Dipl.-Ing. W. Wilke, Dr.-Ing. J. Wotschke. Heraus- 
gegeben von Professor Dr. R. Durrer, Ziirich und Dr, phil. G. Volkert, Séllingen/Karls- 
ruhe, Mit 188 Abbildungen. VII, 409 Seiten 1953. Ganzleinen DM 72,— 


Inhaltsibersicht: J, Allgemeine Metailurgie der Ferrolegierungen: A. Begriffs- 
bestimmung und Wirkungsweise der Te foe ‘Von Prof. Dr. E. Schwarz v. Bergkampf. — 
B. Geschichtliche Entwicklung. Von Dr. phil. C. Volkert. — G. Die Ausgangsstoffe. Von Dr. phil. GC. Volkert. 
— D. Physikalisch-chemische Grundlagen der Herstellungswerfahren flir Ferrolegierungen. Von ‘Prof. Dr. 
E.Schwarzv. Bergkampf. — E. Die metallurgischen Arbeitsweisen: Das Verhiitten im Hochofen. Von 
Prof. Dr. R. Durrer. Das elektrothermische Verhitten mit Kohle als Red.-Mittel. Von Dr. phil. C. Volkert. 
Das elektrosilicothermische Verhiitten. Von Dr. phil. G. Volkert. Das metallothermische Verhitten ohne 
Zufubr elektr. Energie. Von Dr.-Ing. VW. Dautzenberg. If. Die 6fen: A. Modellbetrachtung des Licht- 
bogenofens. Von Prof. Dr. E.Schwarzv.Bergkampf. — B,. Die elektrische Seite der Ferrolegierungs- 
éfen. Von Dipl.-Ing. A. Driller. — ©. Mechanischer Aufbau der Anlagen. Von Dipl.-Ing. W. Wilke. — 
i. Elektroden: Amorphe Elektroden. Von W. Mareth. Graphitierte Elektroden. Von W. Mareth. 
Siderbergelektroden, Von Dr. W. Freigang. —1V.8pezielle Metallurgie der Ferrolegie-~ 
rungen: A. Ferrochrom: Erschmelzung tiefprozentiger Chromlegierung im Hochofen, Von Dipl.-Ing. 
K. Kintzinger. Erschmelzung von Ferrochrom im Elektroofen. Von Dr. phil. CG. Volkert. Alumino- 
thermische Erschmelzung von Chromlegierungen. Von Dr.-Ing. W. Deutzenberg. — 8B. Ferro-Kobalt. 
Von Dipl.-Ing. W. Hilgers. — GC. Ferromangan und Manganlegierungen: Erschmelzung yon Ferromangan- 
carburé im Hochofen. Von Dipl.-Ing. K. Kintzinger. Erschmelzung von Ferromangan-carburé im Elektro- 
ofen. Von Dr. phil. G. Volkert, Silikothermische Erschmelzung yon Manganlegienpung im Elektroofen. Von 
Dr. phil. G. Volkert. Aluminothermische Erschmelzung von Manganlegierungen. Von Dr.-Ing. W. Dautzen- 
berg. — D. Ferromolybdan; Erschmelzung von Ferromolybdin im Elektroofen. Von Dr. phil. G. Volkert. 
Silikothermische Erschmelzung von Ferromolybdan, Von Dr. phil. G. Volkert, — E. Ferronickel. Von Dr. 
G. Fiore. — F. Ferrophosphor. Von Dr.-Ing. J. Wotschke. — G. Ferrosilizium: Erschmelzung tiefprozentiger 
Ferrosilizium-Legierungen im Hochofen, Von Dipl.-Ing. K. Kintzinger. Erschmelzung von Ferrosilizium- 
Legierungen im Elektroofen. Von Dr. phil. G. Volkert. Caleium-Silizium. Von Dr, phil. G. Volkert. — 
H. Ferrotantal und Ferroniob.. Von Dr.-Ing. O. Résner. — J. Ferrotitan: Erschmelzung der Titanlegierungen 
im Elektroofen..Von Dr. phil. GC. Volkert. Aluminothermische Erschmelzung von Titanlegierungen. Von 
Dr.-Ing. W. Dautzenberg. — K. Ferrovanadium. Von Dr. 0. Smetana. — L, Ferrowolfram. Von Dr. 


phil. G. Volkert..— M. Ferrozirkon-Silizium. Von Dr. phil. G. Volkert. — V.Ausbiick. Von Prof. 
, & 


Dr. R. Durrer. — Sachverzeichnis. “al 


Kurzes Handbuch der Brennstoff- und Feuerungstechnik. vou Dr.-Ing. 
Wilhelm Gumz, VDI. Zweite, verbesserte‘ Auflage. Mit 168 Abbildungen. IX, 
583 Seiten. 1953. i Ganzleinen DM. 45,— 


Inhaltsibersicht: Einleitung, — I. Physikalische Grundgesetze: Kennzeichnung von Festkérpern, 
Schittungen und Stéuben. Mechanik gasférmigers;Kérper. Warmetheorie. Die ‘Warmeibertragung. — II. Die 
Brennstoffe: Zusammensetzung und Eigenschaften der festen Brennstoffe. Die Veredlung fester Brennstoffe auf 
mechanischem, thermischem und chemischem Wege. Zusammersetzung und Eigenschaften der flissigen Brenn- 
stoffe. Gewinnung und Veredlung der fliissigen Brenn- und Treibstoffe. Zusammensetzung und Eigenschaften 
gasférmiger Brennstoffe. Die Veredlung gasférmiger Brennstoffe. — II. Verbrennung: Verbrennungsrechnung fir 
feste und fliissige Brennstoffe. Taupunkt und Schwefelsiurebildung im Rauchgas. Sonderfalle der Verbrennung. 
Ermittlung der Luft- und Rauchgasmengen aus dem Heizwert und der LuftiiberschuSzahl oder dem CO,-Gehalt. 
Die Verbrennungstemperatur. Das It-Diagramm. | Die ‘Warmeverluste. — IY. Vergasung: Vergasungsrechnung 
fiir feste Brennstoffe. Schachtofenberechnung. — V. Verbrennungs- und Vergasungsvorgange: Statik und Dy- 
namik der Verbrennung und Vergasung. Der physikalische Vorgang der Verbrennung. Der Ziimdvorgang. Der 


Verbrennungsvorgang in der Feuerung. Der Vergasungsvorgang im Gaserzeuger. — Anhang. — Namene und 


Sachverzeichnis._ bt 


Vergasung fester Brennstoffe. Stofftbitanz und Gleichgewicht. Eine Dar- 
stellung praktischer Berechnungsverfahren, Von Dr.-Ing. Wilhelm Gumz, Colum- 
bus, Ohio. Mit 4 Abbildungen, VII, 102 Seiten, 1952. Steif geheftet DM 9,60 


Inhaltsibersicht: I. Einleitung. — II, Grundsatzliches zur Berechnung yon Vergasungsvorgingen. — 
III. Bemerkungen zur Kinetik der Vergasungsvorginge. — IV. Das Gleichungssystem und seine Lésungen. — 
V. Verhalten des Schwefels im Gaserzeuger. — VI. Zahlenbeispiele. — Formel- und Tabellenanhang. — Namen- 
und Sachverzeichnis. 
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Diesem Heft liegen je 2 Prospekte des Springer-Verlages, Berlin. Gottingen - Heidelberg und des 
Springer-Verlages, Wien, bei. 
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